Mesoscopic phenomena in multiple light 

scattering 



Academisch Proefschrift 
ter verkrijging van de graad van doctor aan 
de Universiteit van Amsterdam, op gezag van de 
Rector Magnificus Prof. Dr. P. W. M. de Meijer 
in het openbaar te verdedigen in de Aula der 
Universiteit op donderdag 11 mei 1995 te 10.30 uur 

door 



lO Marcus Cornelis Wilhelmus van Rossum 

geboren te Utrecht 

< 



o 
in 



I 

o 
o 



promoter: Prof. Dr. Ir. A. M. M. Pruisken 
co-promotor: Dr. Th. M. Nieuwenhuizen 



commissie: Prof. Dr. B. L. Altshuler 
Prof. Dr. P. F. de Chatel 
Prof. Dr. D. Lenstra 
Dr. I. V. Lerner 
Dr. J. M. Luck 
Prof. Dr. B. Nienhuis 



paranimfen: Kees Diks & Rik Kop 



Van der Waals-Zeeman Laboratorium 
Faculteit der Natuur- en Sterrenkunde 
Valckenierstraat 65-67 
1018 XE Amsterdam 



Dankwoord 



Het boekwerk dat voor u ligt, is niet het resultaat van mijn inspanning alleen, maar was 
slechts mogelijk door de hulp van anderen. AUereerst wil ik daarvoor mijn begeleider 
Theo Nieuwenhuizen bedanken. Bijna al het gepresenteerde werk is het resultaat van 
onze plezierige samenwerking. Zijn snelle inzicht in de fysica en de rest van het leven 
is opmerkehjk. Zijn betrokkenheid bij mijn promotie en verdere loopbaan waardeer 
ik bijzonder. Mijn promotor Prof. Aad Pruisken wil ik bedanken voor zijn inzet en 
behulpzaamheid. Leerzaam was het nauwe contact met de experimentele groep van 
Prof. Ad Lagendijk. Van hem en de andere groepsleden heb ik veel opgestoken. 

Mijn ouders en mijn broer ben ik erg dankbaar voor hun niet aflatende steun en 
voortdurende interesse in de voortgang van dit werk. 

Liefste Karen, ik zou je kunnen bedanken voor het proeflezen van een groot gedeelte 
van de tekst, maar behalve dat misschien een verkeerde indruk zou wekken, staat het in 
geen enkele verhouding tot hetgeen je werkelijk voor me hebt betekent in het afgelopen 
jaar. 

Ook gaat mijn dank uit naar Mevr. Wouters-Hellebrekers voor een coUectie prachtige, 
oude natuurkundeboeken, die behalve historisch interessant, ook nog bruikbaar bleken. 

The hospitality of Prof. Ilya Polishchuk and Dr. Alexander Burin during my visit to 
Moscow are acknowledged with warm feelings. 

Ich danke Michael Schreiber und Etienne Hofstetter herzlich fiir die numerische 
Daten unserer gemeinsamen Arbeit. 

Nee, ik ben juUie niet vergeten: Meint van Albada, Johannes de Boer, Hanna Brands, 
Marco Brugmans, David van Coevorden, Hilde Fleurent, Jan de Gier, Floor Goede- 
mondt, Rogier Groeneveld, Jan Groenewold, Tom Hijmans, Bert Holsbeek, Jiri Hoog- 
land, Wim Koops, Rik Kop, Mark en Ron Kroon, Ad Lagendijk, Jom Luiten, Alec 
Maassen van den Brink, Peter Molenaar, AUard Mosk, Nathalie MuUer, Boris Nieuwen- 
huis, Peter den Outer, Pepijn Pinkse, Irwan Setija, Tijn Smit, Rudolf Sprik, Ruud 
Vlaming, Willem Vos, Jook Walraven, Gerard Wegdam en Diederik Wiersma. Door 
juUie waren de afgelopen vier jaren zo levendig, leerzaam en prettig. 

Een heel bijzondere plaats nemen de AIO's en OIO's van de Spectroscopic der 
Verdichte Materie in. Betere coUega's zijn niet voorstelbaar, noch bij menig mensamaal 
en cappootje, noch bij de inmiddels legendarische AlO-nights. 

De geslaagde omslag en titelblad van dit boek zijn van de hand van grafisch ontwerper 
Karel van Laar, waarvoor ik hem heel dankbaar ben. 

Tot slot wil ik de geinteresseerde lezer graag verwijzen ik naar de Nederlandse samen- 
vatting achterin. 

Mark van Rossum 
Amsterdam, maart 1995 



List of publications 



Works on which this thesis is based: 

• Chapter 4: M. C. W. van Rossum and Th. M. Nieuwenhuizen, Influence of skin 
layers on speckle correlations of light transmitted through disordered media^ Physics 
Letters A 177, 452-458, (1993) 

• Section 2.6 and section 3.2: Th. M. Nieuwenhuizen and M. C. W. van Rossum, 
Role of a single scatterer in a multiple scattering medium^ Physics Letters A 177, 
102-106, (1993) 

• Section 2.6: P. N. den Outer, M. C. W. van Rossum, Th. M. Nieuwenhuizen 
and A. Lagendijk, Locating objects with diffuse light in 'OSA proceedings on Ad- 
vances in optical imaging and photon migration', R. R. Alfano ed., 297-302 (SPIE, 
BeUingham, 1994) 

• Chapter 6: J. F. de Boer, M. C. W. van Rossum, M. P. van Albada, Th. M. 
Nieuwenhuizen and A. Lagendijk, Probability distribution of multiple scattered light 
measured in total transmission^ Physical Review Letters 73, 2567-2571, (1994) 

• Chapter 7: Th. M. Nieuwenhuizen and M. C. W. van Rossum, Intensity distribu- 
tion of waves transmitted through a multiple scattering medium^ Physical Review 
Letters 74, 2674-2678, (1995) 

• Chapter 5: M. C. W. van Rossum, Th. M. Nieuwenhuizen and R. Vlaming, Optical 
conductance fluctuations: diagrammatic analysis and non-universal effects^ to be 
published in Physical Review E (1995) 

• Chapter 6: M. C. W. van Rossum, J. F. de Boer and Th. M. Nieuwenhuizen, Third 
cumulant of the total transmission of diffuse waves, to be published in Physical 
Review E (1995) 

Other works on disordered systems: 

• Th. M. Nieuwenhuizen and M. C. W. van Rossum, Universal fluctuations in a 
simple disordered system, Physics Letters A 160, 461-464, (1991) 

• M. C. W. van Rossum, Th. M. Nieuwenhuizen, E. Hofstetter and M. Schreiber, 
Band tails in a disordered system in 'Photonic band gaps and localization', C. M. 
Soukoulis ed., 509-513, (Plenum Press, New York, 1993) 

• M. C. W. van Rossum, Th. M. Nieuwenhuizen, E. Hofstetter and M. Schreiber, 
Density of states of disordered systems, Physical Review B 49, 13377-13382, (1994) 



Contents 



Dankwoord 


V 


List of publications 


vii 


Introduction 


1 


1.1 


Diffuse light 


1 


1.2 


Light and other waves 


2 


1.3 


Anderson localization 






1.4 


Not quite Anderson localization 


5 


1.5 


Outline of this thesis 


5 


Diffuse Light 


7 


2.1 


Mesoscopics 


7 


2.2 


Amplitude properties 


9 




2.2.1 Wave equations 


9 




2.2.2 Propagators 


10 




2.2.3 The t— matrix 


. 10 




2.2.4 The self-energy 


12 


2.3 


Bulk intensity transport: the diffusion equation 


14 




2.3.1 The incoming and outgoing beams in the diffusion approach . . 


. 19 


2.4 


Radiative transfer 


. 20 




2.4.1 The M 7^ case 


. 25 


2.5 


Transmission quantities 


. 25 


2.6 


Scattering from one extra, static scatterer 


29 


2.7 


Most crossed diagrams 


32 


Interference of diffusons: Hikami vertices 


33 


3.1 


Calculation of the Hikami vertices 


33 




3.1.1 The H4 four point vertex 


33 




3.1.2 Six point vertex: Hq 


35 




3.1.3 Manual for Hikami boxes 


36 


3.2 


Beyond second order Born 


36 


3.3 


Corrections to the conductivity 


38 


3.4 


Appendix 


38 




3.4.1 Beyond the diffusion limit 


. 39 



viii 



CONTENTS ix 

4 Short and long range correlations: Ci, C2 41 

4.1 Physics of the correlation functions 41 

4.1.1 Angular resolved transmission: speckle 41 

4.1.2 Total transmission correlation: C2 42 

4.2 Calculation of the Ci correlation 44 

4.3 Calculation C2 correlation 45 

4.3.1 Focus 47 

4.3.2 Contributions from disconnected diagrams to C2 49 

4.4 Comparison with literature 52 

5 Conductance Fluctuations: C3 55 

5.1 Introduction 55 

5.2 Long range diagrams 58 

5.3 Divergencies in the diffusion approximation 60 

5.4 Cancellation of short distance contributions 64 

5.4.1 Short distance behavior of long range diagrams 65 

5.4.2 Extra short distance contributions 69 

5.4.3 Interference vertices without partner exchange 70 

5.5 Application to optical systems 70 

5.6 Discussion and outlook 75 

6 Third cumulant of the total transmission 79 

6.1 Introduction 79 

6.2 Intermezzo: Cumulants and moments 80 

6.3 Cumulants of the probability distribution 81 

6.3.1 The second cumulant 83 

6.3.2 The third cumulant 83 

6.4 Corrections for the experimental situation 86 

6.4.1 Influence of incoming beam proflle 86 

6.4.2 Influence of internal reflection 87 

6.4.3 Contributions from disconnected diagrams 88 

6.5 Comparison with experiments 91 

6.6 Discussion 91 

7 Full distribution functions 95 

7.1 Classics 95 

7.2 Eigenvalues of the transmission matrix 96 

7.3 Distribution of total transmission 98 

7.3.1 Influence of beam proflle 101 

7.4 Speckle-intensity distribution 103 

7.4.1 Joint distribution 106 

7.5 Concluding remarks 106 

7.6 Higher moments of the conductance 107 

A List of symbols 109 

Nederlandse samenvatting 117 



// y aurait trop d'obscurite, 

si la verite n'avait pas des marques visibles. 



Blaise Pascal, Pensees (Br. 857) 



1 

Introduction 

1.1 Diffuse light 

Light scattering phenomena are seen at many places in nature. It is often the case that 
hght from a source undergoes no or only a few scatterings before it reaches us. One of 
the many examples of such a situation is the sun we see on a clear day. The presence of 
only a few scattering events is also called the "ballistic" regime. 

In this work, however, we will deal with light that has been scattered many times, 
also called diffuse light. This typically happens if the light penetrates a medium with 
many randomly placed scatterers. The intensity then diffuses through the medium. An 
important length scale is the mean free path. It gives the length-scale on which a plane 
wave is converted into diffuse light. Whereas in the ballistic regime the mean free path 
exceeds the system thickness, in the diffuse regime the mean free path is much shorter 
than the sample thickness. As an everyday life example of diffuse light, one can think 
of light in the fog or in a glass of milk, or one can try to look through this very piece of 
paper. 

Although on earth diffuse light is common, astrophysicists initiated the first detailed 
studies in this field. With the starlight seen on earth, they try to reconstruct the 
spectrum and intensity of stars. Therefore, they need to know how the light from the 
stars penetrates the interstellar clouds. Doppler shifts and absorption bands make this 
a complicated problem. Chandrasekhar [1] and Van de Hulst [2] did important work in 
this field. 

Closer to earth there is by now a vast range of applications of diffuse light scattering: 
With remote sensing, one extracts from satellite data the health of vegetation, or one 
can investigate the diffuse reflection of ice and oceans. The reflectance of laser light 
in the sky can be used to study the composition of the air (this technique is called 
"lidar"). There are also medical applications: using diffusely scattered light one tries to 
locate cancer cells or one measures the oxygen content of blood. The harmlessness of 
this technique is of course of great importance, yet extracting the information from the 
scattered light is a complicated problem. 

The typical laboratory systems used to study diffuse light are small slabs. These 
slabs are made of a substance similar to white paint, and are some hundred microns 
thick. The slabs are illuminated by a laser source and the transmission or reflection 
is measured. Very small mean free paths can be obtained in these samples, which is 
advantageous for the study of corrections to the diffusion theory. We will apply our 
theory to slab geometries and compare it to the experimental data from such slabs. 
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Very precise measurements are possible for these systems. 

1.2 Light and other waves 

Many principles of diffuse transport of light also apply to other diffuse wave phenomena. 
Diffuse transport as well as the corrections on diffusion are also studied for neutrons, 
sound waves [3, 4], microwaves and electrons [5]. Scattering of sound or seismic waves is 
studied, for instance, to understand propagation of shock waves caused by earth quakes. 
Oil companies use sound when trying to locate oil. On smaller scales one uses ultrasound 
to detect cracks or other faults in, for instance, metal products. Finally, also phonons 
in solid state fall in the category of diffuse wave scattering. 

Closely related to light scattering is microwave scattering. In that case the wave- 
length is a few centimeters, and a typical sample can consist of a "box" of say, one 
meter in size, filled with metal spheres sized a few centimeters. These length scales are 
of course convenient for experimenting. Presently, the research on microwaves is mainly 
carried out by the group of Genack at Queens College in New York. 

Electronic measurements are typically performed on samples sized a few microns or 
less, so called micro or nanostuctures. Using etching techniques samples with a wide 
range of shapes and properties can be manufactured. One calls the diffuse regime for 
electrons also the metallic or Ohmic regime. 

All systems have their particular advantages and disadvantages. The electronic sys- 
tems, studied very intensively, have the advantage that localization (see below) can be 
easily achieved. An applied magnetic field enters conveniently as an extra tuneable 
parameter. Disadvantages in electron systems are scattering mechanisms that destroy 
the phase coherence. One such mechanism is phonon scattering, which can be reduced, 
however, by cooling the sample to, say, 0.2 K. Also the electron-electron interaction can 
obscure the interference effects we are after. Nevertheless, their effect on localization is 
quite interesting, see, for instance, the book edited by Efros and PoUak [6]. 

Electromagnetic waves, notably microwaves and light, have the advantage that both 
photon-photon interaction and photon-phonon interaction are nearly absent. They al- 
low for "clean" experiments on "dirty" samples. Both angular resolved and angular 
integrated measurements are possible for these systems, resulting in three quite differ- 
ent transmission quantities. This stands in contrast with the electronic systems where 
usually only the total conductance is measurable. A disadvantage of electromagnetic 
waves is that they are not easily localized, as the scattering is often weak. 

Most of the phenomena we are interested in, occur in all systems mentioned. The 
reason for this is that, up to some variations, the same equation, namely the scalar 
wave equation with a random potential, describes the systems. Therefore the study of 
analogies between different fields has received some popularity in recent years. For the 
analogies between optics and electronics see, for instance. Van Haeringen and Lenstra 
[7, 8]. The study of the analogy between the scattering of light and sound dates back 
somewhat more years already [9]. 

The wave character, which is present in all these systems, makes this transport 
different from particle diffuson. The (self)-diffusion of a particle is a subject on its own, 
and is ruled by exclusion, interaction and entropy effects. The diffusion of waves is 
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different as such effects are mostly absent. Instead, interference is the most important 
process. It leads to the diffusion, yet also to corrections on the diffusion theory. 

1.3 Anderson localization 

Anderson pointed out in 1958 that in multiple scattering media the interference processes 
may become so strong that normal diffusion vanishes [10]. It is not simply that the 
diffusion constant tends linearly to zero as scattering increases. In that case localization 
could only be reached in the limit of infinitely strong scattering. Rather, the return 
probability of the intensity becomes higher and higher, reducing the diffusion constant in 
a stronger fashion. The diffusion constant can become zero, which means that the wave 
cannot escape anymore from its original region in space. The intensity then typically 
decays exponentially over one localization length. Yet, in the diffuse regime the wave 
function extends up to infinity. There are thus two different phases, the diffuse, metallic 
regime and the localized, insulating regime. In three dimensions a phase transition 
from the extended to the localized state can occur. In one and two dimensions the 
states are almost always localized. Yet for a finite sample the localization length can be 
much larger than the system size, in that case the states appear to be extended and the 
conductance does not vanish. The localized state is solely the result of the interference 
of the waves scattered by the disorder.^ 

The precise behavior near the transition is of particular interest and still not fully 
understood. The standard diagrammatic technique, which we will use, works well for the 
description of diffusion and the first corrections. It is, however, not suited for the study 
of the transition. Therefore, various techniques were developed to study the behavior 
near the transition and the phase transition itself. We mention the most important 
ones. 

An important step was the scaling theory of localization put forward in 1979 by Abra- 
hams et al. [12]. It states that near the localization the only parameter of importance is 
the dimensionless conductance g (which is the conductance measured in units of e^/h). 
The scaling of (/ as a function of sample size was studied earlier by Thouless [13, 14] 
and Wegner [15]. Abrahams et al. extended those ideas and derived the renormalization 
group equations. They concluded that in one and two dimensions there is no real phase 
transition; the states are always localized. Only in three dimensions a phase transition 
can occur. The states far from the band center are localized, whereas the states near the 
band center with high energy are extended. A change in the Fermi energy, for instance 
by doping, can thus cause the transition from a metal to an insulator. 

It is well known that classical diffusive transport is described by the ladder diagrams. 
The so called most crossed diagrams are the next most important diagrams in the diffuse 
regime, they describe the leading interference terms. By self-consistently summing the 
ladder diagrams and the most crossed ones in the diffuse regime, VoUhardt and Wolfle 

^It would be wrong to think that this is the only scenario for a metal-insulator transition in electron 
systems. In fact, electron systems allow for a whole range of possible transitions between the conducting 
and insulating regimes. Due to the band structure, a maximally filled band can cause insulation, while 
all electron states remain extended. Another mechanism in electron systems is the Mott-Hubbard 
transition, where the electron-electron repulsion hinders the electron hopping, see Mott [11]. 
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found a transition to an insulating phase [16, 17, 18]. The most crossed diagrams 
represent loops of the intensity (as we will discuss below), and the approach can be looked 
upon as a self-consistent one loop resummation. Although a self-consistent approach 
will certainly not include all diagrams of higher order, the method works fine even 
near the transition, because the first few higher order correction terms vanish. Kroha 
preformed numerical calculations for the transition on a lattice using the self-consistent 
technique [19, 20]. Good agreement between the self-consistent theory and numerical 
simulations is found, confirming that the theory contains the most important diagrams, 
even relatively close to the transition. 

It is also possible to perform an exact average over the disorder within a field theoretic 
approach. The resulting non-linear action is not very useful. Yet one can subsequently 
integrate over the fast fiuctuations and consider only the "slow" variables of the system, 
i.e. variables with small spatial variation. This yields the so called non-linear sigma 
model [21, 22, 23, 24], or even a super-string model [25]. With this action it is possible 
to generate systematically all corrections to the diffusion process (in 2 -\- e dimensions), 
thus allowing for an explicit foundation of the scaling theory. The same approach 
can also be used to generate all scattering diagrams contributing to a certain physical 
process. However, for more difficult interference processes also this theory becomes quite 
elaborate, as proving the cancellation of divergent terms has to be done at the level of 
the standard diagrammatic approach. 

A rather recent approach is the random matrix theory. The basic assumption here is 
that the total scattering matrix of the system, as it is very complicated, can be described 
just by putting random elements respecting the symmetries of the system. Surprisingly 
this method works well, and just from the assumption of the random ensemble, it pre- 
dicts many features of the systems correctly and in a simple way. Its is, unfortunately, 
only applicable to quasi one-dimensional situations, and it also misses some delicate 
effects (backscatter cone, the naive value of the universal conductance fluctuations, ...) 
which are present in both experiment and a diagrammatic approach. 

Of course all the theoretical results need to be checked. Yet, specially the study of the 
critical behavior is difficult in experiments: either the system is clean, but the critical 
region is far away (in optical measurements) or the critical regime is easily reached, 
but unwanted interactions spoil the situation (in electron systems). Therefore, one also 
uses numerical data, provided by computer simulations to check the theories, see e.g. 
relevant articles in Ref. [26]. 

Although some five years ago there was hope that the Anderson transition might 
be reached for light scattered on disordered samples, this turned out be based on a 
misconception. It was observed that in time-resolved transmission measurements, the 
average transmission time was very long. This was interpreted as an indication that the 
mean free path was very small. This could mean that the light was close to localization. 
It turned out, however, that the long transmission times arise since the light spends 
much time inside the scatterers. This discovery, made by Van Aldaba, Van Tiggelen, 
Lagendijk and Tip [27], showed that Anderson localization was still out of reach. 

A recent development is the construction of so called photonic bandgap materials. 
In analogy with electron systems, where the lattice parameter is of the order of the 
Fermi-wavelength, one makes a structure with a lattice parameter comparable to the 
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wavelength of the light ( ~ 0.5/im). Like in the electronic case, this introduces a band 
structure. The band-structure for electromagnetic waves is more complicated, as the 
vector nature comes into play. Still there are some structures (although less numerous 
than for scalar waves) that exhibit a bandgap [28]. A small distortion of the gap by an 
impurity could introduce a localized state. This would then be a nice playing ground 
for studying atoms, which can no longer decay at a certain frequency. It probably also 
allows for technical applications such as very efficient light emitting diodes (LEDs)[29]. 

1.4 Not quite Anderson localization 

As the localized state is hard to reach for scattered light, we will concentrate on pre- 
cursors of localization. As stated, the simplest leading process is given by the so called 
ladder diagram or diffuson. It brings the diffuse transport of intensity. In the well 
understood "classical" or "Boltzmann" regime, one only considers these processes. The 
situation becomes more interesting if scattering is stronger, and less probable processes 
start to play a role. In the case of very strong scattering, these interference processes 
will finally lead to localization, yet already in the diffuse regime they can be detected. 
Such effects are also called mesoscopic effects. One can roughly distinguish two types: 

1) The interference of the intensity with itself, thus creating loops. This leads to cor- 
rection effects such as negative magneto-resistance in electronic systems [30, 31], and 
the corrections to the diffusion constant. This is also known as weak localization. In 
the optical case such effects are often too small to be seen: their relative contribution is 
small and, in contrast to the electronic case, the magnetic field couples only very weakly. 
Nevertheless, in refiection measurements these processes can be seen as the enhanced 
backscatter cone. 

2) Another possibility is interference of one intensity or light beam with another one, 
as occurs, for instance, in the correlation of two transmitted beams. In this work we 
concentrate on such processes. They are responsible for the observed correlations and 
the non-trivial distribution functions. 

There are various books on localization and mesoscopics that, though mainly con- 
centrating on electron systems, are also of interest here. We mention the proceedings 
Refs. [5, 26, 32, 33, 34]. Classical waves are discussed particularly in Refs. [35, 29] 

1.5 Outline of this thesis 

In this work we focus on interference effects that occur in the diffuse transport of light. 
In particular, we study the correlation of two or more light beams traversing the sam- 
ple. They can have a slightly different frequency and/or different incoming angle. The 
intensities will turn out to be correlated much stronger than one would expect from 
simple classical diffusion arguments. As a result, the intensities do not propagate fully 
independently, but more collectively, causing large fluctuations and yielding interesting 
distribution functions. These effects become of order unity near the localization tran- 
sition. Although our system is far from localization and most effects are quite weak, 
the experimental detection turns out to be very well possible. We do not focus on the 
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transition to the localized state itself or the associated higher order loop effects, instead 
we perform a diagrammatic expansion in the diffuse regime only including the leading 
processes. 

In the following chapter we study the diffuse transport of light. The basic results 
presented there are long known. We, however, include the precise treatment of internal 
reflection and absorption, important in experiments. We also investigate the possibility 
of locating objects hidden in diffusely scattering media. This is relevant for medical 
applications in order to detect, for instance, tumors. 

In chapter 3 we go beyond the diffusion approach. We introduce the so called 
Hikami vertex, which describes how two diffuse intensities can interfere with each other. 
Whereas diffuse transport corresponds to long-range processes, this interference hap- 
pens on short distance scales, typically one mean free path. The number of interference 
vertices will act as an expansion parameter in the theory. The effects rapidly become 
smaller if the number of interference vertices increases. In the following three chapters 
we do explicit diagrammatic calculations with maximally two of these vertices. In chap- 
ters 4 and 5 we study the correlation between two diffuse beams. As mentioned, due 
to the interfence effects two diffuse beams do not propagate independently through a 
diffusely scattering medium, but have some correlation. Which correlation process is 
leading, depends on the transmission quantity considered. We distinguish three corre- 
lation functions that can be observed in three different measurements: 1) correlations 
in the speckle pattern, 2) correlations in the total transmission, and 3) the correlations 
in the conductance of the sample. The last are for electron systems also known as the 
universal conductance fluctuations. In chapter 6 we extend our calculation of the total 
transmission correlation to the correlation between three beams, or, equivalently, to the 
third cumulant of the total transmission. Although for optical systems this is a weak 
effect, it has been measured recently. In chapter 7 we extend our analysis to arbitrar- 
ily many beams. This allows us to reconstruct the distribution of intenstities on the 
outgoing side for both the angular resolved transmission and the total transmission. 

Many effects calculated and discussed in this thesis were flrst observed in experi- 
ments. Often we flnd remarkable agreement with the experiments. On the one hand, 
this is due to the relatively "clean" experimental situation, and the high quality of the 
experimental data. On the other hand, the close collaboration with the experimental- 
ists at the University of Amsterdam has been very fruitful, as it enabled us to focus on 
experimentally relevant situations and to incorporate the dominant corrections. 



2 

Diifuse Light 

In this chapter we introduce a general setting for basic aspects of diffuse transport of 
hght. We introduce the amphtude Green's function, the averaging over disorder, the 
diffuse intensity and the Schwarzschild-Milne equation or radiative transport equation 
for multiple scattered light. The content of this section is known from diverse works in 
the literature. This is also the chapter where we define our system and introduce most 
approximations . 

2.1 Mesoscopics 

We consider a slab of thickness L, and width W. The slab boundaries are at 2; = 
and z = L (this choice of coordinates refiects the astrophysical history of the problem). 
Most of the time we deal with a three dimensional slab, such that the cross-sectional 
area is A = W X W . Yet we can also consider a quasi one or quasi two-dimensional 
slab (then the width is much smaller than the thickness, yet much larger than the mean 
free path). The slab contains randomly placed scatterers. The scatterers in the slab are 
often static, as is the case for dry white paint samples. It is also possible to suspend 
the scatterers in a liquid. For the short time that the light needs to pass through the 
sample, the scatterers can be considered as fixed at random positions. However, on the 
time scale of, say, milliseconds or seconds, they change their position through motion 
in the fiuid. Integrating the signals over a longer time thus probes different sample 
configurations, and an ensemble average is easily obtained. 

The scatterers used in the experiments are often Ti O2 particles or latex spheres 
(both are a main ingredient of white paint). For simplicity we model the individual 
scatterers by point scatterers, that scatter isotropically, i.e. they scatter equally strong 
in all directions. We consider a low density of these point scatterers and assume that 
their positions are uncorrelated. These are very crude approximations indeed. Firstly, 
the density of scatterers is not low in the experiments, as in order to reduce the mean 
free path, the packing fraction in the experiments is often around 40%. The isotropic 
scattering approximation originates from studies on electronic systems, where scattering 
is dominated by s-wave scattering. However, isotropic scatterers do not exist for vector 
waves and small scatterers are Rayleigh scatterers. Yet this is also not appropriate for 
the experiment, as the diameter of the scatterers is somewhere about the wavelength 
of the light. This choice of diameter again reduces the mean free path, as the scat- 
terers are now resonant and scatter very strongly. Moreover, the shape of the Ti O2 
particles is irregular (potato-like). Clearly the scattering is not isotropic and hard to 
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model at all. Nevertheless, for many aspects the precise individual scattering process is 
not so important. In practice, the effects of realistic scattering are usually coded in an 
experimentally determined transport mean free path that replaces the scattering mean 
free path in the theoretical results for isotropic scatterers. Other theoretical parame- 
ters are left unchanged. Therefore, we do not distinguish between the transport mean 
free path and the scattering mean free path. Indeed, experimental comparison under 
this assumption, when available, is good. Only recently this assumption has also been 
founded theoretically [36]. 

The diffuse character of the light is guaranteed if the thickness is much larger than 
the mean free path £, so that on the average many scattering events occur if the light 
traverses the slab. Such a situation is termed optically thick. If the slab is (optically) 
thin, the system is no longer in the diffuse regime but in the ballistic regime. An 
expansion in the number of scattering events becomes then more appropriate than the 
diffusion approximation that we use here. A hazy air or a piece of smoked glass are 
good examples of systems in the ballistic regime. 

The mean free path in its turn is supposed to be much longer than the wavelength 
of the light. This makes it possible to speak of a light ray between the scatterers. If, 
instead, the mean free path and the wavelength become of comparable scale, localization 
is expected to set in. This is the loffe-Regel criterion, which states that localization 



where A is the wavelength of the light. However, in this work we are far from localization 
and the parameter l/[ki) is small, where we introduced the wave- vector k = 2tt/X. In 
the experiments we typically have l/[ki) ~ 0.01. 

We assume that the spectrum of the light, usually a narrow laser band, is not changed 
by the scattering. (In principle, there is a Doppler shift if the scatterers are moving, 
but this is in our case negligible.) Furthermore, we allow for weak absorption, the decay 
length due to absorption, Labs, is assumed to be much larger than the mean free path. 
We assume that absorbed light really has disappeared from our observation. It can, for 
instance, be emitted at a wavelength that is not detected by the detector or be converted 
into heat. 

The mesoscopic, diffuse regime is in summarized in the inequalities 



For electron systems there is an extra condition that L <^ L-mc^ where L-mc is the in- 
coherence length. This incoherence length is somewhat comparable to the absorption 
length. The phase coherence is lost beyond this length scale, but in contrast to ab- 
sorbed photons, the de-phased electrons still contribute to the transport. If L ^ L-mc^ 
the interesting interference effects will average out. Basically we are now dealing with 
particle diffusion, see Kaveh [39], who also discusses some other interesting variations 
of the inequalities (2.2). 

Furthermore, in practice there is often a different index of refraction inside the sample 
and the outside. This causes diffuse intensity propagating outwards to be reflected at 
the surface, changing the intensity as a function of depth. We do incorporate a refractive 
index mismatch in our approach. 



occurs if [37, 38] 





A < £ < L, Labs- 
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2.2 Amplitude properties 

In this section we introduce the wave equation and its Green's functions, which we will 
need for the description of intensity transport. The Green's function also forms the 
basis for our subsequent diagrammatic expansion in the next chapters. The techniques 
employed are well known and are, for instance, treated in the books of Economou [40], 
Mahan [41] and Abrikosov, Gorkov, and Dzyaloshinski [42]. 

2.2.1 Wave equations 

The Maxwell equations inside a medium are generally given by [43] 

^ ^ _ „ laD 4^. 

1 dB 

V-B = 0, VxE = — . 

c ot 

As the medium is not conducting, the free charge pf and free current j/ vanish. Since 
we only deal with nonmagnetic materials, we take /i = 1 and we find 

-^'E = -^^. (2.3) 

It looks as if the components of E are now decoupled, which would allow for studying 
them independently. Yet they are still coupled, as we still have V • D = 0. Never- 
theless we make a scalar approximation for the electromagnetic field of the light. For 
bulk quantities this is justified since the initial polarization is scrambled after a few 
scattering events. The precise consequence of this approximation on the possibility of 
localization [44] and on skin-layer effects is still unknown. For other applications, such 
as acoustic waves and spinless electrons (that applies to electron scattering for which 
the spin direction is conserved), the scalar property is immediate. 

The time-independent scalar wave equation at a frequency uj reads 

VV(r) + ^e(r)V'(r) = 0, (2.4) 

where c is the vacuum speed of light. We shall consider a slab with dielectric constant 
Co and density n of small spheres with dielectric constant €2 and radius Oq, located at 
random positions r,-. In the limit of point scatterers we obtain 

4 

e(r) = eo + -^a;^(e2 - eo)^^(r-r,), {0 < z < L), (2.5) 

= ei, {z <0; z> L), (2.6) 

where ei is the dielectric constant in the surrounding medium, for instance the glass 
container or air. The vacuum value of the dielectric constant is unity. ^ The wave 



^The average index of refraction is caused by atoms acting like small dipoles. They change the 
dielectric constant but cause no scattering. 
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numbers in the surrounding medium and in the random medium are, respectively, 

ki = — -v/e7, k = mki = W— A;i, (2-7) 
c y ei 

respectively. The relative index of refraction m relates to the dielectric constants as 
m = ^eo/ei. 

Note that the wave equation resembles the Schrodinger equation for electrons in a 
random potential 

- — V2*(r) + [y(r) - E]*(r) = 0. (2.8) 
2m 

Note that in Eq. (2.4) the frequency enters into the potential term. If one considers time 
independent problems the equations and their solutions are similar. This similarity is 
the reason for the analogy between light scattering and electron systems. However, 
differences show up in the time evolutions of both systems. 



2.2.2 Propagators 

In our diagrammatic approach we first need the bare Green's function Gq. It is defined 
as the "inverse" of the wave operator in a medium without scatterers 

i-V' -e)Go{r,r') = 5{r-r'). (2.9) 

For later use in the bulk we have neglected boundary effects, or, equivalently, we set 
Co = ei. As then the scatterer-free problem is translationally invariant, the Green's 
function is easily found in the momentum representation as 

Go(p)= [p'-A;2-^0]"^ (2.10) 

We call this the retarded Green's function; likewise, is called the advanced Green's 
function. The bare Green's function describes the propagation of an amplitude through 
a non-scattering medium. It has a pole at p = k -\- iO. The real-space function only 
depends on |r — r'|, by Fourier transformation^ we find 

f d^p e'P'(''~''') Qtk\r-r'\ 

^"^"■'"■'^ = J (2^)3p2-F-^0 = 4^|r-r'|' ^^"^^^ 

which indeed obeys Eq. (2.9). The 1/r decay is just a geometrical factor describing 
spreading in the 3d space. Apart from this there is no decay. 



2.2.3 The t-matrix 

We now include scattering by introducing first one single scatterer. The propagation 
between two points in space can now occur with an arbitrary number of scattering events 
from this single scatterer. This we depicted in Fig. 2.1. The sum of this series is called 

■^In this work we define tlie Fourier transform as /(r) = J j|^^/(p) exp(ip • r), /(p) = 
/d3p/(r) exp(-ip-r). 
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GO + GOV GO + GO V GOV GO + ... 
+ V + V V + V — V V + 



Figure 2.1: The propagation in the presence of one scatterer. The dressed Green's function 
G contains an arbitrary number of scatterings from this scatterer. 

the dressed Green's function G. In real-space this diagram corresponds to 

G(r,r') = Go(r,r') + |driGo(r,ri)y(ri)Go(ri,r') 

+ |dri |dr2Go(r,ri)y(ri)Go(ri,r2)y(r2)Go(r2,r') + ... (2.12) 

Already in 1908 Mie has derived the exact solution for scattering of vector waves from 
spherical particles of arbitrary size [45], but these results are too involved for our cal- 
culations. We rather model the scatterers with the isotropic point scatterers, thus we 
have the much simpler form V{y) = Sir — Yo)V . For our problem the bare scattering 
strength reads 

V = A^k^al{e2leo-l)l'i. (2.13) 

Thus 

G(r, r') = Go(r, r') + Go(r, v^)VG^{v^, r') + Go(r, v^)VG^{v^, v^)VG^{v^, r') + . . . 
= Go(r,r') + Go(r,ro)tGo(ro,r'), 

where the t— matrix is defined as the sum of all possible scatterings from a single scat- 
terer. Hence ^ 

t = V + VGoV + VG0VG0V + .... (2.14) 

This series is the Born series. For disordered electron systems one often deals with 
weak scattering and the diagrammatic expansion can be carried out in second order 
Born approximation, i.e. t = V -\- VGqV . The real part of the return Green's function 
diverges in three dimensions, since 

Go(r, r') ^ ^ + 1^ + 0(|r - r'|). (2.15) 
47r|r — r I 47r 

As a result, the point scatterer approximation is unphysical. A cut-off length could 
be introduced to solve this problem. As in the second order Born approximation the 
renormalized expression is multiplied by the small factor V^, it then leads to a result 
that is much smaller than the linear term in so the real part of G can just as well 
simply be left out. Yet, the imaginary part should be kept as it describes the scattering, 
resulting in 

+ (2.16) 

47r ^ ' 



^In quantum mechanics textbooks ^pp' denotes the scattering amplitude for scattering of a plane 
wave with momentum p into one with momentum p'. Here we can work in coordinate space represen- 
tation. The effective (re-summed) scattering matrix 5[r)5[r')t is then diagonal, and yields ^pp' = t. 
Though, strictly speaking, t is only the value of the matrix elements, it is commonly again called 
"t-matrix" . 
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In experiments on optical systems, however, efficient scattering is achieved by taking 
resonant scatterers. This means that higher order terms in the Born series cannot be 
neglected and one has to calculate the full Born series. By summing Eq. (2.14) we find 

V 

(2.17) 



yGo(ro,ro) 



for which we could use the same cut-off procedure, i.e. omitting the divergent real part 
of Gq. However, physically this is not a well motivated step. Moreover, it would be nice 
to keep the resonant properties. This can be done by taking the point scatterer limit, 
while keeping the total scattering strength constant. Here we follow Nieuwenhuizen 
et al. [46]. We introduce a cut-off, corresponding physically to the inverse scatterer 
diameter. The regulated return Green's function is denoted Gq and reads 

~ '\P\<i/ao (2^^^^r^^^ " 2^ + 4^- ^^-^^^ 
This yields for the t-matrix 

^ = T^^—YE- (2.19) 

Due to the A;— dependence of V this t— matrix still contains one resonance, occurring if 
the real part in the denominator of Eq. (2.19) vanishes, or A;^ = 37r/[2aQ(e2/eo — 1)]. Note 
that in electronic systems V is A;-independent, and the t— matrix has no resonance. 



2.2.4 The self-energy 

Let us now consider the presence of many identical scatterers placed at random positions. 
The potential now becomes a sum of point scatterers. Connecting scattering events on 
the same scatterer with a dashed line we obtain Fig. 2.2. This brings for the dressed 
Green's function 

G = Go + GqSGo + G^oSGqSGo + • • • 

= (Go ' - = (- _ A;2 - s) . (2.20) 

Before averaging over the locations r,- of the N scatterers we have, to leading order in 
their density, S(r) = J^iLi <^(i" ~ The average over the scatterer positions r,- then 

yields the space-independent result, 

S = nt, (2.21) 

where n = N/V = N/[AL) is the scatterer density. We will use this approximation 
for the self-energy, also known as the independent scatterer approximation, throughout 
this work. We consider no recurrent scattering. In recurrent scattering the amplitude 
returns to a previous scatterer after scattering from another one, such as the rightmost 
self-energy diagram in Fig. 2.2, see for a treatment of such processes Van Tiggelen et al. 
[47]. This process indeed seems improbable in the diffuse regime. Another improvement 
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= + X + —ii X— + —ii X X- + 

+ — ){ 9< — + a )<— 

+ —ii * 

+ X 

+ — * — -d — *- 
E = X + * — k + — — x + ^ X k + 



Figure 2.2: First five lines depict the average Green function (thick line) in a medium 
with many scatterers. The thin lines are the bare Green's functions Gq. We depict the 
scattering potentials by crosses. The dashed lines indicates that the scattering occurs from 
the same scatterer. Next two lines define the self energy. The circles represent the i— matrices 
of the different scatterers. The self-energy S contains irreducible diagrams only. In our 
approximation the self-energy contains only one i— matrix; the rightmost diagram for the 
self-energy, corresponding to recurrent scattering, is thus neglected. The bottom line gives a 
general expression for G in terms of S. 
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on our approximation would be the self-consistent calculation of the t— matrix. This 
leads to a rich behavior in which band-gaps in the density of states can appear [48]'*. 

For our purpose such refinements would complicate the problems too much. The 
problem is that in order to respect conservation laws a more detailed approach on this 
amplitude level also has to be made on intensity level, i.e. for two connected amplitude 
propagators, and also on vertex level (four or more propagators). This would be horrible 
because the number of diagrams then increases dramatically. Thus we stick to our simple 
approximation of the self-energy. Later we will see some examples of errors that occur if 
different approximations are made on different levels. These errors are a manifestation 
of the breaking of conservation laws. 

In the bulk (i.e. the part not under the influence of the boundaries) the average or 
dressed retarded Green's function thus reads 

G(p,p') = G(P)W 

G{p) = (p^-k^-ntY\ (2.22) 



Compared to the bare Green's function, the pole of this dressed Green's function, shifts 
to 

p = VA;2 + nt ^[k + nRe{t)/2k] + i[nlin{t)/2k]. (2.23) 

The real part of the pole shifts a little but this can be neglected, or can be absorbed 
in the deflnition of k. More important is the small but flnite imaginary part. It causes 
an exponential decay. (This essentially is the difference in the optical properties of milk 
versus water.) Indeed, the unscattered (or coherent) intensity is given by the absolute 
value squared of the averaged Green's function. It decays as 

-|r-r'|nlm(t)/fc |r-r'|/4x 

G*{r, r')G{r, r') = — = — -. (2.24) 

loTT^Ir — loTT^Ir — 

We deflned the extinction mean free path as 

The extinction mean free path describes the decay of coherent intensity either through 
scattering or absorption. Indeed, our formalism is capable to describe absorption as well: 
If the dielectric constant €2 of the scatterers is complex, the potential V in Eq.(2.13) has 
an imaginary part, and there will be an extra imaginary contribution to, e.g., Eq. (2.16). 
The leading absorption comes from the linear term in V . The scattering part comes 
from second and higher order terms in the Born series. Eq. (2.24) is a manifestation of 
the Lambert-Beer law. 



2.3 Bulk intensity transport: the diffusion equa- 
tion 



■*The local density of states p(r) relates to the dressed amplitude Green's function as p(r) = 
ImG(r, r)/7r. See Refs. [49, 50, 51] and references therein for detailed calculations of p{r). 
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Figure 2.3: Left: an example of an actual scattering process; a retarded (thick line) and an 
advanced amplitude (thin line) come from the left and share the same scatterers and the same 
path through the sample. Right: schematic representation of the average process, the diffuson. 
We will use this representation in the following chapters. The scattering events are not drawn 
for simplicity. 

The intensity of a wave is the square of the absolute value of the amplitude. Squaring 
the average propagator, one finds the unscattered intensity that decays exponentially 
over one mean free path, as we saw in Eq. (2.24). In the diffuse regime the mean free 
path is small as compared to the system sizes and the unscattered beam gives thus 
only a small contribution to the transport. Examples of the coherent part in diffuse 
scattering in daily life are the weak sun to be distinguished in the fog, or the contours 
of the filament inside an opaque bulb. In the rest of this work we will neglect the small 
coherent contribution to the transmission. 

A different process is responsible for diffuse transport. Intensity transport in trans- 
mission is dominated by ladder diagrams or diffusons, which can only be obtained by 
averaging after taking the product of G with G* . A diffuson is made up by pairing one 
retarded and one advanced propagator that visit the same scatterers in the same chrono- 
logical order. In other words, the diffusion of intensity corresponds to an amplitude and 
its complex conjugate following the same path along the scatterers. We will see that 
in the idealized case no phase difference between the propagators builds up, such that 
averaging does not suppress this process. ^ The diffuson diagram is depicted in Fig. 2.3 
and Fig. 2.4. The diffuson is built of the ladder kernels tied after each other. The kernel 
cont and one common scatterer. We will define the diffuson below to start 

and end on a scatterer. 

The consistent choice of the intensity vertex follows from the intensity conservation. 
Conservation laws are in a field-theoretic approach given by the Ward identities. They 
are the field-theoretic extensions of the Noether currents in classical mechanics. The 
Ward identity for the intensity conservation, as derived by VoUhardt and Wolfie in this 
context [17], says that the irreducible intensity vertices are generated by cutting the 
propagators in the self-energy diagrams and then hipping one part down. This is the 
way we constructed the upper line of Fig. 2.4. It is then clear that taking only a first 
order Born approximation in the self-energy (S = nV)^ would not yield diffusons at 
all (there would be no propagator to cut). In the second order Born approximation 
the scatterers are described by "bare potentials", leading to a vertex factor VV. For 
the full Born series the vertex factor consists of two t— matrices, i.e., the product tt, as 

® Similarity to meson physics, where quark and anti-quark are paired together on long distances, was 
noted by Prof. Dr. J. Smit. 
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illustrated in Fig. 2.4. 

Although the scattering mean free path is incorporated in the diffuson equation, 
it is useful to define it on forehand. In second order Born approximation it equals 
£ = 47r/(n|yp). For the full Born series we define it as 



f = ±L. 

ntt 



(2.26) 



The scattering mean free path £, in contrast to the extinction mean free path £ex, comes 
solely from scattering and not from absorption. The two relate as £ex = where a is 
the albedo of the single scatterer. With the definition of the t— matrix Eq. (2.19) we see 
that 



f Im(l/V) 



-I -1 



(2.27) 



If the scatterer-potential is real, there is no absorption, corresponding to a = 1, this is 
also called conservative scattering for obvious reasons. The case a = means that the 
scatterer absorbs all light that hits it. Equivalently we can derive the relation 



ktt 
4tt 



alm(t). 



(2.28) 



For a = 1, which corresponds to no absorption, this is also known as the optical theorem. 
It says that intensity lost in the decay of the coherent beam is seen back in the scattered 
intensity. If a < 1 it is partly absorbed. 

The diffuson is the sum of processes with one, two, three, etc., common scattering 
events. We now calculate the building block of such processes: the ladder kernel. It 
describes one single common scattering event. In this section we assume a thick slab 
and use the diffusion approximation for the boundary behavior. In Fig. 2.4 we show 
part of the diffuson diagram. The ladder kernel U corresponds to the expression 



U = ntt 



(2.29) 



Here we included the momentum difference q and the frequency difference Alo. After 
Fourier transformation they give the spatial and the time dependence of the diffuson. 
In the diffusion approximation variations of the intensity are small. This corresponds 
to expanding the Green's in small momenta and frequency difference 



G{p + q;u + -Au) 



o „ . A; 

c 



4, = /75^G^p)G"(p). 



p + 2p • q + g 

^ G+(-2p-q-q^ + Auj-)G^ + 4(p-qfG^ 

c 

with G = G(p;c<;). We define the "standard" integrals 

d^p 



(2.30) 



(2.311 
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p+q ; co+Aco/2 p'+q p"+q 




Figure 2.4: Detailed picture of a part of the diffuson. On the top line we have drawn one 
possible scattering contribution to the ladder kernel. We depicted the scattering potential with 
'x'. By its definition the full ^— matrix, the circle, contains all such repeated scattering events. 
The Green's functions are bare in the upper line. The bottom line is part of the product of 
ladder kernels. The sum of diagrams with an arbitrary number of kernels leads to the full 
ladder or diffuson. Note that the Green's function are dressed now. 



We calculate these integrals in section 3.4; their pole at p ^ k dominates them. We 
expand in second order of q. Due to symmetry of the integral, terms linear in q cancel, 
and we are left with 

U = ntt(^In + 2Auj^l2i-q'l2i + ^k'q'hi^ (2.32) 

= (l-^£^(g^ + «^-^|Acu)). (2.33) 

The g^/21-term was neglected as it is proportional to l/[ki). The absorption length 
Labs gives the decay length of the intensity due to absorption. However, the inverse 
absorption length k = L~^^ is more convenient for us. It relates to the albedo a as 

KHy3 = l-a, (2.34) 

when absorption is small, i.e. /«£ ^ 1, or Gi ~ 1. The geometric sum of the kernel, the 
ladder-sum, gives the full diffuson 

°° ■ nff 

£(q) = ntt^t/^ = ^— -, (2.35) 

j=0 

or, in the approximation of Eq. (2.33), 

[g2 + «2 + zO]£(q) = ^. (2.36) 

Thus C obeys the well known diffusion equation^ written here in momentum space. Note 
that we have multiplied Eq. (2.35) with one extra factor nti to let the ladder end on a 
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scatterer (conform Fig. 2.4). In taking the lowest value j = 0, we also include the single 
scattering event. We have defined the reduced frequency difference as 

n = (2.37) 

with the diffusion constant D = ci/3. 

It was shown that D = ^ve£ instead of |c£, where ve is transport speed [27, 52, 46]. 
This energy transport velocity is not derived here. It is the result of the resonant 
scatterers, which makes the t— matrix strongly frequency dependent. If now G and G* 
have a somewhat different frequency, they "see" a different t— matrix. This results in 
a much lower diffusion constant in dynamical problems than expected from D = ci/3. 
The physical picture is that the light gets trapped inside and close to the resonating 
scatterer. The discovery of this effect limited the excitement that, concluding from 
the diffusion constant, the mean free path of light might be close to the wavelength in 
experiments. According to the loffe-Regel criterion, Eq. (2.1), localization should then 
be within reach. However, the small diffusion constant was the consequence of ve^ which 
can be much smaller than c , and the actual mean free path was still large. 

Let us return to the diffusion equation 2.36. Let us by C-mt denote the corresponding 
diffusion propagator from r to r'. In an infinite medium it obeys 

[_V2 + ^ + zO]£i„t(r) = ^Sir - r'). (2.38) 

In the slab geometry we can perform a Fourier-transformation in the x, y-directions, and 
the diffusons then obey 

- ^A„t(^, z'; M) + M^CUz. z'- M) = ^5{z - z'). (2.39) 
Here we have defined the decay rate 

= g2 + + zO, (2.40) 

where Q is the 2D transversal momentum^. The complex parameter M describes the 
exponential decay C ~ exp(— Re(M)2;) of the diffuse intensity in the 2;-direction. This 
decay can occur either as a result of absorption, or as a result of de-phasing caused by a 
frequency difference or a transverse momentum difference. In analogy with the range of 
the interaction particles of the fundamental forces, we sometimes call this the mass. In 
our previous derivation we assumed that |M|£ ^ 1, also below we will work only in first 
order of ML For most calculations we will first consider the M = case, to estimate 
the magnitude of the effect. 

By realizing that the diffusion equation (2.39) is just a wave equation with complex 
frequency, one sees that the solutions to the diffusion equation are a linear combination 
of hyperbolic sines and cosines. The solution of the diffusion equation in a slab geometry 
can also be obtained using the method of "image charges" as known from electrostatic 
problems. We then fulfill the boundary conditions by summing the original and mirror 



'We denote 2D vector with capitals, 3D vectors with small bold faces 
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solutions of the infinite slab. Here, however, that method is more involved. Introducing 
a yet undetermined "extrapolation length" Zq^ we can write [53, 54] 

12^ [sinh + M^o cosh Mz] [sinh M(i:-/) + M^o cosh M(i:-/)] 
(M + M^z^o) sinh ML + 2M^zo cosh ML 

(2.41) 

where we assumed that z < z\ otherwise z and z' must be interchanged on the r.h.s. 
The diffuson C describes the diffuse propagation from one point in the slab to another. 
With roughly equal indices of refraction inside and outside the sample, the extrapo- 
lation length Zq is a couple of mean free paths and thus the terms involving Zq yield 
contributions of the order i/L. For optically thick samples [L ^ i)^ this is negligible 
and one has 

^ , , , 127rsinhMz smh M(L - z') 

Consistent with this expression the boundary conditions become in that case £int(0, z') = 
0, Cint{z,L) = 0. 

The case of amplifying media is also of some practical interest. The diffuse intensity 
then increases at each scattering (this can be done by pumping the sample with light 
of another frequency). The albedo is now larger than one, or, equivalently, k becomes 
imaginary. If Q = and = 0, the diffuse intensity has a negative curvature and 
consists of sines and cosines instead of hyperbolic functions. On the top of the enhanced 
backscatter cone (see below) there appears an extra sharp peak [55, 56]. 



2.3.1 The incoming and outgoing beams in the diffusion ap- 
proach 

The diffusion equation does not hold if the intensity gradient is steep, i.e. g£ ~ 1. It can 
therefore not describe properly how the intensity couples into and out of the medium. 
Nevertheless, in a heuristic manner one makes often the following assumptions in the 
diffusion approach, see Ishimaru [57]: 

• The diffuse intensity from an outside plane wave source is assumed to be given 
by substituting z = i m Eq. (2.41). As if all diffuse intensity originates from this 
z = i plane. Likewise, the coupling outwards is obtained by taking z' = L — L In 
this way two extra propagators, the incoming: 2; < 0, < 2;' < L; and the outgoing 
Q<z<L,z'>L are constructed from the internal diffuson (0 < z^z' < L). 

• The boundary conditions are [54] 

= zo. (2.43) 

«=0+,L- 

Solution (2.41) fulfills these conditions. 

The extrapolation length Zq is determined by the refiectance at the surface. In the 
next section we show how we do this properly. Here we treat the diffusion approach. 
With non-reflecting boundaries the boundary condition on the diffuse intensity on the 



C{z,z'-M) 
dX{z,z';M) 



20 



DIFFUSE LIGHT 



surface is the absence of diffuse intensity propagating into the medium. Physically this 
is immediately clear, as there can be no diffuse source outside the medium. In the index 
matched case Zq is 2£/3. This well known result is, for instance, derived in the book of 
Ishimaru [57]. 

However, in realistic systems there can be internal reflection. One of the flrst treat- 
ments of internal reflection in the diffusion approach was given by Lagendijk, Vreeker 
and De Vries [58]. Here we briefly summarize the result of Zhu, Pine and Weitz [54]. 
The reflection of a scalar wave on an index mismatched boundary is [59] (deflne /i as 
the cosine of the angle of the wave to the 2;— axis) 



//i- J/i2-l + l/m2\ 
m = V ' (2.44) 



where m = rio/rii is again the ratio of the refractive index of the sample, no, and it's 
surrounding, rii. Considering the current of diffuse intensity at the boundary of the 
slab, one flnds [54] 

zo = 2£l±^ (2.45) 
3 1 - 2i?i ^ ^ 

where the angular averaged reflectance for the intensity and current, respectively, are 

Ri= f d/i fiRijj) i?2 = / d/i ii^R{ji). (2.46) 
Jo Jo 

If there is no internal reflection, one has m = 1 and R{n) = 0, now Ri and R2 vanish 
and one retrieves Zq = 2/3£. For m ^ 1 the extrapolation length increases. If the index 
mismatch gets larger and larger, the gradient of the intensity near the surface will be less 
steep. Then the diffusion approach, which is good for the bulk, works also flne near the 
surface. One can also say: if the boundary strongly reflects, it is invisible for the diffuse 
light from inside the slab, since the diffuse intensity inside the slab is everywhere the 
same. Therefore, the bulk behavior applies everywhere. The correctness of the diffusion 
approach for this case was conflrmed on an analytical basis in Ref. [59]. 



2.4 Radiative transfer 

Although the diffusion equation describes the transport deep in the bulk accurately, it 
breaks down near the surface as the intensity is not longer slowly varying as function 
of z. Nor is the angular average representative, as the incoming waves have undergone 
only few scatterings. The precise behavior near the (reflecting) surface has to be derived 
from the Schwarzschild-Milne integral-equation. This equation has been known in the 
literature for a long time, especially in the astrophysics community [2, 1]. Particularly 
in the book of Van de Hulst the diffuse intensity is studied in this approach for very 
many situations, including absorption, the i ^ L regime, and anisotropic scattering. 
Yet in those works the average indices of refraction of the scattering medium and its 
surrounding are assumed to match. This is fully justifled for interstellar and meteoro- 
logical clouds, which are very dilute and hardly refract. (The change in the imaginary 
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part of the dielectric constant in going from one medium to another is important for 
describing the scattering, but in those systems the change in the real part is negligible.) 
Nieuwenhuizen and Luck [59] extended the transport equation to the case where there 
is a mismatch between the index of refraction in the scattering medium and outside, 
causing internal reflection. We flrst consider the M = case, see Eq. (2.40). 

For the moment we take our slab semi-inflnite, flUing the space z > 0. Consider a 
plane wave with unit flux impinging on the sample under angle 6a 

tjj^Jr) = ^ ^^QaP+^k^ cosB^z ^ ^ ^ ^2.47) 

\/ Aki cos 8a 

where p = (x,y) is the transversal coordinate and Qa = A;i sin ^a(cos sin (/)a) is the 
two-dimensional transverse momentum of the incoming beam. At the z = 0— plane 
part of this beam is transmitted, and part is specularly reflected. Requiring a smooth 
behavior of if) gives 



V'- 



1 2p, 



iQap+iki cosOaZ Pg Pa ^iQap-iki COsOaZ 



{z<0) 



Pa+Pa 

} ap+ik cos e'aZ-Z I (21 cos e'a) ^ y Q^^ (2.48) 



^/Ak^~COs9^ Pa + Pa 

where we deflned 

pl = kl-Ql Pl = e + nt-Ql. (2.49) 

The angle 6'^ is the angle of the refracted beam with respect to the 2;-axis. 

Inside the slab the unscattered part of the intensity decays over one mean free path, 
see Eq. (2.24). The source of diffuse intensity is the flrst scattered intensity. The flrst 
scattered intensity is the intensity "just after" it has hit the flrst scatterer. We choose 
this deflnition because the next scattering event can now also easily be chosen to end 
on a (different) scatterer. This is consistent with our previous derivation of the diffuse 
intensity, see Fig. 2.4. In leading order of (l/ki) the source S equals 



Saiz) = nttmr)\ 



2 



^-z/(i,a) ^ 4^ j.^^^^ ^-z/ii,a)^ (2.50) 



AkilJai{Pa + Paf Akjla^ 

where /i^ = cos 6'^. The source is proportional to the intensity transmission coefflcient 
of the boundary between the two dielectrics 



{Pa + Paf [p + Jll^-l + l/m2]2 



which obeys r(/i) + R{n) = 1. If /i^ < 1 — we see that = 1, corresponding to 

the case of total internal reflection. 

From Sa the twice scatterered intensity is generated, ntt /dr'|G(r, r') pS'a(r'). Like- 
wise, we obtain the intensity scattered three times, etc. The geometric sum of flrst, 
second, third,..., times scattered intensity is again the multiple scattered or diffuse in- 
tensity This is exactly the same sum, or the same diagrams, as we used in the 
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derivation of the diffuson equation, but the assumption of slowly varying functions is 
not necessary and we treated the coupling into the medium properly. The diffuse in- 
tensity now is calculated from the first scattered intensity by the Schwarzschild-Milne 
equation 

C'iz) = Sa{z) + dz'Msuiz, z')C'{z'), (2.52) 

which is a self-consistent transport equation. By inserting £ = at the r.h.s. and 
iterating the solution, one recovers the original ladder sum. 

The kernel Msm describes the transport between two scatterings 

Msm(^, z) = Mb{z- z') + Ml{z + z') + Ml{2L -z- (2.53) 

that decomposes into the bulk term Mb and into layer terms describing internal refiec- 
tions at the interfaces at 2; = and z = L. The bulk term is given by integrating an 
intensity in the p' = (x', y') plane at z' propagating unscattered towards a point in the 
plane at z. By putting p' and z' temporarily zero one obtains. 



Mb(z,0) = 47r jdxdy\G{z,p;0, 



dxdy — — = / —e-'i'i", (2.54) 

^ Att{x^ + y^ + z^) io 2/i ' ^ ' 

Thus 

Mb{z - z') = £ l^e-l-'l/^- = \e,{\z - .'I), (2.55) 

with El being the exponential integral. Notice again the exponential decay of the 
intensity between scatterings. 

Similarly, we derive the layer kernel. It is proportional to the refiection coefficient 
R[p). Now the argument in the exponent is the total path length from z' to the surface 
at 2; = 0, and back again to the final point z. The kernel therefore reads 

Ml{z + z') = r ^[1 - r(/i)]e-(^+^')/^^. (2.56) 
Jo Zp 

Here we assumed specular refiection at the surface: pin = Pout- 

The precise solution of Eq. (2.52) can be obtained numerically (see Refs. [2, 60] for 
details). We have plotted it in Fig. 2.5 for a relatively thin slab [L = 4£). We took the 
data from table 17 in Ref. [2]. The albedo was unity, and there was no index mismatch. 
The form of the solution near the incoming surface is quite different for the three cases 
drawn: /i = 1, or perpendicular incidence; p = 0.1, the angle with the z— axis is then 
84 degrees; and diffuse incidence uniformly distributed over all angles. At the outgoing 
surface all solution behave the same: there is only a small deviation from the straight 
line crossing zero at L -\- Zq. 

Generally we can split the solution of Eq. (2.52) in a special and a homogenous 
solution, Cs and Ch respectively, which correspond to the solution with and without 



^For non-specular reflection one has in general: / Jd/i' R[iJ, ij')e ^1^^ ^ 1^ ^. One recovers spec- 
ular reflection by putting R{iJ,, p,') = S{p — p')R{p). Lambert reflection corresponds to R{p, p') cx p. 
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Figure 2.5: The solution of tlie Scliwarzscliild Milne equation as taken from the table 17 in 
Van de Hulst [2] for slab of thickness L = 4i with no internal reflection, a = 1. The precise 
solution near the incoming plane strongly depends on the cosine of the incoming angle, ^u. The 
bulk behavior is apart from a multiplication factor the same, they all extrapolate to L -\- zq 
(thin dashed line). 

the source S", respectively. For thick slabs the intensity propagator in the region, say 
< z < L — can thus be written as 

C''{z) = aCH{z) + Q{z). (2.57) 

Apart from a linear component there are exponential corrections in these solutions, as 
is clear from Fig. 2.5. Far from the surface, these corrections of order e~^/^ vanish. We 
can impose the asymptotic behavior jCs{z) — > const. ^ jCh{z) ^ {z -\- Zq) and the solution 
behaves for, say, 5i < z < L — 5i , as 

47rr(/ia)ri(/ia) / , X 
^ (^) = kiAfi, - "(^ + ^o) • (2-58) 

On the other hand, the propagator on the outgoing region, say 5i < z < L — obeys 

C{z, M) = (3Ch{L - z,M) = (3{L- z + zo) . (2.59) 

Matching of the two solutions in the middle region of the slab determines a and (3. In 
the bulk of the slab (at least few mean free paths away from the boundaries) the diffuse 
intensity thus h simple behavior 

4TTT{^a)Ti{^a) L + Zp - Z 

Note that in this region the intensity again satisfies the diffusion equation V^£?j = 0. 
Indeed, using Taylor expansion of £, it can be shown that the Schwarzschild-Milne 
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equation reduces to the diffusion equation if the diffuse intensity is slowly varying with 
respect to the mean free path, i.e. deep inside a thick slab. We only use the more 
complicated Schwarzschild-Milne equation to fix the parameters of its solution, i.e. Ti 
en Zq. In principle one could use the precise form of solution resulting from the Milne 
equation, but this would be much too complicated for further calculations. For thick 
slabs the contributions from the boundary corrections are exponentially small anyhow. 

The Ti in Eq. (2.60) is the limit intensity of a semi-infinite slab [59]; of course, it 
depends on the incident angle. For matched refractive indices (m=l) and perpendicular 
incidence the numerical value is ri(l) = 5.04 [59]. In the calculations of the correlation 
functions in the next chapters, this pre-factor will drop out as we always normalize to the 
average. The extrapolation length Zq in Eq. (2.60) is more important as it infiuences the 
functional form of the diffuse intensity. For isotropic scattering and matched refractive 
indices one has [2] Zq = 0.7104£. However, it becomes larger when there is a mismatch 
of the indices of refraction between the scattering medium and the surroundings. The 
parameters Zq and Zi [zi is defined as Zi = Tii) are calculated numerically in Ref. [59] 
for various values of the index mismatch m. For given m both Zq and Zi take values 
comparable to the mean free path. Note that the diffusion approximation does quite 
well in the estimate of Zq = (2/3)£ as compared to the more detailed Schwarzschild- 
Milne approach that yielded Zq = 0.71£. In the limit of large index mismatch, m — > 0, 
or m — > oo, the extrapolation length tends to infinity. The estimate of Zq using the 
diffusion approximation of Zhu et al. [54] becomes then asymptotically exact [59]. The 
asymptotic value is already quite good at moderate index mismatch, 3/2. 

However, the index of refraction of a diffuse medium is experimentally and theo- 
retically not easily determined. The sample typically is a dense packing, say 30%, of 
irregular shaped particles suspended in a liquid or solid. The average index of refrac- 
tion of a binary mixture is a classical problem [61, 62, 63]. Den Outer and Lagendijk 
[64] have this studied problem both experimentally and theoretically by considering the 
maximum width of the enhanced backscatter cone, see section 2.7. It turns out that for 
a typical system the average index of refraction is best described by independent Mie 
scatterers. It also possible to measure the refractive index by Brewster angle measure- 
ments, although this method is only sensitive to a very thin surface (few wavelengths) 
of the sample, which may be not representative for the rest of the sample [65]. Both 
the treatment in the diffusion approximation and the Schwarzschild-Milne equation as- 
sumed specular reflecting boundaries. Due to the roughness of the surface, however, 
this is experimentally not fully justifled. The surface of typical samples scatters over 
some 40 degrees [66]. The opposite limit, diffuse surface scattering, is called Lambert 
reflection. It is studied for a full reflecting boundary in the book of Kagiwada et al. [60]. 
In Ref. [2] also the thickness dependence of Zq was studied, yet for already moderate 
values of Ljl ~ 4, Zq is close to its value for inflnitely thick slabs. Concluding from 
the above remarks, we consider the extrapolation length Zq as an extra experimental 
parameter. For small Zoj L we have 





2.5. TRANSMISSION QUANTITIES 



25 



2.4.1 The M ^ case 

Similar to the M = case, one can consider the M-dependent solution of the Schwarz- 
schild-Milne equation. For our optically thick samples [L ^ i) somewhat away from 
the surface, we have homogeneous and special solutions of equation (2.52) 

^^(^'^) = m + £ ' 

^^^"'^^ = — kij;!: — ' ' ^2-^2^ 

where we used the analogy with the diffusion equation for the bulk behavior. The 
extrapolation length Zq depends on M. Van de Hulst and Stark treated the dependence 
of the extrapolation length on the albedo [67] (where a different definition of Zq is used). 
Again the Mi <^ 1 simplifies the situation here. In a way similar to the treatment of 
large index mismatch in reference [59], it can be shown that in the limit of large index 
mismatch Zq i ^ ^/M^ as m — > oo, where Mi ^ 1. A useful interpolation formula is 
[68] 

.oi(M;m) = ^Ji%^. (2.63) 

We denote 2;o,i(0;m) as Zq^i. We have checked this formula for M = n numerically to 
be correct within a few percent up to Mi = 1/2. 

Using the same matching procedure as before, and assuming M— independence of 
T[fj,)^ we find the incoming diffuse intensity needed for our further calculations 

.a( . _ 4^r(/i,)ri(/i,) sinh(Mi: - Mz) + Mzp coshjML - Mz) 

' kiAfia [1 + M^zi)smh{ML) + 2Mzocosh{ML) ' ^ ' ^ 

This formula is similar to the diffuse intensity used in the diffusion approximation above, 
but the parameters Zq and Zi are now precisely prescribed. 

2.5 Transmission quantities 

Untill now we considered an incoming wave from a fixed direction. In optical systems 
there is also the possibility to sum over incoming or outgoing directions. We distinguish 
three quantities 

• The angular resolved transmission Tab 

• The total transmission = J2i Tab 

• The conductance T = J2ab Tab 

The calculation of the auto-correlations in these transmission quantities will be a large 
part of this study. In this section we just calculate the average values. The used 
formalism is also known as the Landauer approach, after the Landauer formula for the 
conductance, [69, 70] 

G='^-^Y.Tab. (2.65) 

ah 
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where one sums the transmission coefficients Tab of waves with unit flux coming in chan- 
nel a and going out to channel b. (In the optical case "one channel" means that one 
incoming or outgoing direction is probed.) In optics the Landauer approach is very 
natural, as the basic quantity is the angular resolved transmission Tat, yet a summation 
over the channels is possible. In optical experiments one carries out this summation 
by placing an integrating sphere, which is a sphere whose inside is coated with white 
paint. It has two small holes, one for placing the sample, and one for the detector. Us- 
ing this technique total transmission measurements were performed by integrating over 
all outgoing directions [71, 72, 73]. Conductance measurements are principle possible 
with integrating spheres on both sides of the sample. (We don't know of published 
conductance measurement in optical systems, neither by using integrating spheres, nor 
by other methods, see Ref. [74] for experimental considerations.) 

In the previous sections we already derived the diffuse intensity inside the slab rele- 
vant for the angular resolved and total transmission measurements, when one incoming 
angle is used. Conductance measurements involve waves coming in from all directions. 
Summing the source term of the previous section, Eq. (2.50) over the channels a yields 
a source for the diffuse intensity 



~~ Jo ZTT 



- d^T{^)e-^l'r (2.66) 



This is again the input in the Schwarzschild-Milne equation (2.52). Deflning 

^ _ 7rr(/ia)ri(/ia) 



(2.67) 



and by comparing the pre-factors, the intensity in the bulk now has the diffusive behavior 
A„(^) = [Ea %] ^fv^- One thus has 



We can use the sum rule 



£r„(z) = eXUz)- (2.68) 



Here we used that ji = ^Jl — Q'^/k'^. The last equation is a sum rule on the diffuse 
intensity, see Eq.(2.30) of Ref. [59]. We flnd 

^ , , 4:k L -\- zq — z , , 

We term this object the incoming total- flux diffuson. Although it has a different prefac- 
tor, it has the same depth dependence as This can be understood as the dependence 
of the solution of the Schwarzschild-Milne equation on the incoming angle only enters 
in terms that describe low order scattering. These are only relevant near the surface. 
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see Fig. 2.5. Deep inside the bulk, these exponential corrections vanish and the only 
angular dependence is in the pre-factor. 

Having calculated the coupling into the sample, we also want the radiation coming 
of the slab on the outgoing side. On the outgoing side, radiation emitted at a point 
r = (p^z) inside the slab propagates to a point {Z^p') outside the sample [Z > L) 
as described by the Green's function of semi-infinite medium with dielectric function 
e(r) = Co for z < L and e(r) = ei for z > L 

G{p, z; p', Z) = - Y. G{z: Z- g)e'«(^-^'), (2.71) 
^ Q 

in which Gi^z] Z] Q) is the 1-d Fourier transform of Eq. (2.22), 

G{Z-Z-Q) = ^_^^P(L-.)^^P(Z-L) ^ 

P + p 

P = ^Jk^-Q^ + nt, p = ^Jkl - Q\ (2.72) 
In the far field [Z ^ L) the total transmitted intensity reads 

/ dV|G(p, z- p\ Zt = \Y. \G{z-^ (2-73) 

■' ^ Q 

Since pz = ki cos 6a = kpa^ the according fiux is 

m = -,Y.kl^\G{z-Z-Q)\' = A ['dpT{p)e-^^-y'^. (2.74) 

A Q OTT Jo 

So the outgoing intensity is generated near the outgoing surface z = as is obvious. It 
leads to a source S{z) = 4:TT^[z)/i in the Schwarzschild-Milne equation. Apart from a 
refiection in z^ this expression differs by a factor 4 from (2.66). The pre-factors would 
be the same if our Green's functions were multiplied by a factor 2; this amounts to the 
same as taking a kinetic term V^/2 rather than V^, such as occurs in electronics in 
units where h = m = 1. The corresponding diffuson is the outgoing total-flux diffuson 
in transmission. It reads 

/ N k z -\- Zc) , _ , 

If z^jL is small we may use 

HouAz) = J J. (2.76) 

The outgoing diffuson gives the diffuse intensity inside the slab that will flnally come 
out at the outgoing side. 

For angular resolved measurements the outgoing intensity in a certain direction b is, 
in analogy with (2.68), 

Cl,{z) = efc£out(^) (2.77) 

We call the incoming and outgoing diffusons, external diffusons, because they connect 
to the outside of sample. This in contrast to the internal diffusons that begin and end 
inside the medium (for instance at an interference vertex). 
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We find the conductance (T) by integrating the outgoing diffuson with source of the 
incoming intensity, Eq. (2.66). Equivalently, one can integrate the incoming diffuson 
with the outgoing source, Eq. (2.74). For this calculation the precise behavior of the 
diffuson near the boundary is important and Eq. (2.75) cannot be used. Instead, we use 
that ^ 

l^dzCout{z) ldiiT{ii)e-'/^^' = 2k/3i. 

The derivation of this sum-rule falls outside the scope of this chapter, see however 
Eq. (2.19) of Ref. [59]. We obtain 

g = {T) = Jd'rCU.z)S{z) 

PAi mi 

2.78 



3^(1 + 2zq) 3(1 + 2zq) 



For the electronic case one has for the average conductance [G) = ^ 3n(L+2zo ) " 
Eq. (2.78) the number of modes is = A;^A/47r, which is also known as the Weyl 
formula. It corresponds to dividing the total beam area into small coherence regions, 
or channels, of area A^. The conductance can then be interpreted as the sum of N 
channels, on the average a channel contributes 4:i/3{L -\-2zo) to the conductance. When 
the diameter of the incoming beam profile becomes comparable to the thickness, the 
diffuse intensity broadens inside the slab. The number of modes becomes then ill- 
defined. (We study such effects in chapter 4.) The number of modes is also one of the 
few parameters where the vector character comes in: The two polarization directions 
travel independently through the slab, so the number of available independent channels 
effectively doubles. In our approach only in 1/A^-correction this effect is important, else 
it just changes the average transmission and does not change our (normalized) results. 

The angular and total transmission coefficients read 

Transmission in the diffusion approximation is proportional to the value of the dif- 
fuson at 2; = £ and z' = L — L The internal diffuson always ends on scatterers. Yet 
for transmission quantities, coupled to a source/detector outside, the diagram should 
end on propagators. As it is not clear how to obtain this in the diffusion approach, on 
often estimates the transmission diagram ending with (four) propagators by dividing by 
[ntiy. The total transmission is thus in the M = case 

(T \ 1^^^ ~" - 3£ 

^ " £3 ^ + 2^0 W y ~4tt{L + 2zo) ^ ^ 

One observes that the dependence on i and L is correct, but the pre-factor is not correct 
and all angular dependence is lost. Note that the region where the diffuse intensity is 
non-zero, seems to exceed the slab on both sides by a length Zq^ see Eq. (2.41), Eq. (2.64). 
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This is seen back in the transmissions Eq. (2.81), but also in Eq. (2.80). The resistance 
of the sample is inversely proportional to the conductance, thus R oc L -\- 2zo. The 
first term is the thickness dependence known from Ohm's law, the second term can be 
interpreted as contact resistances at both sides. 

For the total transmission and the conductance measurements one integrates over 
the outgoing surface. If on the incoming side, two or more different plane waves are 
present, they can combine into diffusons with non-zero transverse momentum, if the 
plane waves have different incoming angles. In that case the outgoing diffusons have 
a phase factor ex-p[—iQ ■ /?), which depends on the position on the outgoing surface p. 
The integration over the outgoing surface yields a delta function S[Q). In other words, 
only the momentum independent diffusons, where the phases precisely cancel, survive 
this averaging. The same holds for the frequency difference: The integrating sphere acts 
also as a time integrator due the random scattering inside it, yielding a factor 5[il). In 
conclusion, the only relevant outgoing diffusons have mass M = n. This observation 
is especially of importance in considering correlation functions. There the diffusons 
interfere, such that the mass of outgoing diffusons can differ from the incoming ones. In 
section 4.3.2 we will also calculate small 1/A^ deviations to this "selection rule". 



2.6 Scattering from one extra, static scatterer 

Following Berkovits and Feng [75], we consider the influence of one additional point 
scatterer on the diffuse intensity. This is a problem of practical importance as researchers 
try to locate an object in a strongly scattering medium by looking at the transmitted 
diffuse beam [76]. * We treat this problem within our diagrammatic approach. The 
problem can be rephrased as: 'how can two diffusons he tied together V In Fig. 2.6 we 
present all relevant contributions to the scattering. The two rightmost lower diagrams 
are self-energy diagrams that are not taken into account in the diffuson L. These 
diagrams were overlooked in Ref. [75]. However, their contributions are important as 
they cancel the leading term of the flrst r.h.s. (lower) diagram. 

The extra scatterer is described by a t-matrix te, with an albedo a^. It has position 
Te = (0,0, 2;e)- The diagrams again calculated using the tables in section 3.4. Let the 
diffuson propagators connected to the extra scatterer have momenta q\ and qi. Both 
momenta are pointing towards the extra scatterer. In lowest order of these momenta, 
we flnd for the lower diagrams of Fig. 2.6: 

A'l = Ue J d^pG(p + ^qi)G*(p - ^qi) / d^p'G(p' - \ci2)G*{p' + ^q^) 

-1 2 11 

~ tete[/i,i - -k^qll2,2 + ^k^lh^s + -k^h,! - -HO/1,2] X [qi ^ qa] 



_ i 



2 



t.t 



" '16^2 



®A different method for locating objects in diffusely scattering media is to measure the time-resolved 
transmission. The light that comes out first has scattered few times and mostly in the forward direction 
and is therefore not diffusely broadened, see references in Ref. [76, 77]. 
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Figure 2.6: Contributions from an additional scatterer to the difFuson intensity. C denotes 
the difFuson without the extra scatterer, the resulting diffuson with contributions from the 
extra scatterer taken into account. Thick lines denote particle propagators, thin lines denote 
hole propagators. The dashed line indicates that both the particle and the hole line interact 
with the extra scatterer. The — matrix of the extra scatterer is denoted by a cross. Its total 
effect is denoted by K. 

And for the second and third term one finds 

K, = K; = te|d3pG*(p)G(p + qi)G(p-q2) 

h,i + ^{ql + ql - qi • (\2)k^ h,i - ^HO(-/2,2 + 2/1,3; 



Stt 



-1 + + 92 - qi • q2)^' - ^^f" + 2(1 



The sum of the diagrams brings K = Ki + K2 + K^'- 



-[qrq2r + 3(1 -ae)] . 



f2.82) 



Corresponding to the first fine of Fig. 2.6, the expression for the diffuson with the extra 
scatterer included, reads in spatial coordinates: 



L{r, r') + -^[5{r' - Ve) + 5{r - re)]I(r, r') 
ntt 

tXf 

48^2 

xL(r2, r')5{ri - re)S{r2 - r^) . 



jdndr2[fVi-V2-'i{l-a,)]L{r,n] 



f2.83) 



We assume that the extra scattering is weak as compared to the total scattering Uetete ^ 
ntt, with Ue = 1/V, V is the slab volume. First we average the position of the extra 
scatterer over the whole slab. We expect to recover the diffusion equation (2.36) for a 
density of n + ne scatterers. Indeed, we find a reduced mean free path ie = Air / [ntt -\- 
netete)- Starting from equation (2.83) and substituting formula (2.36) we find to leading 
order in 



Le{q) 



L{q) + 2^L{q) + [q'f - 3(1 - a,)]^^4^L'{q) 



1? 



ntt 

[q^ + zOe + «e 



48^2 



(2.84) 



with Oe = f^^/^e, consistent with our definition of = —3Auj/£ve- Further, /€g£e 
3[ntt[l — a) + netete(l — ae)]/47r, involves a weighted sum of albedos. 
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The situation with the scatterer fixed is more interesting. We consider the transmis- 
sion of a plane wave through a non-absorbing medium. The transmission is (roughly) 
T = ^ '^^j^'^'* \z=L (we are only interested in functional form here). We calculate Le with 
the method of images [76]. Neglecting the second and third upper r.h.s. diagrams, the 
transmission in near field reads 

^ . ^ i2n + l)L- Ze 

Tip) = +2q£ y: 



L (p2 + [{2n + l)L - z^Yfl^ 



-OO 

^ ^ ^i^o (P^ + [(2n + l)i: - -eP)^/^ ' ^'-''^ 



in which p = (x, y), and 



Two cases are to be considered. If the extra scatterer does not absorb, only the p— term 
is present. This is analogous to a dipole in a static electric field between two capacitor 
plates. In the transmitted beam as function of p we see a wiggle if Ze > else we 

see a dip in the transmission. If the extra scatterer absorbs, the term is dominant 
and the scatterer acts as a drain for the intensity. The result is a dip in the transmitted 
beam. This would be equivalent to a (negative) charge in electrostatics. 

It is instructive to compare our approach with calculations within diffusion approx- 
imation. Den Outer et al. performed experiments where they located a pencil or glass 
fiber in a diffusive medium [76]. They found that the experiments are well described 
by a diffusion approximation. Diffusion was assumed both in the medium and inside 
the scatterer. The scatterer was for instance a sphere with radius i?, inverse absorption 
length /€2, and different diffusion constant D2- Consider the situation of weak uniform 
absorption and almost equal diffusion constants inside (-D2) and outside the sphere {D). 
Den Outer et al. find then: q = -kIR^{L - z^) l?,L and p = R^{D - D2)/3LD. We think 
of our physical scatterer as a sphere of radius R with a density Ue = S/AttR^ of extra 
scatterers. The extra scattering is weak as compared to the total scattering, if R is not 
too small, viz. R ^ [iteteY^^. The diffusion constant inside the sphere is D2 = ^ve^c- 
With this identification, our p and q values agree to leading order with the results of 
den Outer et al. 

In reference [75] the self-energy contributions presented in Fig. 2.6 are not taken 
into account. As a result the unphysical result is found that an extra scatterer (without 
absorption) act as a source of intensity. Obviously a Ward identity has not been satisfied 
in their approach.^ 

The vertex describing the infiuence of one extra scatterer can also be used to tie 
two diffusons together. This is sometimes useful, for instance when calculating the 



^Again one easily checks that cutting propagators in the dressed amplitude propagator yields the 
intensity vertices that we used and that respect the conservation laws. 
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Figure 2.7: Most crossed diagrams 



transmission the incoming diffuson can be attached to the outgoing one. For this purpose 
we take te = t. The vertex now obeys 

j dz2C{zi,Z2)KC{z2, zs) = -— j dz2[d^,^jCizi, Z2)]C{z2, zs) = jC{zi, Zs), (2.87) 
which is indeed the wanted behavior. 



2.7 Most crossed diagrams 

Another class of diagrams is also of some importance to us. They are the so called most 
crossed diagrams. We have drawn some elements in Fig. 2.7. These diagrams describe 
the enhanced return probability of diffuse intensity. They are unique to diffusion of 
waves (a particle cannot split itself in two, each part following a different path). The 
interference doubles the return probability. This process is the most important for the 
eventual Anderson localization, as was shown in Refs. [16, 19, 20]. In an electronic 
system, the presence of a magnetic field changes the phase of the loop. The return 
probability is lower, increasing the conductance. This effect is known as the negative 
magneto-resistance, see [31, 30]. 

The angular resolved measurements in optical systems allow for a direct measure- 
ment of the most crossed diagrams. If a beam impinges on a sample, the backscattered 
intensity in the incoming direction is almost twice as high as in other directions. The 
intensity as function of angle has a conical form with a width l/[ki). The extra inten- 
sity comes from the most crossed diagrams. The effect is also known as the enhanced 
backscatter cone. It provides a direct measurement of the most crossed diagrams. Due 
to the presence of single scattering contributions the peak value is somewhat lower than 
twice the diffuse background. 

Measurements of this effect were done by Kuga and Ishimaru [78], Van Albada and 
Lagendijk [79] and Wolf and Maret [80]. For theoretical work, see [81], and Refs. [76, 
59, 82, 83]. In this work, however, we do not treat the backscatter cone and only 
occasionally encounter most crossed diagrams. 



3 

Interference of diffusons: Hikami 
vertices 

After we have studied the diffuse transport in the previous chapter, we introduce in 
this chapter the other main ingredient for our calculations: the Hikami vertices. These 
vertices describe the interaction between diffusons. The Hikami vertices are named after 
Hikami, who used them in 1981 [24] for the calculation of weak localization corrections 
on the conductance. Interestingly enough their original introduction dates back two 
years earlier from the work of Gor'kov, Larkin, and Khmel'nitskii [84]. In the vertices 
diffusons exchange amplitudes, leading to the correlations of diffusons. 

3.1 Calculation of the Hikami vertices 

In this part of the calculation we again assume that the diffusons are slowly varying. 
The technique for calculating diagrams is well known [42]: Firstly, draw the diagrams; 
secondly, write down a momentum for each line; thirdly, use momentum conservation 
at each vertex; and fourthly, integrate over the free momenta left. 

The first step, writing down all leading diagrams requires some care. In Fig. 3.1 we 
have drawn the diagrams in the Hikami box in second order Born as already presented by 
Gor'kov et al. [84]. In second order Born one neglects scattering more than twice from 
the same scatterer. Therefore, the (dashed) interaction line connects two potentials. 
Yet we have not drawn all diagrams, one can imagine. Firstly, the box is attached to 
a diffuson ending on scatterer (which is consistent with the way we have defined the 
diffusons in chapter 2). As a result, common scatterings between propagators on the 
same leg are not allowed within the box. Secondly, diagrams with two or more scatterers 
and parallel dashed lines are sub-leading. They contain loops with two propagators of 
the same type, i.e. integrals like /d"^pG(p + q)G(p), which are sub-leading. Thirdly, 
diagrams with two crossed dashed lines are also sub-leading. Finally, we again work 
in the independent scatterer approximation. It is important to consider all corrections 
to the first term of the figure, as there is again a cancellation of leading terms. These 
cancellations are imposed by energy conservation. 

3.1.1 The four point vertex 

Let us calculate the vertex. As an example we calculate second r.h.s. diagram of Fig. 3.1 
in second order Born. As usual in the diffusion approximation where small derivatives 
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INTERFERENCE OF DIFFUSONS: HIKAMI VERTICES 



1 q1 




Figure 3.1: The Hikami four point vertex. It describes tlie excliange of amplitudes of two 
incoming diffusons 1 and 3 into two outgoing diffusons 2 and 4. The dots linked with the 
dashed line denote the dressing with an extra scatterer. The solid lines are dressed amplitude 
propagators. 



are assumed, we expand the Green's for small momentum and frequency difference. We 
use the same expansion as in the derivation of the diffusion equation. 



w 



to 



p +2p-q + g - A(jj— - nt 



id 



G + (-2p • q - + Alo--^)G' + 4(p • q)^^^ 



f3.11 



with G = G{y>]u). The momenta, numbered according to Fig. 3.1, point towards to 
vertex. We number the Green's functions also according to the figure. To lowest order 
in we have 



(6) 



f d"^p 

J ?27r)3^^'^^ + qi,tu + Ui)G3{p - q4,cu + cu3)G*(p,cu + U4] 



d^p' 



/dp 
— — ^Gi(p' - q2,CU + CUi)G3(p + q3,UJ + UJ3)G*2{p,UJ + UJ2] 



h,! + -iriql + 93 - q2 • q3)/4,i + -(^i + ^3 - 26^4)73,1 

3 c 



X [/2,l + . . . 



f3.2) 



where the /— integrals are given in the appendix. We eliminate nV^, in second order 
Born it holds that nV^ = Air/i. We take absorption terms into account in lowest order, 
i.e. only in the leading contributions, the /2,i terms. This yields 



+■ 



16^A;2 48^A;2 
A similar calculation gives 



q2-q3-qi-q4 + ^g,^ + K^) + 



-it 



[(jJl+(jJ3-(jJ2-(jJ4)- (3.3) 



(a) 



t t 
+ 



8^A;2 24^A;2 
3it 

+ ,0 (^1 + CU3 - 0^2 - CU4), 



qi • q3 - q2 • q4 - ^ q- + 



32^Fc' 

-t t 
+ 



16^A;2 48^A;2 
-it 

+ (^1 + CU3 - CU2 - CU4). 



qi • q2 - q3 • q4 + ^ + 



(3.4) 
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Figure 3.2: Diagrams contributing to tlie interaction of six difFusons: Hq. We did not draw 
possible rotations of the three rightmost diagrams, in total there are sixteen diagrams. 

Furthermore we define the reduced frequency of the legs as, Oi = —{iOi — O2 = 

— {iOi — U2)l O3 = — (cu3 — U2) I D and O4 = —{u^ — D. The sum of the three 
diagrams yields 



-2qrq3 - 2q2-q4 + Y^{ou + «i + 



^(Eq,). (3.5) 



This is the main result of this section. Note that the leading, constant terms proportional 
to £^ jk"^ have cancelled. In generalizing the earlier results, of for instance Hikami [24] 
or Stephen [85], we have included frequency and absorption terms. 

It is important to keep track of the terms. When the vertex is attached to a 
diffuson, the g| terms yield according to the diffusion equation, (g| + + iVlj)Cj = 1, a 
delta-function in real space. If we attach external diffusons to the Hikami box, this delta 
function is roughly one mean free path from the surface, see section 2.3.1 and can be 
neglected. This simplifies the expression for the Hikami box. As a result the absorption 
part and the frequency dependent part vanish. After the cancellations the expression 
is: 

-t 

^4(qi,q2,q3,q4) = TTTTTlqi' ^3 + f\2 q4) • (3.6) 



3.1.2 Six point vertex: H^, 



Below, we also need one higher order diagram: the six-point vertex ife, six diffusons are 
connected to this diagram. We depicted it in Fig. 3.2. This diagram was also already 
calculated by Hikami [24]. Also here, the dressings of the diagrams have to be added 
to the bare diagrams. Taking rotations of the depicted diagrams into account, there 
are 16 diagrams in second order Born approximation. It is not allowed to dress the 
bare six-point vertex (leftmost r.h.s. diagram in Fig. 3.2) with a scatterer that connects 
two opposite propagators. This dressing gives also a leading contribution, even if the 
dressing done with an arbitrary number of scatterers, but the resulting diagram is the 
same as the composed diagram with two four point vertices. Therefore, it should not 
be included in ife, but it will enter in diagrams with vertices. A calculation similar 
to the one above gives [24], 

Hq = — [qrq2 + q2-q3 + q3-q4 + q4-q5 + q5-q6 + qe-qi 



967rA;'^ 

^ 1 I 

+ E(q' + 2^^' + 

i=l 



<^(EqO- (3.7) 
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Apart from the generalization which includes frequency and absorption, Hikami's origi- 
nal expressions can be recovered from this using momentum conservation. 

Other forms of the Hikami-boxes are sometimes more convenient, especially the 
cancellation of short distance terms depends on the explicit form. That is to say, the final 
result is of course always the same, but the intermediate results can be much simpler. We 
will encounter such a situation in chapter 6. Using momentum conservation, c[i = 0, 
the different forms can be transformed into each other. A useful formula is 

=2^q,.q, + ^g,2 = 0. (3.8) 



3.1.3 Manual for Hikami boxes 



In our calculations we do not use the form of Eq.(3.5) for the Hikami boxes. Instead we 
use a Fourier transform in the z— direction, because in our slab geometry the [q^, Qy^ z) = 
(Q^z) representation is the most convenient. Now all become differentiations. The 
Hikami-four point box is 



96^A;2 



i=l 



(3.9) 



The differentiations work on the corresponding diffusons and afterwards Zi should be 
put z. Formally one has 



H4{z; Zi, Z2, Z-i, Z4)Ci{zi)C2{z2)C3{z3)C4{z4) 
= [dzrd,, + . . .]£i(zi)£2(^2)>C3(^3)A(^4)L,.= 

but most of the times we just write this as 

H4{Z)C,{Z)C2{Z)C3{Z)C4{Z). 



(3.10) 



(3.111 



The vertices thus give the spatial derivatives of the diffuson, that is to say the fluxes of 
the diffusons. The observation is helpful to estimate the influence of internal reflection, 
which lowers the spatial derivatives of the diffusons. Note that all leading terms in H4 
and He dropped out; this is again closely related to the conservation laws [86]. 



3.2 Beyond second order Born 

Going beyond second order Born approximation in the calculation of the diffuson dia- 
grams resulted in a simple replacement of the mean free path. We recall that it became 
i = ntt/47r, instead of £ = n|yp/47r. The diagrams of the diffuson were slightly different 
but this was only a superflcial problem, as the same calculation remained applicable. 
However, for the vertices the situation is more subtle. Previous section we only worked 
in second order of the scattering potential. Therefore, we calculated the box up to the 
same order, leading to the flrst three contributions of Fig. 3.3. As we work in all orders 
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Figure 3.3: The Hikami four point vertex, tlie Hikami box, beyond tlie second order Born 
approximation. Instead of tliree, tliere are now eiglit diagrams to be calculated. The resulting, 
however, is apart from a renormalization of the mean free path the same as in second order 
Born approximation. 

of the potential all contributions have to be calculated. To leading order in 1/A;£, these 
are the eight diagrams depicted in Fig. 3.3. 

The calculation goes very similar to the second order Born above. The only difference 
is that instead of the t— matrices t and t occur [t on G's, t on G*'s). Already for the 
cancellation of the leading term, the three diagrams are not sufficient, but one needs all 
eight diagrams. One finds 



(h) 



2qi • qs - 2q2 • q4 + J^i^l + ^1 + ^' 



f3.12) 



This is the same expression as found in the second order Born approximation ! Note, 
however, that now the definition of the mean free path is different: i = ntt/4:TT instead 
of i = nV"^ jAiT. One can conclude that we only renormalize the mean free path if we 
go from second order Born to the full Born series. This can also be shown on a much 
more general level in the non-linear sigma model [87]. The number of diagrams increases 
substantially as compared to second order Born. Most of the times this is not a problem 
as once we did the calculation, we only need the sum of diagrams, Eq. (3.12). (The 
calculation in chapter 5 is a different story as we will undress the boxes again. There we 
will stick to second order Born.) For the six-point vertex we have not checked whether 
the second and full Born approximation give the same result; there are (at least) 64 
diagrams (excluding rotations)! 
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Figure 3.4: The first order correction to tlie conductivity in tlie Hikami formalism. On tlie 
riglit we liave written out tlie vertices in detail in second order Born approximation. On the 
second line the first r.h.s. diagram is redrawn, one sees that it is a loop effect. 

3.3 Corrections to the conductivity 

In the original works of Gor'kov et al. and Hikami the vertices are introduced to calculate 
to weak-localization corrections to the conductivity. We already said that the most 
crossed diagrams give the largest correction to the conductivity. They increase the 
return probability of the intensity above the value obtained from diffusion theory. In the 
Hikami vertex formalism we have drawn these diagrams in Fig. 3.4. Hikami calculated 
these corrections up to second order, i.e. with two Hikami-boxes. In this work we do 
not consider this type of higher order loop effects. We suppose that we are far from 
localization, so that we can restrict ourselves to the leading processes. Then the vertices 
only show up in the interaction between separate diffusons. 



3.4 Appendix 

We calculate the integrals 

d^p 



h,i ^ j -^G'=(p)G*'(p), (3.13) 



where 

G(p) = [p2 - A;2 - G*(p) = [p2 -A;2 (3.14) 

For the simplest integral, Ii^i we find 

The sum of the residues yields l/47r(/i + /!), and the optical theorem brings 

—init — t) n Imt ntt 1 , , 

l^+¥= h - = ^^ = ^ = -„- 3.16 



Therefore, 
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Or, with absorption, /i + /i = [£(1 — /€^£^/3)] ^, yielding 

h 1 = — : r = — 1 - -K^l^ 

' 47r(/i + 71) 47r V 3 

The Ik^i can be found from the Ii^i integral as 

-Id, 



-1 d 



With differentiation we now find without absorption 

h,3 - 

h,4 - 
h,3 - 



1 t 




1 ^, f 


47r(/Lj+7i) 47r ' 


-^2,2 — 


87r/Lj/Lj(/Lj+/Lj)^ SttA:^ 


1 -if 




1 ^: if 




^2,1 = 




1 






167r772(/Lj+77)3 ~ 167rfc2 5 


^3,1 = 


167r/Lj2(;Lj+77)3 ~ 167rfc2 


3 -Si^'' 




3 3i^'' 


327r/Ljl772(/Lj+77)4 327rfc3 ' 


-^3,2 = 


327r/Lj277l (/Lj+77)4 327rfc3 


1 i^'' 




1 ^, -J^" 


32^77" (/Lj+77)4 327rfc3 5 


h,l = 


327r/Lj3(/Lj+77)4 ~ 327rfc3 


1 




1 ^, -f 


167r/Lj77''(/Lj+77)5 ~ 167rfc4 5 


h,2 = 


167r/Lj377(/Lj+77)5 ~ 167rfc4 


3 3^5 













We will also use 



d^p 



p-qfG'=(p)G*'(p) = ^fcV/- 



that has a factor 1/3 from the angular average. 
With absorption there is generally a prefactor 



oc (1 - K^ei?,) 



\n-\-m — l 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



where we have assumed that n£ <^1. However in our approximation this prefactor only 
has to be included for terms not proportional to g or 0, as we work in first order of ML 



3.4.1 Beyond the diffusion limit 

In the previous part we assumed that g£ ^ 1, in real space this assumption means that 
variations of the diffusons are small on the length scale of one mean free path. (This 
is the same approximation that we made for the diffusons ). Yet this inequality is not 
always fulfilled, sometimes variations are large on mean free path. This is can happen if 
we integrate over the distance between two Hikami boxes, this we will see in chapter 5. 
The diffusion approximation results there in a divergent term. (In principle the diffusion 
approximation does also not hold near the boundary, but we won't consider these small 
effects.) We will need the more general integrals /^", which we define as 

4T(?) ^ / ^,GHv)G*\v)G-{v + q)G*"(p + q). (3.21) 
With the properties 
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Table 3.1: The short distance diagrams can be factorized in these integrals, defined in 
Eq.(3.21); A, is defined in Eq. (3.23). 



In the calculation of the diffuson and the Hikami vertices in the previous sections we 
expanded the integrals in but if we are after contributions for g ~ l/£, this expansion 
is not allowed. We can still assume that g ^ A;, as we never will need the physics on 
length scales comparable to the wavelength, but comparable to one mean free path only. 

The integrals needed in the calculation of chapter 5 are given in Table 3.1. We define 
the angular average Ai as 

The integral for i = 1 yields the diffuson kernel Ai = arctan(g£)/g£. The internal 
diffuson now reads 

We recover previous results by expansion in (qi). This yields C-mt = 127r/(£"^g^), in 
agreement with the bulk solution of (2.39). 

The integrals relate to the integrals of the previous section when we expand again 
in qi 



^kl — ^{k+m),{l+n) + g9 [rrml(^k+m + i),{i+n+i) 

+ ^m{m + l)I^k+m + 2),{l+n) + ^n{n + l)I^k+m),{l+n+2) 



(3.25) 



In the non-linear sigma model the vertices in the expansion beyond q^ terms are often 
called higher gradient vertices, as they correspond in real space to higher order deriva- 
tives. It was noted by Altshuler et al. [88, 23] that they become important in the higher 
moments of the conductivity, see also chapter 7. 
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Short and long range correlations: 



4.1 Physics of the correlation functions 

A nice feature of optical mesoscopic systems above electronic systems is the possibility 
to measure several transmission quantities. This depends on whether integration over 
incoming and/or outgoing directions is performed. In electronic systems one usually 
measures the conductance of the sample. This is done by connecting incoming and 
outgoing sides to a clean electron "bath". All electrons are collected and all angular 
dependence is lost.^ In optical systems one usually measures angular resolved, but an- 
gular integration is also possible. The correlation functions of the different transmission 
quantities, their magnitudes, their decay rates, and even the underlying diagrams are all 
very different. In this chapter we discuss the angular transmission and total transmission 
correlation functions. 



4.1.1 Angular resolved transmission: speckle 

If a monochromatic plane wave shines on a disordered sample we see large intensity 
fluctuations, so called "speckles", in the transmitted beam. The speckle pattern is 
wildly fluctuating (as function of the frequency of the light or the outgoing angle). 
Simply shining a laser on the (rough) wall already produces such a speckle pattern. The 
typical diameter of a spot at the outgoing surface is A. In this chapter we are interested 
in the correlation between speckle patterns of two different beams. The beams may have 
different incoming angles, different frequencies or different positions. 

There are short, long and "inflnite" range contributions to the correlation function of 
the speckles. In a recent review by Berkovits and Feng [90] the different correlations and 
their physical interpretation are discussed. We denote the transmission from incoming 
channel a (wave coming in under angles (pa) to outgoing channel b (waves transmitted 



^An exception is the recent measurement of Gao et al. [89]. They injected electrons into a two- 
dimensional sample using a so called quantum point contact. A quantum point contact is a very small 
connection between two metals, such that only one electron mode can be selected. The electrons were 
also collected from the sample with a quantum point contact. A fluctuation very similar to the speckle 
pattern was seen. 
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SHORT AND LONG RANGE CORRELATIONS: Ci, C2 



into angles 6^^ (j)b) as Tab- the correlation functions can then be classified as [90, 91] 



Feng et al. [91] originally put this expression forward as an expansion in l/(/, the Ci 
being of order one, the C2 of order l/(/, and the C3 of order Ij ■ We present a sketchy 
picture of the diagrams of the different correlations in Fig. 4.1. The largest part of 
transport are the independent diffusons, yet correlations are present, mixing diffusons 
with different frequencies or angles. 

The unit contribution in Eq. (4.1) just comes from the uncorrelated product of 
average intensities. The Ci-term in the correlation function is the most important one 
if the angular resolved transmission Tab is measured. The Ci is namely unity when both 
the incoming directions a and c, and the outgoing directions h and d are pairwise the 
same, and Acu is zero. If one changes the frequency of the incoming beam, the speckles 
of the outgoing beam deform. Eventually the speckle disappears and the correlation 
vanishes exponentially. This effect is the short ranged Ci contribution. Also of interest 
is the case of the angular Ci correlation. This time we keep the frequency fixed. If we 
change the angle of the incoming beam, the speckles change due to two effects. First, 
the speckle follows the incoming beam. This is also known as the memory effect, which 
one can even see by the naked eye. Secondly, the speckle pattern also deforms, i.e. it 
de-correlates, this is again the Ci. The Ci is a sharply peaked function non-zero only 
if the angle of the incoming and outgoing channel are changed by the same amount. 
Diagrammatically, the Ci factorizes in two disconnected diagrams, see Fig. 4.1. Just as 
the unit contribution in Eq. (4.1), it is equal to the product of two averages. 

Theoretical studies of the Ci were first done by Shapiro [92]. Experimental work on 
the Ci as function of angle, explicitly showing the memory effect, was done by Freund 
et al. [93]. Albada et al. studied the frequency dependence of the Ci [71]. Later, effects 
of absorption and, especially, internal reflection were studied and turned out to be of 
practical importance [94, 68]. 

4.1.2 Total transmission correlation: C2 

The C2 gives the long range correlations in the speckle pattern. Consider again a 
single direction in, single direction out experiment. Yet this time we look at the cross 
correlation of two speckles far apart. As there is an angle difference in the outgoing 
channel, though not in the incoming beam, Ci is now absent. Instead one sees a much 
weaker correlation. As the frequency shifts, this correlation decays algebraically. This 
is the C2 contribution, it describes correlations between speckles far apart. In a single 
channel in, single channel out experiment it is only possible to see the weak effects of 
these higher order correlations in very strongly scattering media [95], i.e. if the sample 
is rather close to Anderson localization. The C2 can be measured more easily in a set-up 
using one incoming channel and integrating over all outgoing channels. By collecting 
the outgoing light, one measures the total transmission = 2^afc- In this set-up 
the C2 correlation function, which has a much smaller peak value but is long ranged, 
contributes for all outgoing angles. The C2 is now the leading correlation as the sharply 



(Tab{'^)Tcd{'^ + Acu)) 
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Figure 4.1: A schematic picture of the different correlation functions present in the transmis- 
sion of two intensities. The arrows distinguish between advanced and retarded propagators. 
Propagators with equal transverse momentum or frequency have the same style of line (dashed 
or solid). The black box contains in principle all contributions to (TabTcd)- The main part just 
factorizes, but correlations are present. Below, we will present the diagrams for the C2 and 
C3 in more detail. 

peaked Ci is overwhelmed. As we already said in section 2.5, only outgoing diffusons 
with no transverse momentum and no frequency difference are leading in the total 
transmission. This way the phases of the outgoing amplitudes are exactly opposite and 
thus cancel. From Fig. 4.1 one sees indeed that it holds for the C2 diagram (outgoing 
lines of similar style pair), but not for the Ci diagram. The long range character arises 
due to interference of the diffuse light paths. The C2 correlation, which still depends 
on the angles of the incoming beams a and c, is of order . The C2 corresponds to 
a diagram where the two incoming diffusons interact through a Hikami-vertex. In this 
vertex the diffusons exchange amplitudes. Rather than thinking of this exchange as a 
real physical process, we consider this more like something left after subtraction of the 
average process, see Fig. 4.1. 

The long range C2 correlation function was first studied by Stephen and Cwilich 
[96]. Zyuzin and Spivak introduced a Langevin approach to simplify the calculation 
of correlation functions [97]. Pnini and Shapiro applied this method to calculate the 
correlation functions of light transmitted through and refiected from disordered samples 
[98]. The C2 correlation functions were measured in several experiments. For optical 
systems Albada et al. performed measurements [71, 72]). Microwave experiments were 
performed by Genack et al. [99], [100]. In several papers effects of absorption and 
internal reflection was studied, see Pnini and Shapiro [101], Lisyansky and Livdan [53], 
Zhu et al. [54], and Ref. [68]. They showed that absorption and internal reflection, 
neglected in the earliest calculations, signiflcantly reduce the correlations. 

Finally, the C3 term in Eq. (4.1) is dominant when the incoming beam is diffuse. 
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Figure 4.2: A contribution of order to tlie C2 correlation. By following the amplitudes, 
one can check that the incoming pairings ij* and ji*, change on the outgoing side into pairings 
ii* and jj*, as is required for C2 diagrams. It is thus not contributing to the C3, which is also 
of order but requires diagrams without amplitude exchange. 

and one collects all outgoing light. Then one measures, just as in electronic systems, the 
conductance g = (T) = ^ Tab- In that measurement contributions to the correlation 
where a and c are arbitrary far apart are dominant. Though C3 is of order (/~^, it 
dominates over the Ci and C2 terms as it has contributions for all incoming and outgoing 
angles. Therefore it is sometimes called the 'infinite' range correlation. In contrast to 
the C2, now also the incoming amplitudes must have opposite phase. This occurs in 
a diagram where the two incoming diffusons interact twice. Note that a loop occurs, 
see Fig. 4.1. Apart from the normalization to the average in the C3, the C3 describes 
just the well known Universal Conductance Fluctuations, or UCF [102, 103, 104]. We 
discuss the C3 in detail in chapter 5. 

In principle the C2 and C3 correlations are present as a sub-leading term in the 
angular correlation function, but there are also other contributions of order 1/g and l/(/^, 
respectively, to the correlation Eq.(4.1). The weak localization correction, see section 
2.7, is an example of a l/(/ contribution to the Ci. We have drawn a contribution 
to the C2 in Fig. 4.2. In a paper by Garcia et al. [100] these corrections are studied. 
Therefore, the interpretation of the C's as an expansion in Ij g is somewhat impractical. 
Rather we define the Ci, C2, and C3 as the leading term in the correlate of Tab, 
and T, respectively. Equivalently, the Ci contains disconnected diagrams; incoming and 
outgoing amplitudes must have the same pairing, the C2 contains connected diagrams 
swapping the initial pairing, and the C3 contains the connected diagrams in which 
incoming and outgoing pairings are the same. 

4.2 Calculation of the Ci correlation 

The calculation of the Ci is rather simple as it is just the product of two independent 
diffusons with exchanged partners. Both diffusons consist of one amplitude from one 
beam and one complex conjugated amplitude from the other beam. Due to momentum 
conservation in each diffuson, the perpendicular momentum difference of the two incom- 
ing beams, is equal to the difference of the outgoing momenta. Therefore, the Ci is 
only non-zero if incoming and outgoing angles are shifted equally during the experiment. 
We calculate the normalized product of a diffuson with 'squared mass' = Q"^ -\- 
-\- iil and the complex conjugated diffuson. Due to outgoing Green functions, Ci is a 
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Figure 4.3: Ci angular correlation function plotted against the scaled perpendicular momen- 
tum difference. Solid line: small zq, no absorption. Short dashed line: absorption (k = 2/L), 
small Zq. Long dashed line: large skin-layers {zq = L/3), no absorption. 



Laplace transform of the diffusons. For small mass [Mi ^ 1 ) the same integrals occur 
in the numerator and denominator. The result can be expressed by the values on the 
boundary z = L. This yields: 



Ci(M) 



Tab{K 

MP 



jC{z = L,M) 



C(z = L, k) 



1 -|- K z^ ] sinh kL -\- 2kzq cosh kL 



[I + M'^zl) sinh ML + 2Mzq cosh ML 



(4.2) 



In Fig. 4.3 we plotted this function for various situations. In thick samples in the limit of 
Mzq <^ 1 and no absorption this reduces to the result already known from the diffusion 
approximation [91] 

ML 



Ci(M) 



sinh ML 



(4.3) 



For large angle or frequency difference (|M|L ^ 1), Ci decays exponentially, as Ci{M) ~ 
g-2Re(M)L_ The spatial correlations of beams can be extracted from the Fourier transform 
of the angular correlation function. 



4.3 Calculation C2 correlation 

The diagram of the long range correlation C2 is depicted in Fig. 4.4, Remember that the 
Ci is zero if incoming and outgoing angle are changed unequally, but the C2 is connected, 
allowing momentum flow from one diffuson to the other. It is the simplest diagram that 
is non-zero if one changes incoming and outgoing angle by a different amount. C2 stems 
from the interaction between two diffusons, which exchange partners somewhere inside 



46 



SHORT AND LONG RANGE CORRELATIONS: Ci, C2 



Figure 4.4: The diagram of the long range C2 correlation function. The shaded box denotes 
the Hikami four-vertex, worked out in Figs. 3.1 and 3.3. 

the slab. The Hikami box describes this exchange. The full expression reads 
1 

C2 = irT.\lrT^\ I dzdzidz2dz3dz4Ci{zi,Mi)C3{z3,M3)H4{zi,Z2,Z3,Z4) X 

{T}{T} Jo 

S{Z - Zi)S{z - Z2)S{Z - Z3)S{Z - Z4)C2{L - ^2, M2)C4L - Z4, M4) , (4.4) 

where we used the fact that also C2 can be expressed as the ratio of the values at the slab 
boundary. We have labelled the incoming diffusons with odd numbers, and the outgoing 
diffusons with the even numbers; a convention used throughout this work. Attaching 
diffusons to the Hikami box allows a simplification of the expression for the box. We 
calculate the C2 by inserting = Q'^ -\- k'^ i: zO, Mi = M3 = n. In the 2;— coordinate 
representation the box then reads, see chapter 3, 

+ d,^d,^ + Q ) ■ (4.5) 

487rA;'^ 

This gives one contribution to the C2, another is obtained by interchanging incoming 
and outgoing beams. But as in a typical C2 experiment one integrates the outgoing 
beam over all directions, this second term averages out. 

The above expression for C2 is similar to the one obtained in the Langevin approach, 
as one can see by comparing formula (39) of reference [98] to (4.5) and (4.4) in this chap- 
ter. Indeed, in the Langevin approach one assumes that there is a macroscopic intensity, 
which describes the average diffusion, and an uncorrelated random noise current super- 
posed. In our approach we have a slowly varying diffuse intensity, while the point-like, 
uncorrelated random interactions originate from the Hikami-box. 

We obtain the resulting general expression for C2 by inserting the appropriate diffu- 
sons and the Hikami-box Eq. (4.5) into Eq. (4.4). The result is rather lengthy and given 
in the appendix. 

Let us first study the top of the correlation function. By its definition the top of 
the correlation functions corresponds to the second cumulant of the total transmission 
distribution function. In the plane wave limit of the incoming beam, all transverse 
momenta are absent, we also neglect absorption for the moment. Mathematically these 
approximations correspond to all M being equal zero. The diffusons are now simple 
linear functions, given by Eq. (2.61) and (2.76). Neglecting internal reflection, one 
obtains the known result for the second cumulant [96] 
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= (r„)-2| I dxdy £ dzH,C,{z)C2{z)C3{z)C4z) 

where double brackets denote cumulants normalized to the average. The correlation is 
proportional to the inverse power of g counts the number of four-point vertices. 

It can be looked upon as the chance that two intensities interfere. Later we see more 
calculations supporting this interpretation. 

We discuss some other simplified cases. Consider the case of angular correlations 
and no absorption but taking the boundary effects into account: 

C2iQ) = ^HQL) , (4.7) 



in which we define the dimensionless function F' 



2 



F2{QL) = [smh2QL-2QL + Qzo{6smh^QL-2Q''L'') 

+AQ'^zl {s\n\,2QL - QL) + QQ^zl sinh^ QL + Q^z^ s\n\,2QL\ x 

[2QL {(1 + Q'^zl) sinh QL + 2Qzo cosh Qlf]'^ . (4.8) 

which decays like IjQ for large Q [i.e. large angles). However, note that it still should 
hold that Q£ <^ 1. Neglecting boundary effects we find the result known from diffusion 
approximation [101]: 

sinh(2gi:) - 2QL 
2QLsmh^{QL) 

In figure 4.5 we have plotted the F2 function for various situations. 

Above, we have studied the correlation as function the two-dimensional momentum 
Q. One also can study the real space correlations, which relate the momentum space 
correlation by a Fourier transform. The functional form in real space is of roughly 
similar shape as Fig. 4.5, with typical decay length L [96]. 

We can also obtain frequency correlations. Without absorption, neglecting boundary 
effects, one finds [72] 



2 / sinh V20 L - sin V20 L \ 
F2 = ^ ' ' . 4.10 

As one can see in the appendix, the functional form changes when including boundary 
effects. 

Without internal reflections (m = 1) at large M, C2 correlation function decays as 
1 /M. Lisyansky and Livdan found that this behavior changes at high internal reflectivity 
[53]. We do not conflrm such a behavior, as at large Zq the decay remains of order 1/M. 



4.3.1 Focus 



Let us now study the influence of the beam proflle on the correlation. If the spot of 
the incoming beam is flnite, amplitudes with different angles, i.e. different transverse 
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Figure 4.5: F2 angular correlation function plotted against the scaled perpendicular momen- 
tum difference. Solid line: small zq, no absorption. Short dashed line: absorption (k = 2/L), 
small Zq. Long dashed line: large skin-layers {zq = L/3), no absorption. 

momenta are present. They can combine into incoming diffusons with perpendicular 
momentum. We suppose that the incoming beam has a Gaussian profile. We decompose 
it into plane waves defined in Eq. (2.48) (for convenience we assume the average incidence 
perpendicular) 

= ^E<^(^«)c, m) = ^e-«^^^/^ (4.11) 

where po is the beam diameter. In order to have two diffusons with a momentum Q and 
— (5, we find that the four incoming amplitudes combine to a weight function 

f d'P.d'Ps ct^iPiWiPi + Q) (piPsWiPs -Q)= e-^oQV4. (4_^2) 

We get the correlation function by integrating the (^-dependent correlation function over 
the momentum with the corresponding weight. 

We again neglect boundary layers. For a Gaussian beam profile we find for the top 
of the correlation [72] 

{{Tl)) = ^J<i'Q e-'^oQy^F^iQL), (4.13) 

with ^2(3;) = [sinh(2a;) — 2a;]/[2a; sinh^ x]. If the incoming beam is again very broad, po ^ 
L, only the term F2{Q = 0) = 2/3 contributes and one recovers the plane wave behavior 
{{T^)) = 2/3g. Note that this agreement is found by identifying the area of a Gaussian 
profile with A = ttPq. The second cumulant decreases as 1/pq at large po. In a real 
space picture evidently the correlation increases if the two incoming channels are closer 
to each other, i.e. if the beam diameter is smaller. At smaller beam diameters the top of 
the correlation is proportional to l//9o, and would diverge in our current approximation, 
in which Qi ~ i/ po ^ 1. In the experiments of De Boer, Van Albada and Lagendijk [71, 
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a) 



b) 



Figure 4.6: The two contributions to tlie second moment of tlie total transmission. In diagram 
(a) the transmission channels are independent; this process is of order unity and is almost 
completely reducible to the mean value squared. Diagram (b) corresponds to two interfering 
channels; this is the second cumulant or C2 diagram, it is of order 1/g. The close parallel lines 
are diffusons. 



72, 105], they kept the focus small in order to minimize the dimensionless conductance 
g and therefore to maximize the fluctuations. 



4.3.2 Contributions from disconnected diagrams to C2 

So far we calculated the leading contributions to the long range correlation. They are 
the connected diagram in Fig. 4.6(b). Yet there are also contributions to the C2 from 
the disconnected Ci diagrams of Fig. 4.6(a). The disconnected diagram corresponds 
to cumulant contributions that are not due to interference inside the slab. It describes 
effects that have little to do with the interference effects we are after. Here we calculate 
the contribution and show that it is small. 

We discuss the contribution of these diagrams only for the top of the C2 correlation 
function; due to the much slower decay of C2 as compared to the Ci the correction 
should be the largest at the top of the €2- 

As a start we use the model of a waveguide. We assume that the disorder couples 
one incoming mode a to all outgoing modes. A waveguide has discrete modes, and for 
the moment we assume that different outgoing modes are uncorrelated. The second 
moment (T^) splits into a connected part {T^)com Fig. 4.6(b) and a disconnected part 
(r^^)dis, Fig. 4.6(a). The total transmission is the summation over all outgoing modes 
Ta = J2bTab- The disconnected part of the second moment is 

N,N N 
{Ta)dis = ^ {TabiTab-,) + ^ {Tab^Tab-,) 
bi ^62 61 =62 

= iv(iv-i)(r„,)2 + iv(r„\) 

= (r„)2 + iv(^„fc)^ (4.14) 

where N is the number of modes supported by the waveguide. For the last equality 
sign we used that the averaged second moment of the intensity speckle is given by the 
speckle distribution function Eq. (7.1), (T^^) = 2{TabY ■ From Eq.(4.14) we see that 
the disconnected diagram {T^)dis does not completely factorize into the average squared 
{TaY . Therefore it contributes to the second cumulant. As shown above, the connected 
part of the second moment of the total transmission {T^)con is proportional to L/N£. 
For the sum of the disconnected and the connected contribution to the second cumulant 
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Figure 4.7: The disconnected contribution to tlie second moment. In (a) tlie two transmissions 
factorize into the average value squared, and does not contribute to the second cumulant. 
Diagram (b) is much smaller, but gives a contribution to the second cumulant. The amplitudes 
making up the two difFusons propagate out in different directions. 



one thus finds 

{{T^)) = {{T^))con + ((r„^))d. = ^ + (4-15) 

which also holds for plane wave case. 

Now we turn to the more realistic situation of a diffusely scattering slab, with a finite 
focus of the incoming beam. The incoming beam broadens in the transverse direction 
by diffusion, changing the above result. To calculate the different contributions to 
the cumulants, we need the intensity distribution at the exit interface at transverse 
coordinates Ri and R2. 

The amplitudes, making up each diffuson, can propagate from the outgoing surface 
in different directions, respectively Pi, P2, P3 and P4. 

^{Ri)e'^'^'^^Ri)e-'^'^'^{R2)e'^'^'^*{R2)e-'^'^\ (4.16) 

To get the contribution to the second moment of the total transmission, one first inte- 
grates over the transversal coordinates Pi and P2 to get the contribution of the whole 
exit interface to the intensity in a certain direction, then one integrates over all direc- 
tions to get the total transmission. To obtain intensities the amplitudes need to be 
paired giving the following possibilities. The first possible pairing of the amplitudes is 
Pi = P2, P3 = P4, see Fig. 4.7(a), and brings 

J dPi J dP2 (*(Pi)**(Pi)*(P2)**(P2)) = /'(O), (4.17) 

where we have defined I{Q) as the transmission by a diffuson with transverse mo- 
mentum Q. Including the incoming Gaussian beam profile, it is proportional to, see 
Eq.(2.42),(4.11) 

|n|e-Q'po/8 



Integrating over all outgoing directions results in 



(4.19) 
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But this is just the factorizing contribution {Ta){Ta), and therefore does not contribute 
to the second cumulant. 

The second possible pairing of the amplitudes, Pi = P4, P2 = P3, does give a 
contribution to the second cumulant. It is the diagram in Fig. 4.7(b) 

j AR,AR2 (*(Pi)**(Pi)*(P2)**(i?2))e'''^(^^-^^)e-'«^(^i-^^)/(Pi - P2)/(P2 - Pi). 

(4.20) 

Subsequent integration over the outgoing directions Pi and P2 yields 



Tl) = f d'Q I\Q). (4.21) 

J\Q\<k 



l\Q\<k 

The integral can be extended to infinity since I{Q) is an exponentially decaying function. 
The contribution of the disconnected diagram Fig. 4.6(b) to the second cumulant is thus 

u^a;;dis- ^2^4/2(0) - N' ^ > 

In analogy to the waveguide, this result describes irreducible contributions from dis- 
connected diagrams. It can be interpreted as the inverse of the number of independent 
speckle spots in transmission at the exit interface [106, 105]. 

We did not take the vector character of the light into account yet. The two in- 
dependent polarizations of each outgoing direction effectively double the number of 
independent speckle spots with respect to Eq. (4.22), thus reducing Eq. (4.22) by a 
factor 2. 



Reflection correlations 

Similar to transmission correlations, also reflection correlations can be studied. We have 
seen before that the transmission corresponds to the diffusons, but in the reflection there 
are two intensity diagrams: The diffusons bring a constant background, the most crossed 
diagrams yield the enhanced backscatter cone, which is only of importance if incoming 
and outgoing angle are the same. This leads to more diagrams and to more peaks in 
the reflection correlation than in the transmission correlation. 

For the Ci one can both pair diffusons and most crossed diagrams, see Berkovits et 
al. for some theory [107]. Indeed in optical experiments two peaks were seen by Freund 
and Rosenbluh [108]. 

The C2 in reflection was treated by Berkovits [109]. Measurement of the C2 in 
reflection by collecting all reflected light is tricky as one does not want to interfere with 
the incoming beam. Nevertheless, the measurement is may be possible as a long range 
component in the angular resolved reflection (this however requires quite small values 
of g). 

For both Ci and C2 in reflection we don't expect a very good agreement with the 
theory. The problem is that the precise behavior of the diffuson near the surface is 
important. The diffusion approximation g£ ^ 1 is not valid anymore, in particular 
if the indices of refraction match (then the intensity gradient near the surface is the 
largest). The calculation of the Ci is probably easily extended, yet for the C2 one also 
would need the Hikami box beyond flrst order in qL 
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4.4 Comparison with literature 

When comparing our findings to results in literature, care has to be taken of the def- 
inition of the beam profile as this changes the prefactor of €2- Also the definition of 
the transmission T is different is most places, and not all results are normalized by the 
transmission as done here in formula (4.1). A diagrammatic analysis has previously 
been carried out by Stephen and Cwilich [96]. Their result exceeds our equations (4.7) 
and (4.9), by a factor 9. Pnini and Shapiro used the Langevin approach to calculate 
C2 in the absence of internal reflections [98]. Their result agrees with ours in the limit 
Zq <^ L. It can also be seen that equation (4.10) coincides with the Langevin approach 
of de Boer et al. [72]. In conclusion, we have shown that the diagrammatic and the 
Langevin approach lead to the same result for the calculation of €2- We calculated 
the short and long range, angular and frequency correlations. Previous approaches are 
extended by the presence of absorption, skin layer effects and scattering beyond second 
order Born approximation. 



Appendix 

Here we present the general expression for the long range correlation function, given by 
formula (4.4) with the diffusons of formula (2.64). The resulting expression is a fourth 
order polynomial in Zq. It is convenient to split M2 and M4 in their real and imaginary 
part, M2,4 = ^(5^ + lb = a ± ib. 



(g2 + K^)2 + 02 cos((^/2), b = ■{/(g2 + ^2)2 ^ Q2 sm{ip/2) 

_n 

g2 + «2 



Ump = 77:— —7. (4.23) 



We flnd 



- ¥lA BL ' ^^-^^^ 



where we have deflned: 



B = [1 + 2z^{a^ -b^) + z^{a^ + b^)]{cosh2aL- cos 2bL) 

+2zo[l + z^ia^ + 6^)](asinh2aL + bsm2bL) + Az^ia^ + 6^)(cosh2aL + cos 2bL) 

Na = —-^ X {[(2a^ - K^) + z^(8a* + 2a%^ -b\^ -3a\^ + K*) + 

2a[a^ — K^) 

z^{a^ + b^){2a^ - 2a\^ + k^)] sm}i2aL + 

[4:Zoa{3a^ - k^) + Az^aiSa'^ + a%'^ - 2a^K^ + k'^)] sinh^ aL} 

Nk = -—J-—-x{[-{2b' + ^) + z',{2a%' + 8b^ + a\' + 3b\' + K^)- 
2b[b'^ + K,'^) 

zl{a? + 62)(264 + 2b^^ + n% sin 261 + 

[4zo6(3&' + «') - AzlbiZy + + 2b^n^ + n% sin^ bL\ 

N. = 2,(,2 :l2)(^2 ^ ,2) X {[-^ - 4i^a%^ + 5a-.- - 56^.^ + + 
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[4zo{-a%^ - + - k*) + 4z^{-3a%^ - 2a\^ + 2b\^)] srf kL} (4.25) 

Without internal reflection we flnd 

_ 37r 1 r(«^ - 2a2) . 

' " 2mA {cos{2bL) - cosh{2aL)) ^ a{a^ - k^) ^^^^^^aL) 

(k^ + 2P) , , , Kia^ + 62) , , , , 
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Conductance Fluctuations: C3 

In this chapter we study the optical analog of the universal conductance fluctuations 
in metals. The optical conductance of a multiple scattering medium is the total trans- 
mitted light of a diffuse incoming beam. This quantity, analogous to the electronic 
conductance, exhibits conductance fluctuations. In the electronic case these fluctua- 
tions are known as universal conductance fluctuations (UCF). The optical conductance 
and its fluctuations are difficult to measure and measurements were not performed untill 
now.^ Yet measurements of this quantity would form a cornerstone in the analogy be- 
tween optical and electronic mesoscopic systems. We again use the Landauer approach, 
which we used also for Ci and C2, see chapter 4. However, when calculating the C3 in 
the Landauer approach, unphysical divergencies show up and showing their cancellation 
turns out to be elaborate. 

First, we review the electronic UCF. In section 5.2 we present the long distance 
diagrams and analyze the divergencies arising from these diagrams in the diffusion ap- 
proximation in 5.3. In section 5.4 we show that the divergencies cancel if we analyze all 
diagrams in detail. We calculate the general form of the fluctuations in section 5.5 and 
applied to optical systems, where absorption and internal reflections may be present. 

5.1 Introduction 

Let us briefly review the situation for mesoscopic electron systems. The electronic 
conductance of mesoscopic samples shows reproducible sample-to-sample fluctuations. 
This was seen in experiments on mesoscopic electronic samples by Umbach et al. in 
1984 [103]. The theoretical explanation was given by Altshuler [104] and by Lee and 
Stone [102]. For reviews on the subject see Lee, Stone and Fukuyama [110] and the 
book edited by Altshuler et al. [5]. The discovery of the fluctuations boosted the held of 
mesoscopics, as it showed that interference effects, are important in electronic systems, 
even far from the localization transition ! 

The fluctuations are a consequence of scattering from static impurities, and hence 
they are static. Their magnitude is independent of the sample parameters such as the 
mean free path, the sample thickness and, most remarkable, the average conductance. 
Therefore, they are called universal conductance fluctuations. The mean conductance 
in the considered regime comes from multiple scattered diffuse electrons. The UCF are 
a consequence of interference of multiple scattered waves, causing correlations between 

^Although attempts are now made in the group of George Maret in Strasbourg. 
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two diffuse paths. The conductance fluctuates when the phases of the waves in the 
dominant paths change. This happens, of course, if one changes the position of the 
scatterers, e.g. by taking another sample. One also may keep the scatterers flxed but 
apply a magnetic fleld or vary the Fermi energy. (In optical system one can vary the 
frequency of the light.) In all these cases one modifles the phases of scattered waves, so 
that different propagation paths become dominant. The fluctuations are much larger 
then one would obtain classically by modeling the system by a random resistor array, in 
which interference effects are neglected. That approach is valid only on a length scale 
exceeding the phase coherence length, where the fluctuations reduce to their classical 
value. 

As both the electronic and optical systems are in the mesoscopic regime governed by 
the same equation, the scalar wave equation, one expects the fluctuations to be present 
also in our optical systems. Unfortunately, despite all advantages of optical systems, the 
optical analog of the UCF has not yet been observed in optical systems. Such experi- 
ments turn out to be difficult. Although the magnitude of the fluctuations is universal, 
they occur on a background of order (/, where g is the dimensionless conductance (in 
optical experiments one typically has g ~ 10"^). The relative value of the fluctuations 
to the background is thus !/(/, so that the C3 correlation function is of order 
typically of order 10~^. For electrons this problem is absent as moderate values of g 
are achievable. This is also the reason that electrons are easier brought near Anderson 
localization, for which g has to take a critical value of order unity. In the electronic case 
the moderate values of (/, combined with very sensitive techniques for current measure- 
ments, have led to many observations of the universal conductance fluctuations. Recent 
optical experiments suggest, however, that the optical analog of UCF should just be 
experimentally accessible: In chapter 6 we discuss the third cumulant present in the 
total transmission. This quantity is of the same order as the optical UCF, namely 
and was recently measured [105]. We expect that similar techniques can be applied to 
measure the optical UCF. Microwave scattering is also interesting as it combines lower 
values of g with many advantages of optical systems. 

We are not only interested in just the size of the fluctuations, but also in the some- 
what more general frequency correlation. A change in the frequency alters the inter- 
ference pattern, just as occurs by changing the magnetic fleld or the Fermi energy in 
the electronic case. One knows from experiments that using the frequency as tunable 
parameter provides a good way for measuring the fluctuations. In contrast to electronic 
systems, in optical systems both angular resolved and angular integrated measurements 
of the transmission are possible. This led us to introduce three correlation functions, 
Ci, C2, and C3, see chapter 4. We recall that the Ci is the leading correlation in angular 
resolved measurements, the C2 is the main contribution in the total transmission cor- 
relation functions. And the C3, flnally, is the correlation function in the conductance. 
The Ci and C2 correlation were successfully calculated using a diagrammatic technique 
based on the Landauer approach, see chapter 4 and, for instance, Stephen [85]. One 
might hope that the calculation of the C3 or UCF in this approach is also straightfor- 
ward. It is well known, however, that the calculation in the Landauer approach is quite 
cumbersome, since divergencies show up on scales of one mean free path when one treats 
the problem on a macroscopic level using diffusons. 
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To avoid these difficulties, one might be tempted to use the Kubo approach, often 
used in electronic systems to calculate the UCF [104, 102, 110]. Furthermore, the results 
for the conductance obtained by Kubo or Landauer formalism should be identical, see 
Fisher and Lee [111] and Janfien [112]. Yet the Kubo approach cannot be applied 
directly to optical systems, since it is not clear how external lines should replace current 
vertices, and how absorption and internal reflections are to be included. Therefore, 
we use the Landauer approach. Technically, the vertices for partner exchange of two 
diffusons, the Hikami boxes, cause the difficulties in the Landauer approach. Each 
Hikami box brings the square of the internal momentum, whereas the current vertices 
are momentum-independent in the Kubo formula. As a result, the integral over the 
internal momentum of the closed loop is convergent in the Kubo formula, while naively 
divergent in Landauer approach. 

We know of two studies of the C3 in the Landauer approach. In the flrst, Kane, 
Serota, and Lee [86] consider electronic systems and consider local conductance fluctua- 
tions. In contrast to the global conductance where one averages over the whole sample, 
divergencies are present. Kane et al. run in to similar problems as we do, they then 
make elegantly use of current conservation to derive an expression for the correlation 
function. Although in optical systems the conserved quantity is not the intensity but 
the energy, their prediction applies to optical systems as well, since it amounts to a re- 
sult for the same sums of scattering diagrams, involving different parameters only. This 
result has not been conflrmed by a direct derivation, however. Moreover, since it relies 
on a conservation law, it is not clear what happens when absorption is present. The 
second study was done by Berkovits and Feng [90]. After giving a clear discussion of the 
problem, they calculated one of the macroscopic diagrams (presented earlier by Feng, 
Kane, Lee and Stone [91]) and subtracted the divergent parts by hand. In this way the 
correct order of magnitude and the qualitative frequency dependence was obtained. 

It is our goal to clarify the situation by calculating the optical C3 diagrammatically. 
We need a complete analysis of all leading diagrams, in which flnally the divergent parts 
should cancel. We specialize to the case where one measures the conductance. In this 
setup the amplitudes of the incoming and outgoing diffusons are exactly in phase. If the 
C3 correlation is measured as a (small) part of the correlation in the angular transmission 
or total transmission, this phase condition need not be full flUed. Other contributions of 
order that are angular dependent are then present. These contributions are both 
the diagrams presented below, with different decay rates for the incoming and outgoing 
diffusons, but also new diagrams contribute, see Fig. 4.2. In conductance measurements 
these complications do not occur. Our approach immediately allows for inclusion of 
effects due to boundary layers and absorption. Our calculations, although specialized 
to optical systems, are valid for any mesoscopic system. 

The C3 correlation function, deflned in (4.1), involves incoming diffusons -C-'n^, and 
outgoing ones C^^^u^. Due to the factorization of external direction dependence, see (2.68), 
(2.77), and (2.78), it cancels from C3. We can write 



where Q is the [d — l)-dimensional transversal momentum. The function F is the 




Q 
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main object to be determined in this chapter. We thus calculate it at fixed Q and 
with external diffusons being total-fiux diffusons. One finds from (4.1) and (5.1) the 
conductance fiuctuations 

Ct(/€,0) = {T{u)T{u + ^u))-{T{u)){T{u + ^u)) = Y.{T)^,{T),aCs 

abed 

= ^F(g,«,o) 

Q 

= F(0,/€,O) quasi 1-d, (5.2) 

= lU /^F(g,/€,0) quasi 2-d, (5.3) 

= lU2/^F(g,«:,0) 3-d. (5.4) 

In contrast to the Ci and the C2, the C3 depends on the (quasi) dimensionality of the 
system. For electronic systems one finds 

{G{k)G{k + hn)) - {G{k)){G{k + hn)) = (^^^ Ct(0, 0). (5.5) 

These results can be extended for other geometries. If the width is comparable to the 
thickness of the slab, the momentum integral discretizes into a sum over transversal 
eigenmodes [110]. The result can be generalized further to arbitrary geometries by 
taking x^y dependence into account, and calculating the diffusons using appropriate 
boundary conditions. 



5.2 Long range diagrams 

At this point we recall the structure of the leading diagrams for the correlation functions 
defined in (4.1), see again Fig. 4.1. For all diagrams there are two incoming advanced 
fields, which we momentarily term i and j, and two retarded ones, i* and j* . The 
first term on the r.h.s. in Eq. (4.1) follows from the diagram where i and i* pair into 
an incoming diffuson, and the same for j and j*. These diffusons have no common 
scatterers, so for this contribution the expression factorizes into a product of averages. 
For the Ci term in (4.1) such a factorization also takes place [91]. However, in this term 
the pairings are ij* and ji*. In the C2 correlation function there are two terms. In 
the first the incoming diffusons have pairings ij* and ji*. These diffusons interfere in 
some point in space, where they exchange partners through a Hikami-box. The outgoing 
pairings are then u*, jj*. It is this term that contributes in measurements of the total 
transmission. 

In electronic experiments where one measures the conductivity, and in optical exper- 
iments where one uses an integrating sphere on the incoming side for creating a diffuse 
beam, the two amplitudes in the diffuson have to be exactly in phase in order to be lead- 
ing. Therefore, the incoming diffusons cannot have a momentum or frequency difference, 
and the pairing must be ii* and jj*. In Fig. 5.1(a) the incoming diffusons interfere some- 
where in the slab. In a diagrammatic language the diffusons interchange a propagator 
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Figure 5.1: The leading contributions to tlie conductance fluctuations, apart from some 
special short distance processes dealt with in section 5.4. The incoming difFusons from the left 
interfere twice before they go out on the right. The close parallel lines correspond to diffusons; 
the shaded boxes are Hikami vertices; q denotes the free momentum which is to be integrated 
over. 



so that the pairing changes into ij* and ji*. Propagation continues with these diffusons, 
which, due to the different pairing can have non-zero frequency difference and non-zero 
momentum. But to be dominant, the outgoing diffusons can also not have a momentum 
or frequency difference. Therefore, somewhere else in the slab a second interference 
occurs. Again exchanging an amplitude, the original pairing, ii* and jj*, is restored 
and the two diffusons propagate out, see Fig. 5.1(a)i. Some other contributions occur as 
well, yet as the incoming and outgoing pairings are always ii* and jj*: In Fig. 5.1(a)ii 
the first incoming diffuson meets an outgoing diffuson and amplitudes exchange. These 
internal diffusion lines meet at a second point where the original pairings are restored. 
Clearly in this process the intermediate paths are traversed in time reversed order. Due 
to time-reversal symmetry they give a similar contribution as previous Fig. 5.1(a)i. In 
Fig. 5.1(b)i a diffuson breaks up such that one of its amplitudes makes a large detour, 
returns to the breaking point and recombines into an outgoing diffuson. The second 
incoming diffuson crosses this excess path of the amplitude, and one of its amplitudes 
follows the same contour as the first diffuson. The fourth amplitude resides and finally 
recombines with its original partner amplitude to form an outgoing diffuson. Finally, 
Fig. 5.1(c) depicts the situation where only one internal diffuson occurs. Its endpoints 
must lie within a distance of a few mean free paths. Because of its local character this 
class does not show up in the final result; we need it, however, in the regularization 
process, since it contains terms that cancel divergencies from the other two classes. 

Due to the diffusive behavior of the internal propagators we call all these the "long 
range diagrams". Their internal lines are diffusons (ladder diagrams), that is to say, 
in these lines there can be an arbitrary number of scatterers. Note that these long 
range diagrams also include terms with only a few scatterers, e.g. one or two. The 
latter contributions are of course not really long range; however they contribute to the 
geometric series that represents the ladder diagram. We state this explicitly, since below 
we will discuss some unexpected problems of these short ranged contributions to the long 
range diagrams. There are also some special short ranged contributions. Due to various 
subtleties, they resist a general treatment; we postpone their calculation to section 5.4. 
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5.3 Divergencies in the diffusion approximation 

The diagrams for the conductance fluctuations contain a loop; the two internal diffusons 
have a free momentum, over which one has to integrate. In Fig. 5.1(a) we denote 
this momentum q. Physically, one expects important contributions to the conductance 
fluctuations if the distance between the two interferences vertices ranges from the mean 
free path to the sample size. Yet in this section we will show that the q— integral for 
the long range diagrams diverges for large momentum, i.e. when the two interference 
processes are close to each other. The standard picture of diffuse transport with diffusons 
and interference described by Hikami vertices, that works so well for loop-less diagrams 
such as the C2 correlation function (chapter 4), and the third cumulant of the total 
transmission (chapter 6), now becomes spoiled by these divergencies. In section 5.4 we 
solve this problem by going back to mesoscopic scales, and considering all scattering 
events. 

The problem becomes clear if we calculate the diagrams of Fig. 5.1. First, we need 
the expressions for the Hikami vertices. To derive the vertices, we expand the momenta 
to leading order in (q^). For the large q£ this is in principle not allowed, but in practice it 
could still work. The four point vertex and the six point vertex are found by summation 
of the bare vertex and its dressings. We recall from chapter 3 

1 ^ 

H4 = h^i-qi • qs - q2 • q4 + - ^(q- + + i^i)], (5.6) 

i=l 

He = -/i6[qi •q2 + q2 -qs + qs •q4 + q4 -qs + q5 -qe + qe -qi 



+ E(q- + k + ^a-)]. (5.7) 

i=l 



We deflned the pre-factors as 

= aLt^^ = 96^- ^^-^^ 

We denote the momenta of the diffusons attached to these vertices by q,-, where the 
diffusons are numbered clockwise on the vertex and their momenta are directed towards 
the vertex. In the actual calculations we use the Fourier transforms in the 2;-direction of 
the vertices. The additional frequency and absorption terms, according to the diffusion 
equation (2.39), together with the q^ terms lead to a source 5[z — z'). For external 
diffusons, we neglect such terms as they bring contributions of the order Ij L. This ap- 
proximation simplifled the calculation of for instance the long range correlation function, 
see chapter 4. For the internal diffusons, however, the source terms are of leading order 
and cause divergencies. They correspond to the situation where the two interferences 
take place within a distance of a few mean free paths. 

As an example we calculate the diagram presented in Fig. 5.1(a)i. This diagram was 
flrst depicted by Feng, Kane, Lee and Stone [91] and considered in detail by Berkovits 
and Feng [90]. These authors pointed out that a short distance divergency appears. 
For the case of external momenta approximately zero, the Hikami-box (5.6) yields 
if4(q, 0, — q, 0) = 2/^4^^, while the internal diffuson has the form £int(?) = 127r/(£"^g^). 
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Omitting the external lines, the diagram Fig. 5.1(a)i then simply leads to 

J ^HK^,,-^.0)ClM = ^/|^.° = = 0). (5.9, 

This is indeed a cubic divergency in three dimensions, and in general, a (/-dimensional 
divergency in d dimensions. As it is arising from the physically innocent situation where 
the two interference vertices are close to each other, we expect that the divergency has 
to disappear finally. 

We now calculate the diagram for the slab geometry. For simplicity we first consider a 
quasi one-dimensional system in which frequency differences and absorption are absent, 
therefore the decay rate or "mass", Eq. (2.7), vanishes, i.e. M = for all diffusons. 
(Beyond quasi Id one would have to take non-zero M = Q and sum over the allowed 
Q.) From Fig. 5.1(a)i one directly reads off its corresponding expression Fa.i 

F^., = J J dzd/£f„(z)F4(^)/:L(^,^')^4(/)/:L(^). (5.10) 

We label the two incoming diffusons 1 and 3, the outgoing ones 2 and 4, and the internal 
ones C{z5^zr) and C{ze^Z8)^ (with 5, 7 at 2; and 6,8 at z'). The real space expressions 
for the Hikami boxes become 

H,{z) = h,[d.,d.,+\{d.,+d.,f-dl-dll (5.11) 

H,{z') = h,[d.,d., + ^-{d.,+d.,Y-dl-dll (5.12) 

in which is the derivative of the corresponding diffuson; after performing the differ- 
entiation 2^1,3,5,7 should be put equal to 2;, while 2^2,4,6,8 should be put equal z'. Keeping 
-2^2,4,6,8 fixed, we obtain for the 2;-integral after some partial integrations 

/ dzH4{z)Cin{zi)Cin{z3)jCint{z5, ZQ)C;„t{zT, Z^) 

Jo 

^ dz[Cint{z,Z6)S{z-Z8) + Cint{z,Z8)S{z-Z6)]C^z) 



+2h4 j AzC;„t{z,ZQ)C;„t{z,Z8)C'^^{z) 
_ 127r/i4 \r2/ X , 127r/i4 \r-i( \ 

— l~'mt\Z8, Zq)L,-^^{Zs} -\ — L,{nt\Z6, Zs}L^^„{Zq) 

+ 2h4 j dzCint{z,Z6)Cint{z,Zs)C'^n{z). 

Here we also used the diffuson equation, which in this simplified case reads d^&m = and 
d^Cint{z^ z') = 12tt5[z — z') /P . Also carrying out the 2;'-integral we find after performing 
again some partial integrations 



Fa.i = j dz' H4{z')Cout{z2)jC-out{z4) 



127r/i4 ^ , X ^2 



"'^int(^87 ^6)'^in (^8) + 



127r/i4 



jCint{z6, Zs)C'^n{z6) + j dzCint{z, Z6)Cint{z, Zs)C'-^{z) 
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= J dzcuz)cijz) 

+ I 'izCint{z,z) [Cg{z)Cl^^{z) + C'^^{z)C'^^^{z) 

+2/::„(^)/:L(^)A„(^)/:o.t(^)] 

+4hl f dz' f dzCUz,z')Cg{z)C'U^'). (5.13) 

We denote the spatial derivative of C{z) by C . All diffusons are simple linear functions 
in this case, yielding 

F^,, = ^5{0)L + ^. (5.14) 

Note that the pre-factors of the diffusons and the Hikami boxes have canceled precisely. 
This is closely related to the universal character of conductance fluctuations in electronic 
systems, see (5.5). 

The term 5(0) is a linear divergency, which is the cause of all troubles. In the three- 
dimensional case one has to take Q ^ 0. The 5[0)L term will occur also for transversal 
momentum 7^ 0, so that the Q-snm yields the cubic divergency 

W'J^Si0)L = W'LS('H0), (5.15) 

as expected from the above bulk consideration. 

We now give the results of all diagrams of Fig. 5.1. We no longer restrict to the 
M = case. We label the expressions according to the diagrams in flgure, F^, Fj, and 
Fc. In the diagrams of Fig. 5.1(a), the decay rates of Eq.(2.7) for the internal diffusons 
are each other complex conjugate, M and M* . In the diagrams of Fig. 5.1(b) both 
internal diffusons have the same decay rate. Using the deflnition of the Hikami vertices 
and the diffusion equation, we obtain 



FaiM) = —m j dz£f„£L 

+ ^Re j Az A„t(^, z- M) [^C'-lCl,, + 3£f„£ 



12 
out 



+10£:„£' „£i„£out - 4zO£,U' 



-^Ahl j Jdzdz'jC,^t{z,z';M)jC,^t{z,z';M*) x 
HM) = ^^(0) / d.£f„£;„, 

2 



+ 2F /d^A„t(^,^;M) [-2£:XutA„/:out + + -CfXlt 



+ hl J Jdzdz' [cUz,z';M) + CUz,z';M*: 
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c/2 (P , , 

— [Ci^{z)Cout{z)] [Ci^{z')Cout{z')] 



Fc{M) = -^^(0) I dzClCl, 



Re j dz Cint{z, z; M) X 



A;2 

[2'^in'^out'^in'^out+'^in'^out + '^in'^out ~ ('^ + ^^)'^in'^out] (5.16) 

where we used the short hand notation that in the single integrals all incoming and out- 
going diffusons have argument z. To obtain the variance of the fluctuations in quasi-one 
dimension, F's is evaluated at transverse momentum Q = 0. The internal momen- 
tum Q enters the equations via the decay rate of the internal diffusons, deflned by 
= Q'^ -\- k'^ -\- iil. In two and three dimensions a sum or integral over the transversal 
momentum has to be done. 

In Fa and Fj, we can distinguish three contributions. In the flrst term, both Hikami 
boxes operate on the internal diffusons, yielding in the diffusion approximation a delta 
function evaluated in zero. The resulting term is independent of all the momenta of 
the external diffusons. One can see from the diffusion equation that a diffuson decays 
rapidly if its momentum becomes large. Terms of the diffusons with few scatterers are 
dominant at large momentum; they cause the divergence. Our present description of 
these processes is incomplete. To flnd the cancellation of this divergency, calculation 
of the long range diagrams is not sufficient, so that the short distance processes have 
to be examined in detail. This will be done in the next section. The second term of 
(5.16a,b) is a single integral, it comes about when the boxes act on one internal and 
on one external diffuson. This corresponds to the case where one internal diffuson is 
almost empty, while the other diffuson contains an arbitrary number of scatterers. In two 
and three dimensions the momentum integral diverges, since for large Q it behaves as 
J dz J d'^~^QjC[ntiQ', ^1^) ^ I d'^~^Q Q~^. Yet when summing the a, b and c contribution 
this term cancels. The sum gives 



Fa{M) + Fk{M) + F,{M) 

+4hl J Jdzdz'Cint{z,z';M)Cint{z,z';M*) 

X [£:^(.)£l(/) + £:„(.)£:„,(.)£:„(/)£:„,(/ 
+ hl J Jdzdz' [cUz,z'-M)+CUz,z'-M* 



X ^ [Cin{z)Cout{z)] ^ [£i„(/)£out(^')] • (5.17) 

The last two terms involve a double integral describing interference vertices at different 
points in space; technically it arises when both boxes act on external diffusons, or 
from terms where they do so after partial integrations. This term is absent in the 
expression Fc, which contains only one 2;— dependence as one sees from Fig. 5.1(c). 
When performing the integral of F over the transverse momentum the double integral 
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term behaves at large Q as Jdz Jdz' J d'^~^QjC?^^{Q; z, z') ~ / d'^~^Q Q'"^. It is thus 
convergent. One expects that finally only this term will survive. It is also the only 
contribution depending solely on derivatives of the external diffusons. 

In the quasi one-dimensional case where absorption and frequency terms are absent, 
the expressions reduce to 



The second part is a well known result for the UCF in one dimension, but a singular 
part is annoyingly present. Before we can obtain the UCF and correlation functions, 
we have to show its cancellation. We expect the unphysical divergency to disappear by 
summing all diagrams, thus we set out finding all leading diagrams. 

Generating diagrams by computer 

Finding the correct and complete set of diagrams by educated guess proves very difficult, 
especially for the low order diagrams. The set is quite large and not very systematic. 
In the article closely related to this chapter [113], we pursue the technique founded by 
Kadanoff and Baym [114, 115] to generate all diagrams. This method provides a system- 
atic way to construct the diagrams in a particular approximation. The approximation 
is made on the level of a generating functional. The method was very successful for 
our purpose and was almost fully done by dr. Ruud Vlaming. The interested reader 
is referred to Ref. [113]. We verified the set of diagrams generated by the Kadanoff- 
Baym approach using a computer program. We developed this program to generate all 
diagrams with six or fewer scatterers automatically. In the appendix we describe this 
method in detail, in short it does the following: After generating the set, the program 
checks it against double counting. Next, momenta are assigned to all propagators, such 
that momentum conservation is obeyed. It then decides whether a particular diagram 
is of leading order in {k£)~^; this is the case if all propagators can have a momentum 
approximately equal to k. The leading diagrams are expressed in terms of the stan- 
dard integrals defined in section 5.4. We did the subsequent analytical calculation and 
summation of diagrams by hand. The program turns out to be especially useful in 
finding the precise set of leading short-distance diagrams, but also the long distance 
diagrams of section 5.2 were reproduced. It should be mentioned that apart from the 
diagrams from the Kadanoff-Baym approach we also found extra diagrams; in section 
5.4.3 we discuss them and we show that they cancel. Both the set of diagrams from 
the Kadanoff-Baym approach as the computer generated set were calculated in second 
order Born approximation only. Going beyond second order Born would increase the 
number of diagrams dramatically, without bringing different physics. The cancelation 
present in second order Born, will in the end also be present beyond second order Born. 



5.4 Cancellation of short distance contributions 





In this section we show that the strong divergence, that occurs if both Hikami boxes act 
on internal diffusons, cancels. Since a diffuson decays rapidly if its momentum becomes 
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large, terms of the diffusons with few scatterers are dominant when the momentum is 
large; they cause the divergence. We thus show the cancellation by considering in great 
detail the short distance processes. 

It turned out that complications arise for diagrams with less than four scatterers. 
The external diffusons may still contain an arbitrary number of scatterers. For simplicity 
we perform the calculation in this section in the bulk and the external diffusons are not 
attached to the diagrams. Since the divergency is independent of the external momenta, 
the cancellation is generally proven at zero external momentum. For simplicity we may 
then consider an infinite system. As a result all factors can be expressed in the d- 
dimensional momentum q = [Q^qz). We first look at fixed value of the (/-dimensional 
internal momentum q and postpone the integration. We will show that under present 
conditions the integrand is zero for all q, so that there is no divergency after integration. 

In the calculation of the diffuson and the Hikami vertices in the previous sections 
we expanded the integrals in Since we are after contributions for q ~ this 
expansion is not allowed, instead the diagrams can be factorized in products of the 
integrals, defined in section 3.4.1. We still assumed that g ^ A;, as we do not 
need the physics on length scales comparable to the wavelength, but comparable to one 
mean free path only. 

5.4.1 Short distance behavior of long range diagrams 

First, look at the long range diagrams with few scatterings. By writing out the dressing 
of the vertices one obtains the detailed structure of the diagrams, given in Fig. 5.2. In 
the diagrams Fig. 5.1(a)i-(c)i, the incoming external diffusons are connected to 1 and 3. 
The diagrams in Fig. 5.1(a)ii-(c)ii are obtained by connecting external diffusons 1 and 2 
to the incoming side. If such a diagram contains only a few scatterers, it can be grouped 
into one of the classes of the long range diagrams. A simple example is the upper left 
diagram of Fig. 5.2, drawn there with six scatterers. Interchanging the external diffusons 
numbered 2 and 3 clearly leads to another topology. However, with only two scatterers 
present, the topology does not change under this operation, and one must be careful not 
to overcount this term. Bearing this important observation in mind, we sum all long 
range diagrams for an arbitrary number of scatterers. The expression for each diagram 
and its combinatorial factor is given in Table 5.1. It is convenient to collect diagrams 
with an equal number m of scatterers that connect different propagators. Scatterers 
on which one given amplitude scatters twice are thus, momentarily, not counted as 
new scatterers. Looking at Fig. 5.2, m equals the total number of scatterers minus the 
number of scatterers indicated by curved lines. After this re-summation we find the 
important result that each class of "long range" terms vanishes for m not equal to 2. 
Denoting the terms of Fig. 5.2 respectively by -Rfc, and Rc we find Ra = Rh = Rc = 
0, if m > 2. If m = 2, the classes cancel against each other, since we then obtain 
Ra = Rb = -Rc/2 = ^Aj. This imphes 

Ra + Rb + Rc = 0. (5.19) 

We have thus shown that the long range diagrams for an infinite system cancel and thus 
cause no divergency after integration over q. It is essential that we count the degeneracy 
of the low order contributions properly. 
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Figure 5.2: The dressings of tlie sliaded vertices of Fig. 5.1 are written out explicitly, this 
shows the detailed structure of the long range contributions to the conductance fluctuations. 
The close parallel lines, only explicitly shown in the first diagram, are the external diffusons. 
Fig. 5.1(a)i corresponds to the case were 1 and 3 are incoming diffusons; Fig. 5.1(a)ii cor- 
responds to the case were 1 and 2 are incoming diffusons. The vertical, diagonal or curved 
lines linking two dots represent a common scattering of two amplitudes. From top to bottom 
the horizontal lines in a diagram are advanced, retarded, advanced and, retarded propaga- 
tors, respectively. The number of scatterers in the internal and external diffusons is arbitrary. 
Topologically equivalent diagrams are not shown. 
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Table 5.1: Table used for the calculation of the diagrams presented in Fig. 5.2 (for zero 
external momenta). The expression for each diagram is given in the second column. Its 
combinatorial prefactor in the other columns for given number of scatterers. The last column 
counts the number of scatterers that are not included as such in the re-summation. There are 
six diagrams with degeneracy lower than expected: cl for n = 1, also al,c2,c5 for n = 2, and 
c6,c7 for ra = 3. 
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Another, more standard way to verify the important cancellation is the following. 
Recall that we look at a bulk situation with fixed internal momentum q, while we put 
the external momenta zero. Beyond the diffusion approximation an internal diffuson is 
given by Eq.(3.24). We also need the Hikami-boxes beyond the diffusion approximation. 
It holds that 

F4(q, 0, -q, 0) = Illiq) + y + y . (5.20) 

Inserting the values from the table and doing similar but longer calculations for two 
other vertices we find 

F4(0,q,0,-q) = 2F4(q,-q,0,0) = -^Ai(l-Ai), (5.21) 

F6(0,0,q,0,0,-q) = -/-Ai(l - Ai)(l - 3Ai). (5.22) 

Now the long range diagrams of Fig. 5.1 can be evaluated and we insert for the diffuson 
Eq.(3.24, that is dnt = ^{l - Ai)-K 

For the internal part one has at fixed q 

2H4q, 0, -q, OYCUq) + 4H,iq, -q, 0, O fCUq) 
+2F6(q,0,0,-q,0,0)£i„t(g) 

= ^{2A? + AI + 2M1 - 3Ai)} = ^(2Ai - 3A?). (5.23) 

It is essential that the denominators 1 — Ai have disappeared from this expression. 
As mentioned, it means that all high order terms cancel, allowing a cancellation of the 
remainder by low order contributions. From Table 5.1 one sees that we have over-counted 
six types of low order terms. The correction to be subtracted is 

£^ 1 
-{Ai + A2 + As)} 




PAiY -i£^(2Ai + A2)-ifAi if if , , 



8^A;2 / 32^A;3 Snk Snk 32^A;3 



/4^y [ if if -P -if if f \ 

^wJ 1^8^ A; 8^ A; ' 16^F ' ^ ^ 8^A; 8^A; 8^F 

= ^(2Ai-3A?). (5.24) 

It indeed exactly cancels expression (5.23), which arose as a remainder of the long range 
terms. 

From the results of present section one may be tempted to conclude that the long 
range diagrams do not lead to divergencies if the correct degeneracies of the low order 
terms is properly taken into account. Though this conclusion is correct, it is too early 
to draw it, since we show in next parts that some further complications arise. 
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Figure 5.3: Leading diagrams witli topologies not yet contained in the long range diagrams. 
They all contain two scatterers. 




Figure 5.4: Different choices for the loop momenta, all yielding leading contributions; thick 
and thin lines depict different large momenta, both with length k. The first diagram is leading 
for any number of scatterers. When only two scatterers occur also the other two diagrams are 
leading. 



5.4.2 Extra short distance contributions 

Unfortunately, we have not yet finished with the calculation as also other terms are of 
leading order at short distances. First, diagrams with a different topology also occur. 
There are five new diagrams, all with two scatterers, see Fig. 5.3. They can be looked 
upon as diagrams from the set of Fig. 5.2(a), but without any scatterers in the ladders. 

Second, diagrams can have more than one configuration for a resonant arrangement 
of the momenta, i.e. several configurations can be leading. The most important con- 
tribution to the integral (3.21) arises if the amplitude loop momentum p is close to the 
wave-vector, i.e. |p| = A; -|- 0(q). With more loops present it is sometimes possible 
to find more than one choice for the momenta of the amplitude propagators. This we 
illustrated in Fig. 5.4, where the diagram of Fig. 5.2(a)i is repeated. The expression for 
the shown diagram has contribution from three different poles. The first arrangement 
of the momenta is dominant for any number of scatterers, this one is part of the set of 
long range diagrams for any number of scatterers. The two configurations on the right 
give an extra leading contribution. This only occurs with two scatterers present; with 
more scatterers the two diagrams on the right become sub-leading, and only the left 
diagram remains. Apart from this diagram, the same thing occurs for the diagrams of 
Fig. 5.2(c)i and Fig. 5.2(c)ii also with two scatterers, and for the diagram Fig. 5.2(c)vi 
with three scatterers. 

Also the diagrams of Fig. 5.3 have more than one pole: The four left diagrams of 
Fig. 5.3 have two resonant momenta configurations, the right one has three. We sum all 
the contributions thus found (41 in total) and denote them as S„ for n scatterers. They 
are only non-zero for diagrams with two or three scatterers, and read 



= (f) [5(A\'^)^ + 8A\'^/°,^ + 12(/°,^)^ + 8/°:^/^ 



2k 



^jA,A2, (5.25) 



^3 = 8(f)'/°:^/°,f/°:° = |jA,A,. (5.26) 



One sees that these special contributions cancel. 
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Figure 5.5: Leading diagrams witliout amplitude excliange witli one internal difFuson. The 
two drawn diagrams cancel against each other. All diagrams of this class cancel. 






Figure 5.6: Leading diagrams without amplitude exchange with two internal difFusons. The 
two drawn diagrams cancel against each other. All diagrams of this class cancel. 



5.4.3 Interference vertices without partner exchange 

Apart from the diagrams generated by the Kadanoff-Baym approach, we found another 
class of leading long range diagrams when we generated the diagrams by computer. 
In these diagrams the two diffusons have scatterers in common, but no amplitude ex- 
changes. They either have one internal diffuson, see Fig. 5.5 or two, see Fig. 5.6. These 
diagrams are actually of similar type as Fig. 5.1. The diagram with one internal diffuson 
is analogous to Fig. 5.1(c)i, while the ones with two internal diffusons are analogous to 
Figs. 5.1(a)i,(b)i. The equivalents of Figs. 5.1(a)ii,(b)ii and (c)ii also occur. The impor- 
tant difference with the previous diagrams, though, is that at the interference vertices 
no amplitudes exchange. In this respect they differ qualitatively from Hikami-boxes, 
where partner exchange does occur. Explicit calculation shows that these terms cancel 
both at zero and non-zero external momentum, this was also already noted by Kane, 
Serota and Lee [86], see their Fig. 5(c). These classes of diagrams can thus be fully 
neglected. This cancellation is due to time reversal invariance. Apart from these two 
classes, the Kadanoff-Baym approach generates all presented diagrams as we checked 
with the computer program. 

In conclusion, the analysis of this section 5.4 shows that when the external momenta 
vanish, all leading terms cancel at fixed value of the internal momentum. Upon inte- 
gration over the loop momentum one still has zero and, in particular, not a divergent 
contribution. For an infinite system the theory is thus well-behaved at short distances. 
The divergence has canceled in a careful study of the short-distance process; renormal- 
ization as known from field theories was not needed. 



5.5 Application to optical systems 

We have now seen that in infinite systems all divergent terms cancel. In realistic systems, 
such as a slab, the same diagrams describe the relevant physics. They have to be 
evaluated with appropriate diffuson propagators. Generally they can be written as the 
bulk term with additional mirror terms. For a slab we gave them in chapter 2. Knowing 
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that the short distance behavior is regular, we can self-consistently consider all scattering 
diagrams in the diffusion approximation. For the long range diagrams this was done 
already in section 5.3, where we evaluated all long range contributions. The discussion 
of previous section has shown that the only new effect comes from the subtraction of 
Eq. (5.24). In the diffusion limit for a quasi one-dimensional system at fixed transversal 
momentum Q it results in a contact term, labeled F^- 



F,{M) = J dzCliz) ^^^L(-) = {^,m J dzCl{z)CUz). (5.27) 

This indeed cancels exactly the leading divergency that remained in Eq.(5.17) for the 
sum of Fa, Ff, and Fc. The milder divergency present in the individual terms F^, F^^ 
Fc in two and three dimensions canceled already by summing them. Diagrammatically 
this one sees as follows: The diagrams responsible for the latter divergence contain 
one internal diffuson (such in Fig. 5.1(c)). For these diagrams one does not have the 
complications of the previous section; double counting corrections, different resonant 
momentum configurations, and extra diagrams are absent. This was confirmed by our 
computer generated diagrams. The only contributions that remain in (5.17) are thus the 
double 2;— integrals. This finite remainder of our theory yields the sought conductance 
fiuctuations. We thus obtain, with again h4 = £^/(487rA;^), 

F(g, /€, 0) = F,{M) + Fk{M) + F,{M) + F^M) 

dz dz' £int(-2, z'; M)Cint{z, z'; M*) 

X [£:^(.)£i(.') +£:„(.)£:„,(.)£:„(.')£:„,(.')] 

+ |//dzd/[£L(z,/;M) + £L(z,/;M*)]^[A„(z)£o„,(z)] 

xi[£i„(/)£out(^')], (5.28) 



/2 



dz 

in which again = (^^ + + iil. This equation is the central result of the present 
chapter. The upper line of (5.28) corresponds to the diagrams of Fig. 5.1(a), whereas 
the lower line corresponds to the diagrams of Fig. 5.1(b). Note that only derivatives of 
external diffusons are present. As compared to the first line, extra terms are present in 
the second one. According to (2.39) we have 

^ [CUz)Cont{z)] = 2C'Jz)C'^Jz) + 2K^CUz)jCont{z) (5.29) 

The /€— terms are extra terms arising when absorption is present. Finally, with Eq.(5.4) 
the value at vanishing transversal momentum gives the variance of the conductance in 
one dimension, integration over the transversal momentum yields the correlation in two 
and three dimensions. 

Using the general result of (5.28) and (5.4) we consider various cases by inserting 
the diffusons derived in chapter 2. First consider the case of fully transmitting surfaces; 
if we neglect absorption and frequency differences this gives 

32 + 2^2^2 _ 2 cosh 2^1 + ^1 sinh 2^1 ,^ 

^^^^ = 2 guwgl ' ^'-''^ 
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C-r(K,a) 




Figure 5.7: The correlation of tlie conductance as a function of tlie frequency for various 
absorption strengths in one dimension without internal reflection. From upper to lower curve: 
no absorption(K = 0), k = 1/L and, k = 2/L. 



C-r(ic,^2) 




Figure 5.8: The correlation of the conductance as a function of the frequency for various 
absorption strengths in a two dimensional slab without internal reflection. From upper to 
lower curve: k = 0, k = 1/L and, k = 2/L. 
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Figure 5.9: The conductance correlation function versus tlie frequency difference for various 
absorption strengtlis in a tliree dimensional slab without internal reflection. From upper to 
lower curve: k = 0, k = 1/L and, k = 2/L. 



which decays for large Q as Q~ 



H 



ere we recover 



— -0.133, 

3 W W 
4C(3)^^ 0.116^ 

1 



quasi Id 



quasi 2d 



3d, 



(5.31) 
(5.32) 
(5.33) 



in which C is Riemann's zeta function. For cubic samples Lee ct al. [110] find in 2D 
0.186 and in 3D a value of 0.296 (the value in ID is of course unaltered). 

We determine also the frequency dependency of the correlation; this is important as 
it determines the frequency range of the light needed to see the fiuctuations. Taking the 
frequency dependency into account we obtain 



F{Q,u) 



4 {M*' - M' + M'M*L coth M*L - M*'ML coth ML) 



+Re 



2 + 2M'^L'^ - 2 cosh 2ML + ML sinh 2ML 
2M^L^ sinh^ ML ' 



(5.34) 



The correlation decays for large frequency differences as 0^"''/^, as was stated by Lee, 
Stone and Fukuyama [110]. 

We were unable to perform the double integral over the position analytically in 
presence of absorption. In Figs. 5.7, 5.8, and 5.9 we show the one, two, and three 
dimensional correlation functions for various values of the absorption. One sees that 
especially the top of the correlation reduces due to absorption. 

Next we applied our theory to the case of partial reflection at the surfaces of the 
sample. We assume an index of refraction m = J CqI Ci ^ 1. For our purpose the internal 
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0-^(0,^) 




Figure 5.10: The influence of internal reflection on tlie one dimensional frequency correlation. 
With zq = 0, L/10, L/5 (upper, middle, lower curve); no absorption. Also here the fluctuations 
reduce. 



0-^(0,^) 




Figure 5.11: Influence of internal reflection on the frequency correlation in 2D; zq = 
0, L/10, L/5 (upper, middle, lower curve); no absorption. 
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Figure 5.12: Influence of internal reflection on tlie frequency correlation in 3D; zq = 
0,L/10,L/5 (upper, middle, lower curve); no absorption. 

reflections code in only one parameter, the injection depth Zq^ see (2.70), (2.75), (2.41). 
For optically thick samples the correlations are now determined hy Zq/ L. It is physically 
clear that the internal reflections lead to a less steep diffuse intensity in the sample as 
function of the depth. The fluctuations are proportional to the space derivatives and 
thus reduce. We present the results in the Figs. 5.10, 5.11 and 5.12 where we plotted 
the correlation function for various values of the ratio between extrapolation length 
and sample thickness. One sees that already for small ratios of Zq/ L the correlation is 
signiflcantly lower than without surface internal reflection. One sees that neither the 
variance (the value at vanishing 0), nor the form of the correlations are universal, they 
are sensitive to absorption and internal reflections. 

5.6 Discussion and outlook 

The universal conductance fluctuations (UCF) of mesoscopic electronic systems have a 
direct counterpart in other mesoscopic systems with multiply scattered classical waves. 
The corresponding average normalized correlation function is the C3. It is known that a 
naive calculation of this object in the Landauer approach is plagued with short distance 
divergencies. In this chapter we have presented a detailed diagrammatic approach to 
the calculation of the C3 correlation function. We flrst evaluated the leading long range 
diagrams. As expected, each diagram, but also their sum, contains a short distance 
divergency. Also a sub-leading divergency occurs, but for a slab geometry this cancels 
automatically. The study of the cancellation of the leading divergency then was the main 
theme of the chapter. Consistency of the approach requires flnding extra contributions 
that exactly cancel the already determined divergency, and therefore serves as proof 
that the set of long range diagrams exhaust all of them. 

We developed a diagrammatic method that systematically generates all leading scat- 
tering diagrams. This set was checked by computer. Many diagrams had to be consid- 
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ered in detail. We summed these diagrams for an infinite system at fixed value of the 
loop momentum. It was seen that also beyond the diffusion approximation the action 
of Hikami boxes on the long range propagators is to eliminate the long range terms, and 
leave only some low order contributions. Moreover, some low order terms of diagrams 
with diffusons have a lower degeneracy than their higher order equivalents. Taking this 
into account led exactly to a cancellation of all terms in an infinite system. Next we 
discussed that some extra classes of leading diagrams occur, but they all add up to zero. 
Thus in an infinite system all diagrams cancel, so that, in particular, no short range di- 
vergency occurs. All short range contributions could be coded in a contact term (5.27). 
It would be interesting to investigate how one derives it in a non-linear sigma-model 
formulation of the theory. Probably this would go along the lines of Serota, Esposito 
and Ma [116], who introduce so called higher gradient vertices (see section 3.4.1) for 
these terms. 

Subsequently we applied the theory to systems of finite size. Here the short distance 
divergencies cancel as well, because the large scale geometry of the system does not have 
infiuence on short range effects. The final result is non-zero as it describes the correla- 
tion function in terms of derivatives of external diffusion propagators of the geometry 
considered; such terms have no meaning in an infinite system. 

Our central result for the correlation function of the conductance is Eq. (5.28). It is 
obtained by adding the result of the long range diagrams (5.17) and the contact term 
(5.27). When there is no absorption, we find agreement with the result of Kane, Serota 
and Lee [86]. All external diffusons are differentiated once. When absorption is present, 
however, their approach is no longer valid. We found that extra terms appear where 
some external diffusons are differentiated twice or are proportional to k^, see (5.29). 

We have applied the results to realistic optical systems. The frequency dependent C3 
correlation function was calculated for the case where a diffuse incoming beam is used 
and all outgoing intensity is collected. It was seen that both absorption and internal 
reflections decrease the correlations by a considerable amount. This is important for a 
quantitative analysis of experimental data. 

Electromagnetic measurements that involve the C3 correlation have been reported 
by Genack et al. [117], who described the data of their infrared experiments by adding 
the Ci, C2 and C3 contributions, but they incorrectly assumed that the C3 is frequency 
independent. The experimental investigation of optical universal conductance fluctua- 
tions is known to be very difficult. One problem is that if the incoming beam has to be 
diffusive, it h low intensity. 

We propose here a different way to measure the same interference effect. Consider 
a laser beam coming in a given direction a and measure the frequency dependent total 
transmission. Such can be done using an integrating sphere [72]. Then repeat the 
measurement for a very different incoming direction c. Each of these two signals will 
exhibit the large C2 correlation function [72]. However, when the directions a and c 
are not close to each other, the C2 will not contribute to the cross correlate. The cross 
correlation of the total transmission is much smaller than the auto-correlates; it just 
represents the typical term of the UCF in (5.5). As this cross correlate is of relative 
order l/(/^, it is of the same order of magnitude as the third cumulant of the total 
transmission. A very precise measurement of that quantity, see chapter 6, was reported 
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p4-p2-p3 

Figure 5.13: Example of the generation and calculation of a UCF diagram: diagram c5 in 
Fig. 5.2. We follow the placing of the scatterers (top line) and the momenta (lower two lines). 
For two choices of the momenta this diagram yields leading contributions if p3 = — p2 + q or 
P4 = P2 + q (where q ~ £ < p,). 

[105]. It may therefore be expected that it is just possible to measure C3, and thus 
essentially the UCF, with visible light. 

Appendix: Generating diagrams automatically 

In this part we discuss the program to generate diagrams in more detail. We used this 
program to generate to short-distance UCF contributions. We did this calculation with 
zero external momenta. Also absorption and frequency differences were neglected. It 
proved very difficult to produce all diagrams by hand for the short distances. With only 
a few scatterers present, the systematic treatment with Hikami boxes breaks down. On 
one hand some diagrams are over-counted, on the other hand new diagrams appear. 
The goal of the program was to produce all diagrams without double counting. 

Placing scatterers 

The first step in the calculation is to place the scatterers. Starting from the upper left we 
choose the two points on the propagators where we place the scatterer. We allow every 
placing in this step. Recursively we construct the tree of diagrams, until the desired 
number of scatterers is reached. 

We thus have constructed a diagram with n scatterers. First thing we do next, is 
to check the topology. First, it is possible that a certain topology is not allowed. For 
instance the case where the same propagator scatters twice on a scatterer directly after 
one another is already included in the self-energy. Yet also scatterers that can be 'pulled' 
into the external diffuson are not allowed, as they are also already considered. Second, 
diagrams also can be of lower order. An example is the class of diagrams where one 
amplitude scatters on scatterer a), then on b), now on c), then again on a) and finally 
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on b). The sub-leading diagrams were thrown away. 

Next we checked whether the topology was seen before. The recursion tree is very 
degenerate. By its construction, we also get permutations of the placing of scatterers, 
but they should only be counted once. If the topology was seen before, we skipped the 
diagram. After these steps we are left with distinct topologies that can not be excluded 
solely on their topology. This set is quite large and is a bottleneck in computation: some 
500,000 diagrams for 6 scatterers. 

Placing momenta 

The other important reduction of the number of leading terms asks for the explicit eval- 
uation of the diagrams. We place momenta in the diagram and calculate the diagram. 
First we fix the momentum in the left part of the uppermost propagator. Because of 
momentum conservation and the fact that the outgoing diffusons carry no momentum, 
the momentum in the left part of the lowest propagator is the same. Now we repeatedly 
use momentum conservation at the external legs of the diagram and at the scatter- 
ers. If after the use of momentum conservation there are still undetermined momenta, 
we introduce a new independent momentum. We repeat this procedure until no mo- 
menta are left undetermined, see Fig. 5.13. Although at first sight this procedure looks 
straightforward, it is possible that, if the new momenta are put randomly, momentum 
conservation is broken in the end. By carefully choosing the order in which we introduce 
new momenta, it is possible to conserve momentum everywhere. 

Now we know all momenta. To have a resonant, i.e. leading, diagram the momenta 
of all propagators need to be around |p| k. It is possible that this is already the case 
after the last step. But otherwise one momentum has to be made small as compared to 
I A; I; this is the intensity loop momentum q. It turned that there are maximally three 
possibilities for the choice. In Fig. 5.13 the two possible choices are p3 = —p2 + q and 
P4 = P2 + c[- For each momentum configuration thus found, we counted the number of 
independent of momenta left. This yielded the expression for the diagram in terms of 
the / integrals of section 5.4 and 3.4.1. If either the number of retarded of advanced 
propagators is zero, the integral has no poles and becomes sub-leading. 

On an IBM RS6000 work station the maximal number of scatterers in the diagrams 
was limited to six. This proved large enough, as already with four or more scatterers 
the special features of the short distance diagrams disappeared, as we saw in section 
5.4. Thus we were able to check this up to two orders higher. 
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Third cumulant of the total 
transmission 

Uptill now we have considered correlations between intensities defined as 

{TabTcd) 



(Tab) (Ted) 



1 + Ci + C2 + C3. (6.1) 



The Ci((5,0), C2((5,0), and C3(0) correlation functions thus describe correlations be- 
tween two intensities only. By their definition the value at zero frequency and zero 
momentum-difference, i.e. the peak value, presents the second cumulant of the distri- 
bution functions of the transmission quantities. This interpretation of the correlation 
functions and the interest in higher moments is rather recent. It is especially of interest 
near the Anderson transition, where interference gets more and more important (more 
Hikami boxes "appear") changing the distribution function. In chapter 7 we calculate 
the full distribution function for the angular resolved transmission and the total trans- 
mission in the regime of moderate g. Yet first, in this chapter we calculate the third 
cumulant of the total transmission. Using the good old diagrammatic approach we 
relate the third cumulant normalized to the average, {{T^)) to the normalized second 
cumulant {{T^)). For a broad Gaussian beam profile we find {{T^)) = ^((^^a which 
is in good agreement with experimental data. 



6.1 Introduction 

The interest in the full distribution functions in mesoscopic systems is rather recent. 
Examples of such distributions are the intensity distribution in speckle patterns for 
classical waves [95, 118] and the conductance distribution for electronic systems [23]. 
The size of the fiuctuations and the shape of the distribution relates to the 'distance' 
from the localization transition. Far from localization, diffusion channels are almost 
uncorrelated and fluctuations are small (except the optical speckle pattern in the angular 
resolved transmission). The correlation between the channels increases if the localization 
is approached. The relevant parameter is the inverse dimensionless conductance l/(/, 
which can be interpreted as the chance that two channels interfere. We express the 
dimensionless conductance in the thickness of the sample L, the mean free path £, and 
the number of channels 

4Ni , , 
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We calculate the number of channels in analogy with a waveguide, where it is unam- 
biguously defined. In the diffuse mesoscopic regime is a small parameter of our 
perturbation theory; experimentally this proves fully justified, as there one has g ~ 10"^, 
[105]. Close to Anderson localization g approaches unity, and fiuctuations increase. The 
central question is how the distributions change as we approach the strong localization 
regime [23, 119]. 

The total transmission is a constant superposed with fiuctuations. In first order of 
g~^ the fiuctuations have a Gaussian distribution [72, 118]. The relative variance of this 
distribution, the top of the C2, is proportional to (/~^, it is thus a factor g larger than for 
the conductance fiuctuations. This sensitivity of the total transmission to interference 
processes and its simple limiting behavior (as compared to the angular resolved trans- 
mission) make it an ideal quantity to study mesoscopic transmission. Recently the third 
cumulant of the distribution was found experimentally by De Boer et al. [105]. In this 
chapter we present the theoretical details of that work. We focus on the Gaussian distri- 
bution and the deviation from the Gaussian due to the presence of the third cumulant. 
The structure of this chapter is as follows. In section 6.3 we discuss the character of the 
probability distribution. In section 6.3.2 and 6.4.1 we calculate the third cumulant of 
the probability distribution. Next we calculate experimental corrections to our result in 
section 6.4.2 and 4.3.2, after which we compare our results with the experimental data 
in section 6.5. 

6.2 Intermezzo: Cumulants and moments 

In the following we need some generalities from probability theory, see e.g. Van Kam- 
pen [120]. To a general probability distribution P{x) one can assign the characteristic 
function </)(A;), defined as 



where (x") is the n*'* moment of the distribution: JdxP{x)x" . From this characteristic 
function one can again reconstruct P[x). In statistical problems one often encounters 
the addition of many independent chances. The central limit theorem tells us that the 
sum will tend to a Gaussian distribution (no matter the distribution of the individual 
variables). Moments are a bit unapt for Gaussian distributions, as each moment is non- 
zero. Cumulants are more appropriate here. The cumulants K„ are implicitly defined 
as 



Expansion with respect to k yields the relation between cumulants and moments 






(6.4) 



n 




(6.5) 



j=i 
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Thus Ki = (x), K2 = {x'^) — and K3 = (x^) — 3{x'^){x) + 2{x)^. For a Gaussian 

distribution 

P{x) = —=e-^t^ (6.6) 
(7^/2^T 

only the first two cumulants are non-zero: Ki = (x) and K2 = c^, where cr^ is the 
variance. 



6.3 Cumulants of the probability distribution 

In this section we introduce the probability distribution of the total transmission of 
scalar waves and we discuss some of its properties. The corrections for vector waves will 
be made in section 4.3.2. We link the moments of the distribution to diagrams. The 
moments of the probability distribution P{Ta) can be extracted as 

{T',) = jdT,P{T,)Tl (6.7) 

In a diagrammatic approach the A;-th moment can be represented by a diagram with 
k diffusons on both incoming and outgoing side. The k = 1 term is the average total 
transmission (Ta), as given by the Schwarzschild-Milne equation in Eq.(2.80). This 
quantity is a single diffuson and is thus independent of channel-to-channel correlations. 
The second moment can be decomposed in the first two cumulants: 



the double brackets denote cumulants normalized to the average. Diagrammatically 
we depict the second moment in Fig. 4.6. The decomposition in cumulants proves 
useful as each cumulant corresponds to a different number of interactions between the 
diffusons. In the first term, Fig. 4.6(a) there is no interference; it factorizes in the 
average transmission squared (apart from a small correction discussed in section 4.3.2). 
The second term, Fig. 4.6(b), is the second cumulant {{T^)). It gives the variance of 
the fiuctuations. Interactions between two diffusons are responsible for the presence 
of this second cumulant. This is just a special case of the C2 correlation function 
studied in chapter 4. The C2 frequency correlation function was defined: C2(Ac<;) = 
{{Ta{io)Ta{io -\- Aid))). FoT OUT casc find that {{T^)) = 6*2(0) and thus corresponds to 
the peak value of this correlation function. 

Similar to the second moment, one can distinguish three different contributions to 
the third moment, 

^ = l + 3((r„^)) + (ra). (6.9) 

We have drawn the corresponding leading diagrams in Fig. 6.1. The first term. Fig. 
6.1(a), again corresponds to the transmission without interference. The second term. 
Fig. 6.1(b), is reducible in a single diffuson and a second cumulant diagram. From the 
figure it is clear that this decomposition can be done in three ways which refiects in the 
prefactor of {{T^)) in Eq. (6.9). The third contribution stands for the third cumulant 
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a) 



b) 



c) 



d) 



Figure 6.1: The three contributions to the third moment of the total transmission. Diagram 
a) corresponds to independent transmission channels; it is of order 1. In diagram b) there 
is correlation but it can be decomposed into the second cumulant; this is of order g~^. The 
diagrams c) and d) are the contributions to the third cumulant, 0{g~^). 



in the distribution and expresses the leading deviation from the Gaussian distribution. 
This is the term we are mainly interested in. It consists of two related diagrams: Fig. 
6.1(c+d). The three intensities can interfere twice two by two, or the intensities can 
interact all three together, with a so called Hikami six-point vertex. Both contributions 
will prove to be of the same order of magnitude. The strength of the effect can be easily 
estimated using the interpretation of 1/g as an interaction probability. By looking at 
the diagram, the third cumulant is proportional to the chance of two diffusons meeting 
twice, thus of the order We find the basic result 

{{T^)) oc {{T^)r. (6.10) 

The rest of the chapter essentially consist of proving this relation and determination of 
the prefactor. Finally we will compare this relation to experimental data of Ref. [105]. 

Because one samples only a finite number of channels in an experiment, the law of 
large numbers predicts a distribution with some non-zero width even if we only would 
consider disconnected diagrams. However, we have shown in section 4.3.2 that this effect 
brings only a negligible contribution to the measured second cumulant; this is also true 
for the third cumulant as we will show below. The large fluctuations are almost only 
due to the interference. 

Note that in our calculations we again neglect diagrams with loops. An example of 
a loop diagram is the C3 contribution to the second cumulant, see chapter 5. The C3 
diagram contains two Hikami boxes, therefore it gives a contribution of order to 
the second cumulant. In general, one easily sees that to create a loop, one needs more 
interference vertices. Therefore, these diagrams are of higher order in l/(/, and we did 
not calculate them. The leading contributions to the cumulants are by far sufficient for 
the description of the experiment of Ref. [105]. 
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6.3.1 The second cumulant 

First we recover the results for the second cumulant as calculated in chapter 4. In 
the plane wave limit of the incoming beam, all transverse momenta are absent and the 
diffusons are simple linear functions. One obtains for the second cumulant 



-4U 



Taking only this second cumulant we find for the moment a Gaussian probability dis- 
tribution in the plane wave limit [118, 72] 




(6.12) 



We also recall the infiuence of the beam profile on the correlation. There the Q- 
dependence of the incoming diffusons should be taken into account. One finds the 
total cumulant by integrating over the momentum with the corresponding weight. For 
a Gaussian beam profile we found 

{{T'J) = ^J<i'Q e-''oQy^F,{QL), (6.13) 

with F2{x) = [sinh(2a;) — 2a;]/[2a; sinh^ x]. If the incoming beam is again very broad, 
po ^ L, one recovers the plane wave behavior {{T^)) = 2/3g. Note that this agreement 
is found by identifying the area of a Gaussian profile with A = ttPq. This definition is 
somewhat arbitrary and other choices are also possible. After fixing this definition, no 
freedom remains and we will see that for the third cumulant a Gaussian profile leads to 
other results than a plane wave. 



6.3.2 The third cumulant 

We now discuss the calculation the third cumulant. As mentioned, there are two pro- 
cesses contributing. One with two four point vertices, that we term {{T^))c and one with 
a six point vertex ((r^))^, where we have chosen the subscript according to Fig. 6.1. 
These processes were, to our knowledge, not studied before in detail. Our calculation 
follows the lines of the second cumulant calculation. 



Interference via two four point vertices 

First consider the diagram in Fig. 6.1(c). We have labeled to incoming diffusons with 
odd numbers, the outgoing ones with even numbers. Two incoming diffusons, Ci and 
C3 meet at a position z. In a Hikami box the diffusons Ci and C3 interfere into C2 
and an internal diffuson C™g. The C2 propagates out, whereas C™g interferes again 
at z' with incoming diffuson £5 into two outgoing ones C4 and Cq. Apart from this 
process, also three other sequences of interference are possible. This means that the 
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diffusons can also be permuted as: (£i, £3, £5, £2, -^4, ^e) — ^ (£3, £5, £1, £4, ^e, £2) — ^ 
(£5, £1, £3, £5, £2, £4)- We denote the sum over these permutations as J2per- As the 
diagrams also can be complex conjugated, there is also a combinatorial factor 2 for all 
diagrams. (Note that complex conjugation for the second cumulant diagram, yields no 
different diagram and thus it should not be counted.) The expression for the diagram 
of Fig. 6.1(c) is now 



((r„3))e = {T^)-''2Y,A j"^ d.z j"^ d.z' H,{z)H,{z')Cr{z)C2{z)Cs{z)x 

per -^0 -^0 

C,{z%,{z')Ce{z')C^^{z,z'). (6.14) 



In turns out that it is useful to rewrite the form of the Hikami boxes as introduced in 
Eq.(3.5) into an equivalent expression, using momentum conservation qa+c[a'+c[fc+c[fc' = 
0. In real space the use of momentum conservation corresponds to partial integration. 
The Hikami box is again simplified using the fact that there are no transversal momen- 
tum terms, or Q terms, for the outgoing diffusons. Also we neglect the source terms 
of incoming and outgoing diffusons. Using the numbering in Fig. 6.1, we obtain 

H,{z) = ^[25^1 + 2a.,a,3], H,{z') = ^[2a,,a,, - dl + qH (6.15) 

Source terms, i.e. g^^— terms of the incoming and outgoing diffusons were again neglected, 
but the source term of the diffuson between the vertices is important. As one sees with 
the diffusion equation (2.39), it brings 

- dlC-^{z, z') + QlC-^{z, z') = ^5{z - z'). (6.16) 
The contribution from the source term, i.e. H4[z') oc —d^ ''fQl^ H4[z) oc dz^dz^ +5^2 5; 



)2 I /02 

is 



-4A y^^j E / Md.A.+dz,dz,)C,C,C, j dz'{-dl+Ql)C}-'C,C,C, 
—f A 

= dz[dz,dz^+dz^dz,+ ... + dz,dz,]Ci...CQ. (6.17) 

4Q7r/e^ Jo 

Although this corresponds to a local process , just one 2;— coordinate is involved, it is of 
leading order. Together with the expression coming from H4[z') proportional to dz^dz^^ 
we find for the total contribution of the process in Fig. 6. 1(c) 

/' dz' C',{z')C,,(z')C',(z')aX'"'{z, z') 

Jo 

- 4g^fe4^r ^3 X [^-1 + dz, dz, + dz, dz, + dz, dz, + dz, dz, + dz, dz, ] 
xCi{z)C2{z)C3{z)C4z)C5{z)Ce{z), (6.18) 
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where C'{z) denotes the derivative towards z of jC{z). Calculated for a plane wave it 
gives 

22 



((T!)). = (6.19) 
which is indeed proportional to (/~^, as predicted. 



Contribution of the six point vertex 

There is another diagram contributing to the third cumulant which is of the same order 
of magnitude as the process calculated above; it is depicted in Fig. 6.1(d). The six 
sided polygon is again the Hikami six point vertex ife, see chapter 3. It can be thought 
of the following way: The use of Hikami box in the previous section assumes that the 
outgoing legs scatter at least once before they propagate out or interfere again. This is 
a reasonable assumption for the outgoing diffusons, but for the internal diffuson C™g it 
also possible that coming from z it directly, i.e. without scattering, interferes again at 
z' . This process is not included in the calculation of the previous section but has to be 
studied separately. The unscattered intensity decays exponentially over one mean free 
path, therefore this process is only important if z and z' are within one mean free path, 
to this diagram. It is not allowed to dress the bare six-point vertex (leftmost r.h.s. 
diagram in Fig. 3.2) with a scatterer that connects two opposite propagators. This 
dressing gives also a leading contribution even if the dressing done with an arbitrary 
number of scatterers, but the resulting diagram is the same as the composed diagram 
with two four point vertices and should thus not be counted. Yet this observation is 
useful to check the combinatorial ratio between the six point vertex and the composed 
diagram: the forbidden dressing can be performed in three ways. As the diagrams can 
also be complex conjugated, there is also a factor 2 for all diagrams. 
In the lowest order of (qi) we found for the six point vertex 

-f 

^6 = 7^^^[qi-q2 + q2-q3 + q3-q4 + q4-q5 + qs-qe + qe-qi + (6-20) 

After a Fourier-transformation in the 2;— direction the six-point vertex yields a contri- 
bution to the third cumulant 

((r„3)), = j\z[d.,d.,+d.,d.,+d.A,+d.A.+d^A.+dA] 

487rA;^ Jo 

xCi{z)C2{z)C^{z)C^{z)C^{z)Ce{z). (6.21) 

Here we used the fact that all outgoing diffusons have zero transversal momentum. 
Therefore all QiQj-ievm.s are absent. In the limit of an incoming plane wave we find a 
contribution to the third cumulant 




(6.22) 



The contribution from the source term, i.e. Eq.(6.17) of the previous section, exactly 
cancels the contribution from the six-point vertex. The cancellation seems plausible as 
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one does not expect short distances properties to be important in the total process. 
Nevertheless, this cancellation depends on the precise form of the Hikami four-point 
vertex in Eq.(6.15). If we use other equivalent forms of the Hikami-box the contributions 
of the single and double integral in Eq.(6.18) shift with respect to each other and a full 
cancellation is not present. Of course, neither the result for Eq.(6.18) nor the final result 
for {{T^)) relies on this choice. The precise mechanism behind this is not clear to us. 
However, using the cancellation we only need to consider the term in Eq.(6.18), which 
comes from H4[z) oc dz^dz^ + dz^dz^, H4^{z') oc dz^dz^ and the permutations. We thus 
obtain for the third cumulant 

r dz' C\{z')U(z')C'„(z')d,C""U, z'). (6.23) 

JO 

In the next section we calculate this expression for several incoming beam profiles. 



6.4 Corrections for the experimental situation 
6.4.1 Influence of incoming beam profile 

Now that we know the leading interference processes, inserting the diffusons gives the 
final value of the third cumulant. We first consider the simple case of incoming plane 
waves. As there can be no transversal momentum difference in the incoming amplitudes, 
all Qi vanish. As a result all diffusons are simple linear functions of z. And we find 
from Eq. (6.23) 

{{T'a)) = = f ((^.'))' Pl^^e ™' ^0 » L. (6.24) 

In practice, however, we deal with a Gaussian beam with limited spot size, infiuencing 
the cumulants in two ways. First, if the spot size decreases to values comparable to the 
sample thickness we have to convolute over a range of incoming momenta, just like we 
did when calculating the second cumulant. Second, the Gaussian profile brings an extra 
geometrical factor as we will show below. 

We need the expression when diffusons with arbitrary momentum are connected 
to the Hikami boxes. Because of momentum conservation and phase condition on the 
outgoing diffusons, the transversal momentum Q-j of the diffuson connecting the two 
four-boxes must equal Q^. The integration over the possible momenta results again in 
a Gaussian weight function. From the definition (4.11) we derive 



j d^Pid^Psd^Ps ci>{P,)r{P, + Q,) cf>{Ps)cf>*{Ps + Qs) cf>{P5)4>*{P5 + Q,) 
= e^^[-pl{Ql + Ql + Ql)l^] (6.25) 
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We use momentum conservation to eliminate also and reduce the integration to two 
transversal momenta. The final result for the third cumulant is obtained by inserting 
the momentum dependent diffusons into Eq.(6.23). This gives 

(6.26) 



wi 



ith 



^ (a;i + a;3)^a;5 cosh(a;i + a;3) (xi - a;3)^a;5 cosh(a;i - Xs) 

F3{xi,X3,X5) = 



per 



{Xi + X3 + X5y{xi + X3 - X^y {Xi - X3 + X5y{xi - X3 - X5) 

(xi + a;3)a;5 sinh(a;i + 3:3) ^ (xi — a;3)a;5 sinh(a;i — 3:3) 



(a;i + a;3 + a;5)(a;i + a;3 - 3:5) {xi - X3 + X5){xi - X3 - X5) 
{xi-\-X3) cosh(a;i + a;3 + 2a;5) (3^1 + 3^3) cosh(a;i + a;3 — 2x5) 



4:{xi+X3 + X5y 4:{xi + X3-X5 



2 



[xi — X3) cosh(a;i — a;3 + 2a;5) ^ [xi — X3) cosh(a;i — X3 — 2xz 



A{xi-X3 + x^y A{xi-X3-x^y 
1^X5 sinh(a;i) sinh(a;3) sinh^(a;5)j . (6.27) 

These two equations are the main result in this chapter. We study again the behavior 
if the beam diameters are wide. In the limit of large beam diameter (po ^ L) one finds 
^3(0, 0,0) = yI, this means for the third cumulant ((T^)) = 4F3(0, 0, 0)/3^2^ or 

16 

((r„')) = — ((r„2))2, Gaussian profile; po > i^, (6.28) 


which differs by a factor | from the plane wave limit Eq.(6.24). This is purely a geomet- 
rical effect, depending on the profile of the incoming beam. In a real space picture one 
understands this effect best. The correlation depends on the distance: it is strongest 
if the incoming intensities are close. Therefore it is not surprising to see the influence 
of the overlap. In next chapter we calculate this geometrical factor for higher orders 
also (the area of a Gaussian beam is deflned differently there). For the experimental 
relevant case that the beam diameter is roughly equal to the thickness, we calculated 
Eqs.(6.26) and (6.27) numerically. It then turns out that the behavior of Eq.(6.28) is 
found for a large range of beam diameters. The increase of the correlation for smaller 
beams, turns out to be roughly the same for both the third cumulant and the second 
cumulant squared. All corrections to (6.28) turn out to be relatively small, as we will 
discuss below. Apart from this advantage, errors in the sample thickness and the mean 
free path cancel by presenting the results as the ratio between the second cumulant 
squared and the third cumulant. 



6.4.2 Influence of internal reflection 

In this section we calculate the influence of internal reflection on our results. We saw in 
chapter 4 that surface reflection decreases the C2 correlation. In Eq.(6.28) corrections 
from boundary reflection cancel partly. We did not calculate the influence of internal 
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reflections for the general case, but only for the case of very broad beams (i.e. only for 
Q independent diffusons). One expects that this behavior may be extrapolated to the 
(5— dependent case. At least, for the second cumulant this is a good approximation [68]. 
The Q— independent diffusons in the presence of internal reflections are given in chapter 
2, Eqs.(2.60,2.75). If internal reflections are absent Zq equals 0.71£ and the corrections, 
which are of order Zoj are often negligible. With internal reflection present Zq increases 
and should be taken into account. In flrst order of Zoj L the second and third cumulant 
behave as 

^2\\ 2 / Zq\ iirrSw 16 15 Zq 



Therefore, the central relation (6.28) h correction 

{{ri'» 16 / 3 z„ 



The experimental determination of the index of refraction of the sample, which gives Zq, 
is difficult [64]. Fortunately, the correction is rather small for the experimental situation 
considered. 



6.4.3 Contributions from disconnected diagrams 

So far we calculated the leading contributions to the second and third cumulant. They 
are the connected diagrams in Fig. 4.6(b) and Figs. 6.1(c+d), respectively. Yet there are 
also contributions to the second and third cumulant from disconnected diagrams. As 
we saw in chapter 4, the diagram in Fig. 4.6(a) gives an additional contribution to the 
second cumulant, and likewise, the diagrams of Figs. 6.1(a) and (b) give a contribution to 
the third cumulant. These disconnected diagrams correspond to cumulant contributions 
that are not (fully) due to interference. They describe effects that have little to do with 
the interference effects we are after. Here we calculate their contribution and show that 
they are small. ^ 

The contribution of the disconnected diagram Fig. 4.6(b) to the second cumulant 
was calculated in chapter 4 

^pr^d^gm _ 1 
((^. - — ^v^) — = -N' ^^-^^^ 

In analogy to a waveguide, this result describes irreducible contributions from discon- 
nected diagrams. It can be interpreted as the inverse of the number of independent 
speckle spots in transmission at the exit interface [106, 105]. 

Extra contributions to the third cumulant 

We apply the same method from chapter 4 for contributions to the third cumulant. 
Following the waveguide argument as above, we find the contributions to the third 
cumulant of the diagrams of Figs. 6.1(c-|-d), 6.1(b), and 6.1(a), respectively, 

(TO) = + + (^-32) 



^This calculation was initiated by Drs. Johannes de Boer for our joint articles Refs. [105, 73] 
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The first r.h.s. term is tlie connected diagram. Clearly the second term, the diagram 
Fig. 6.1(b), gives a much larger contribution to the third cumulant than the Fig. 6.1(a), 
as the diagram Fig. 6.1(b) is already enhanced by some interference. In the following 
we only consider this diagram. 

As one sees from Fig. 6.1(b) there are three possibilities to combine the three diffusons 
into two connected diffusons. and a single diffuson. Attaching outgoing directions to 
the amplitudes at the exit interface gives (see Fig. (6.2), 

*(i?i)e'^i^i **(i?i)e-'^^^i *(i?2)e'^''^^**(i?2)e-'^^^^*(i?3)e'^^^^**(i?3)e-'^'=^^ (6.33) 

Obviously there are six possibilities to pair the outgoing directions into intensities. In- 
tegrating over the transversal coordinates i?i,i?2 and R3 gives the following six contri- 
butions, 

/con(0)/con(0)/(0) 
/co„(Pl - P3)/co„(i'3 - ^0/(0) 
/con(0)/con(i'l " P5)I{P5 " ^1 ) 
/con(0)/con(i'3 " P5)I{P5 " ^3) 
/co„(Pl - P,)Icon{P3 - Pl)I{P, - Ps) 

/co„(Pl - P3)Icon{P3 - P,)I{P, - Pi), (6.34) 

where Icon denotes the transmitted intensity coming from the connected part of diagram 
Fig. 6.1(b). The third cumulant is in the discrete mode model given by 



{TaY 



{TabTab'Tab") " ^'^{Tab) '^{TabTab') + '^'^{Tab)" 



b,b'b" b b,b' 



(6.35) 



Similarly to the second cumulant we insert all pairings of Eq. (6.34) and find for the 
contribution from disconnected diagrams 

((^a ))dis = 77^ / d'Q,d'Q2d'Qs [/co„(0)/co„(gi)/(-gi) + Icon{Ql)Icon{Q2)I{Q3)] • 

{Tar J 

(6.36) 

The integrand is dominated by its first term, which is depicted in Fig. 6.2. We explicitly 
evaluated this term. The calculation follows completely the line of the second cumu- 
lant calculation, yet it is slightly more complicated as it has a f^— dependent outgoing 
diffuson. 

The outgoing diffusons have, from top to bottom in Fig. 6.2, momenta 0, Qi, and 
— Qi, respectively. The incoming diffusons therefor have momenta Qi — and —Qi, 
respectively, where Q is an extra free momentum. The weight factor for the incoming 
diffusons is 

exp[-|(g? + g2 + (gi-g)2)] (6.37) 

For fixed momenta the connected part of the diagram is 

f^^dzCiz; Q)C3{z; Q, - Q)H4z)C2{z; 0)C,{z; Q,) 

{J-a) -Jo 
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Figure 6.2: The leading disconnected contribution to the third cumulant, the box symbolizes 
again the Hikami vertex. 



where we have defined x = \Q\L^ Xi = \Qi\L^ and y = \Qi — Q\L. As for the connected 
diagram there is again a weight factor from the Gaussian incoming beam profile. The 
integration over the free momenta Q2 and Q3 in Eq. (6.36) gives a factor [irk'^Y . Com- 
paring the result thus found for Eq. (6.36) with the last two terms in third cumulant in 
the waveguide model, Eq. (6.32), we define N* as 



The number N* is inversely proportional to the contribution of disconnected diagrams 
to the third cumulant for a Gaussian profile. 

Polarization effects 

We did not take the vector character of the light into account yet. The two independent 
polarizations of each outgoing direction effectively double the number of independent 
speckle spots N' . Thus for an incoming plane wave of unit intensity with fixed polar- 
ization the total transmission and the conductance g are twice as large as they would 
be in the scalar case. As we work with normalized cumulants, this effect reduces the 
value of the second cumulant (oc !/(/) and the value of the third cumulant (oc 
Therefore, one immediately sees that the relation {{T^)) oc ((T^^))^ (Eq. 6.28) is not 
affected. The vector character does reduce the correction Eq. (4.22) by a factor 2. For 
the experimental data of Ref. [105] we list the number of modes N' as well as N* in 
Table 6.1 (including the doubling). The difference between N' and N* is small. Yet due 
to the diffuse broadening of the beam, N' is much larger than the number of modes in 
the waveguide N = k'^AjAir (typically some 30 times in the experiment). 

Summarizing the previous sections we have included three corrections. We first 
obtained the result for very broad Gaussian beams, in the large Ljl limit, Eq. (6.28). 
The first correction was the infiuence of a finite beam diameter, it changes the diffuse 
intensity from linear into an exponential decaying, see Eq. (2.42). This correction is 
contained in Eq. (6.27) and Eq. (6.26). The presence of internal reflections also changes 



— [Ax^xy sinh Xi -\- Xi[x^ -\- y^ — xl{x -\- y) — xxi{x -\- y)] sinh(a; -|- y) 

-\-xi[—x^ + y^ — xl^—x + y) — xxi{—x -\- y)] sinh(a; — y)] x 
\[x'^ — 2x\{x'^ -\- y^) -\- x'^ -\- X2 -\- 2x^y^) sinh x sinh Xi sinh yj , 




(6.39) 
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the spatial dependence of the diffuse intensity resulting in a correction Eq. (6.30). The 
third correction is of another nature, it is the only process that does not come from 
interference, but from disconnected diagrams. Only this term depends the number of 
modes, which in the vector case is twice as large as in the scalar case. 

6.5 Comparison with experiments 

We reproduce the data set found experimentally in Ref. [105] in Table 6.1 and Fig. 6.3. 
The experiments reported there were performed with seven different samples. The 
experimental setup and measurement technique used is extensively described in Ref. 
[72]. Samples consisted of 36 vol.% rutile Ti02 pigment on a transparent substrate. 
The extrapolation length was estimated from the effective index of refraction to be 
Zq l.lfim. The absorption length ia was determined to be ~ 70/im. Different values 
of the conductance g were probed by taking several sample thicknesses and by varying 
the beam diameter. The fluctuations in the total transmission were measured by varying 
the wavelength of the light. 

For a very broad beam we found the simple relation (6.28) between second and third 
cumulant. A weighted least square flt to 



of the raw experimental data yields a prefactor 2.9 it 0.6. However, as discussed above, 
there are three corrections to be made. First, if the beam width becomes comparable to 
sample thickness we have to perform the integrals (6.13) and (6.26): If the beam width 
reduces, g decreases accordingly and both cumulants increase in absolute size. Yet the 
precise increase is somewhat different resulting in a somewhat smaller prefactor in Eq. 
(6.28). We corrected each data point individually for its flnite focus, mapping it to the 
inflnite focus case. The third cumulant was multiplied by a factor that ranged from 1.03 
to 1.13, as L/po ranged from 0.41 to 6.3, see Table 6.1. This is the largest correction; it 
changes the prefactor some 10%. Secondly, we corrected for internal reflections according 
to Eq.(6.30). The third correction comes from the disconnect diagrams. After all these 
corrections the data should again obey the law: {{T^)) = 3.2 ((T^^))^. We have plotted 
the results in Fig. 6.3, where the points are the corrected data points and the line is the 
theoretical prediction. A least square flt gives 



Note that there is no adjustable parameter. We flnd that there is good agreement 
between experiment and theory. All corrections are minor as compared to the error in 
the data and the data can even be described rather well disregarding the corrections. 
Inspecting the flgure one might be tempted to use a linear flt, but in Ref. [105] it was 
shown that this flt is statistically very improbable. 



const.{{Tl)f 



(6.40) 





6.6 Discussion 



We have calculated the third cumulant of the distribution of the total transmission and 
compared it to the data of Ref. [105]. The second and third cumulant are a conse- 
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Figure 6.3: The third cumulant plotted against the second cumulant. The points are the cor- 
rected experimental data (see text). The line is the theoretical prediction, {{Tj^)) = 3.2{{T^)y; 
we introduced no free parameters. 

quence of interference between diffuse channels. The third cumulant is proportional to 
the second cumulant squared. We also found a non-trivial dependence on the profile 
of the incoming beam used. We calculated the cumulants for arbitrary beam diameter, 
but the infiuence of a finite focus on the ratio is rather weak. Also boundary refiections 
were included. Our calculations confirm that the main contributions come from dia- 
grams with interference processes, i.e. connected diagrams, as we have shown that the 
contributions from disconnected diagrams is small. The experimentally found ratio of 
the third cumulant versus the second cumulant squared is well described by our theory. 

The extension of the calculations to higher cumulants is straightforward. The ra- 
th cumulant will contain (n — 1) Hikami four-point vertices. So the contribution is 
((T")) oc g^~". Also corrections and cancellations from higher order vertices are present, 
but it is clear that the calculation becomes very laborious at large n. We follow another 
approach to this problem in the next chapter. 
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Table 6.1: Sample thickness (in microns), beam width (in microns), second cumulant, third 
cumulant for the different samples as taken from Ref.[105]. Next are the beam diameter 
correction factor on the third cumulant (see section 6.5), and the number of modes N' and 
N*. Together they yield the corrected cumulants (plotted, not listed). 



Sample 


Beam 


Second 


Third 


Beam 


Number 


Number 


thickn. 


diam. 


cumulant 


cumulant 


corr. 


of modes 


of modes 


L 


Po 


(xlO-4) 


(xlO-^) 


factor 


N' 


N* 


30 


77 


0.36±0.01 


0.014 ± 0.035 


1.03 


388000 


300000 


12 


26 


0.97 ±0.03 


-0.03 ±0.25 


1.04 


46800 


36800 


22 


32 


1.24 ±0.04 


0.68 ±0.28 


1.06 


88600 


72100 


30 


33 


1.57 ±0.04 


1.30 ±0.46 


1.07 


119000 


99300 


53 


35 


1.80 ±0.03 


0.91 ±0.53 


1.10 


241000 


212000 


30 


26 


1.90 ±0.03 


0.92 ±0.56 


1.09 


94900 


81300 


45 


33 


1.90 ±0.05 


1.33 ±0.43 


1.10 


187000 


163000 


53 


26 


2.18 ±0.03 


1.77 ±0.59 


1.10 


208000 


189000 


170 


27 


2.69 ±0.06 


2.02 ±0.82 


1.13 


1420000 


1430000 


78 


28 


2.74 ±0.03 


2.43 ±0.62 


1.11 


396000 


372000 


30 


17 


4.82 ±0.10 


9.1 ±3.3 


1.11 


71900 


64400 


30 


10 


8.01 ±0.36 


5.3 ±6.4 


1.11 


60200 


56400 



7 

Full distribution functions 

In this chapter we calculate the full distribution of the total transmission and the an- 
gular resolved transmission. Both can be mapped on the eigenvalue distribution of 
transmission matrix. This distribution is known for moderate values of g. 

7.1 Classics 

Already from the autocorrelation functions we know that interference is important at 
moderate values of the dimensionless conductance g. For instance the C2 correlation 
is proportional to Ijg and was clearly seen in experiments. But let us first discuss the 
'classical' situation, which is the limit of large g. Then correlations between intensities 
are absent, and we recover well known distribution functions. For the angular resolved 
transmission, or speckle intensity, there is the Raleigh law [106] 

P(r„fc) = -l-e-^"^/<^"^) (7.1) 

\J-ah) 

Anyone who ever saw a laser speckle pattern on the wall, will remember the wild pattern 
with both very bright but also dark spots. This is exactly due to the Rayleigh law. 

The Rayleigh law can also be derived in the context of the ladder diagrams [92]. 
The n— th moment of the speckle intensities is made up of n amplitudes and n complex 
conjugated amplitudes. As stated before, there is no constriction on the pairing of 
two amplitudes into a diffuson for this type of measurement. Therefore the amplitudes 
{a^, ttj, 0(3, . . . , a*} have n\ possibilities to pair with the amplitudes {ai, 02, 0(3, . . . , 
Thus 

(T;,) = n!(r„,)" (7.2) 

which is easily checked to correspond with the Rayleigh law. Yet, as we know from the 
example of a laser beam refiecting on the wall, the Rayleigh is more general. We need 
no multiple scattering at all. Indeed, Rayleigh's original derivation goes as follows [106]: 
The field ?/> on a given position on the outgoing side is the sum of many fields, and 
therefore real and imaginary part have each a independent Gaussian distribution. The 
intensity / being the amplitude squared has thus an exponential distribution 

p _ g-[(Re^)H(W)^]/2<T2 _ g-//(/> 

Recently Kogan and Kaveh addressed the question of the speckle statistics of (opti- 
cally) very thin samples [121]. In very thin samples a considerable part of the light does 

95 



(7.3) 
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not scatter at all, but is coherent. In the limit of zero thickness, only coherent light, 
the distribution becomes a delta-function. In the intermediate regime a simple counting 
argument yields the distribution, which intra-polates between the delta distribution and 
the Rayleigh law. 

In the mesoscopic regime interference modifies the speckle distribution. The leading 
correction was derived by Shnerb and Kaveh [122]. Genack and Garcia observed a 
deviation from the Rayleigh law at large intensities [95, 123]. A crossover to stretched 
exponential behavior was derived by Kogan et al. [118]. 

For the total transmission and the conductance the distributions are similar in the 
large g limit. For both quantities the outgoing diffusons should be frequency and mo- 
mentum independent as we explained in chapter 4. As there is no interaction between 
the diffusons in this limit, the same applies for the incoming diffusons. Thus for both 
quantities only one type of pairing is possible. This yields in principle a (^-distribution 
if one probes an infinite number of channels. In practice of course only limited num- 
ber of channels is probed, and the law of large numbers predicts a narrow Gaussian 
distribution. 



7.2 Eigenvalues of the transmission matrix 

In mesoscopic systems the observables are random quantities and are therefore not al- 
ways characterized by the mean values, but their entire distribution functions are of 
interest (see [124] for a toy model). This is particularly prominent in the distribution of 
eigenvalues of the transmission matrix. Assuming that all eigenvalues contribute equally 
to the conductance, one expects a Gaussian distribution. As there are N eigenmodes, 
the eigenvalue distribution should be peaked gjN = j-. But this picture proves wrong. 
The eigenvalues have a "bimodal" distribution peaked around and 1. This was first 
put forward by Dorokhov [125] and later also by Imry [126]. The eigenvalues T„ of the 
transmission matrix tH can be expressed as r„ = 1/ cosh^(L/(^„). Not the eigenval- 
ues but the inverse localization lengths l/^„ are uniform distributed, see Pendry et al. 
[127] and Stone et al. [128]. In all derivations quasi-lD was assumed. Recently, how- 
ever Nazarov [129] showed that it is not only true in quasi-lD, but under very general 
conditions. It implies that: 

Or 

iT.Ui'y) = (T.Ti)=,j;^^T'. (7.5, 

This result was first derived using random matrix techniques by Mello et al. [130]. 
We plotted the distribution in Fig. 7.1. Note that few eigenvalues have the mean value 
i/ L <^ 1. Note that, the normalization of the distribution is ill-defined, the distribution 
should be understood in the sense that all its moments are well defined. This problem 
can be avoided by taking the proper lower boundary of the integral instead of 0. The 
eigenvalues are expressed as 

T„ = cosh-'{L/Q (7.6) 
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Figure 7.1: The bimodal eigenvalue distribution. The average value has only a small proba- 
bility, instead the eigenvalues are almost all either or 1. 



The channel localization length ^„ ranges (roughly) from i to infinity. The minimal 
value corresponds to the decay of unscattered intensity (whereas the upper limit cor- 
responds to the case where there is no scattering). Therefore the minimal value for T„ 
is Tinin = cosh~^(L/£) exp{—2L/i) and the distribution Eq. (7.4) can be normalized. 
Although this cut-off is important for the normalization of the distribution, its influence 
on the momenta is very small and we refrain from the correction. The moments of the 
distribution are 



1.3. 



2n 



Such that (r„) = f , (r„2) = 2/3f , and (T^) = 8/15f . The distribution Eq.(7.5) is vahd 
only in the regime where g is not too small. Loop effects near the Anderson transition, 
where (/ ~ 1, will change it, as we will discuss below. 

The physical meaning of this curious distribution function, is the concept of open and 
closed channels. An eigenmode of the transmission matrix is, according to this distri- 
bution, either essentially blocked or, with a smaller probability, essentially conducting. 
This picture was conflrmed in various computer simulations [127, 131]. Instructive is 
the simulation of Oakeshott and McKinnon [132], where the conducting channels are 
made visible in real-space; it can then be seen that the current is carried by only a few 
channels. Even in the conducting regime most channels are localized. When approach- 
ing localization the localization length of more and more channels becomes smaller than 
the sample size, untill flnally the last channel is localized, and the full sample becomes 
insulating. It also provides us with a very nice picture explaining why the correlations in 
mesoscopic systems are so large. If all the channels are conducting equally, fluctuations 
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of the channels average out by the law of large numbers. Yet if only a few channels 
are conducting, fluctuations in one channel will be clearly seen. Yet this reasoning is 
not always right, consider for example the case of the UCF. In the old picture that all 
channels contribute equally to the conductance, the variance in the conductance Ni/ L 
is caused by the flniteness of the number of channels N. In the new picture of either 
closed or open channels the number of open channels A^gff is much smaller 

N£ 

iVeff X 1 = — « TV (7.8) 

The conductance now fluctuates according the law of large numbers with a variance: 
(Y^nTn) = which is incorrect as we know from chapter 5. Below, we will see why 
the argument went wrong. 

In practice it is very difficult, if not impossible, to measure the eigenmodes and 
eigenvalues directly. As they are the eigenmodes of the very large random transmis- 
sion matrix, they have a very complex structure. Nevertheless it has some important 
observable consequences. First it was shown by Beenakker and Biittiker [133] that the 
shot-noise of electronic conductors universally reduces by a factor 3 because of this dis- 
tribution. Shot-noise only occurs in electronic systems and is thus not of much interest 
for us. Secondly, we will show the total transmission and speckle intensity distribution 
function are fully determined by the eigenvalue distribution function. 



7.3 Distribution of total transmission 

We first calculate the probability distribution of the total transmission. Already from 
the distributions without interference effects, see section 7.1, we saw that the total 
transmission has a simpler distribution (in the sense of cumulants) than the speckle. 

We neglect all corrections for absorption, skin layers, and disconnected diagrams. 
Consider again an incoming plane wave in direction /i^ (/U denotes again the cosine 
with respect to the 2;— axis). The wave is transmitted into outgoing channel h with 
transmission amplitude tab = 2k^JJi^^Gah and transmission probability Tab = \tah\^ ■ 
The average total transmission is obtained by summing all outgoing channels, see chapter 
2 

{T.) = {Y.Tab) = "-^. (7.9) 
We find the average conductance by 

9^UTa) = ^, (7.10) 

thus (Ta) = eaQ^ while one also has {Tab) = ^a^bQ- 

We consider the j^^ cumulant of T^. In a diagrammatic approach this object has 
j transmission amplitudes tab and an equal number of Hermitian conjugates = t*^. 
The explicit calculation for the second cumulant C2 (chapter 4), and third cumulant 
(chapter 6) showed that the leading diagrams are connected, yet have no loops. 

Let us fix the external diffusons in the term t^i, tl A, ■ ■ -t^i^.tl.^. Contributions to 

aoi via a02 o,0j oja 

the sum over bi only come from diagrams with outgoing diffusons that have no transversal 
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Figure 7.2: Schematic picture of tlie leading diagrams for the fourth cumulant total transmis- 
sion. Although there are 6 diagrams, there are only three topologies on intensity level. On the 
left the thick lines represent tab, the thin lines ij^^; this type of diagram simplifies the counting. 
On the right the more conventional picture using Hikami boxes. 



momentum. These are the diagrams where the hues with equal bi pair into diffusons. 
We indicate this pairing of the outgoing diffusons with a underlining bracket 

The outgoing diffusons are fixed now. For the incoming side it is convenient to keep 
for instance the t^fe's fixed, while the tj^^'s are permuted among them, see Fig. 7.2. This 
leaves us j\ possibilities for the incoming side. The number of Hikami boxes in this 
process "misfits". The number of misfits for diagrams that are connected and loop-less, 
is J — 1, giving (j — 1)! possibilities. Diagrams with less misfits are disconnected, diagrams 
with more contain more vertices and are thus sub-leading. (Of course just as happens 
for the third cumulant there will be local corrections from, for instance the six point 
H(y vertex, we do not consider them explicitly here). Examples of allowed diagrams are 
cyclic permutations of t),^ with respect to the outgoing side. 

Next we factor out the incoming and outgoing diffusons and group the remainder of 
the diagrams into a skeleton K . Making use of Eq. (2.68) we obtain: 

(ri)co„ = e^„(j-l)!|dridr;---dr,dr;Xi„(ri)£out(r;) 

• • •£i„(r,)£o„t(r;)K(ri,r;, • • • ,r„r;.). (7.12) 

The integral just describes 

(Tr(tt^)-') = Tr ^ (^aifci4ia2^a2fc2 ' ' ' ^ajfcj^jai ) 

6l ,a2,b2,---,aj,bj 

= X] (^aifci4ia2^a2fc2 ■ ■ ■ ^a,fcj4,ai )• (7-13) 



There is only one way to attach incoming and outgoing diffusons to K. The sums over 
the indices lead exactly to the total-flux diffusons in Eq. (7.12). We thus flnd 

(ri)co„ = (j - l)!4(Tr(ttt)^), (7.14) 
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which is the crucial step in the derivation. Normalizing with respect to the average, we 
introduce Sa = Taj {To). The generating function of the connected diagrams is easily 
calculated 

J = l 

= glog^ ijl + xlg + ^j^g^ . (7.15) 

Since the cumulants are solely given by connected diagrams, the distribution of Sa follows 

as 

— eXp[xSa-^conix)]. (7-16) 

-ioo Zm 



In Fig. 7.3 we present the distribution (7.16) for some values of g. At moderate g 
we clearly see the deviation from a Gaussian. Let us examine some properties of the 
distribution. For Sa near unity and large g we can expand $ up to order x^, recovering 
the Gaussian behavior found by Kogan et al. [118]: 

P(3„)^^/^exp[-^(3„-l)^] (7.17) 

This corresponds to including only the C2 contribution. The integrand in Eq. (7.16) has 
a branch cut from x = —g to x = — 00. For < the contour can be closed to the right 
half-plane and P{sa) vanishes. For limiting values of Sa we use a saddle point analysis. 
The saddle point is found by the condition ^ [xSa — ^con{x)] = 0. Thus we find 



log i^^l + x/g + ^x/g) 

■^a — I — — — 

yx/g{l + x/g) 

The r.h.s. of Eq.(7.18) diverges if x approaches —g (from above) and decreases mono- 
tonically for larger x. Thus for large 1, we find the saddle point near x = —g. By 

inserting the saddle point one finds a simple exponential decay 

P{s,) ^ exp(-^s„ + g—), Sa > 1. (7.19) 

The saddle point also dominates the shape for small Sa <^ 1 (and large g). This time 
x/g has to be large to solve Eq.(7.18), i.e. Sa ~ 9(log4a;/(/)/2a;. One finds essentially a 
log-normal growth: 



P{sa) ~ exp 



g g ( z z ^ ^ 

4-4rsi:+^°s^°si:-' 



(7.20) 
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Figure 7.3: Intensity distribution of the total transmission, in units of ^3^/47r, versus the 
normalized intensity Sa for an incoming plane wave. With g = 2,4, 16, and 64 (upper to lower 
curve). 

7.3.1 Influence of beam profile 

So far we have considered the case of an incoming plane wave. We already mentioned 
that in optical systems a Gaussian intensity profile is more realistic. For the third 
cumulant we saw a non-trivial dependence on the incoming beam profile, suggesting 
that also higher cumulants are sensitive. For perpendicular incidence the incoming 
amplitude is V'in(i") = X)a </'(9a)V'in('")5 where V'in('") plane wave of Eq. (2.48), 

and where 



We consider the limit where the beam is much broader than the sample thickness [po ^ 
L) but still much smaller than the transversal size of the slab [po <^ W). If we would 
consider a smaller beam diameter, the incoming transverse momenta, which are of order 
1/po, become of the order of 1/L. The diffusons will then become the well known cosh- 
functions, chapter 2. Here the momentum dependence of the diffusons can be neglected; 
apart from a geometrical factor, the situation is identical with the plane wave case. Due 
to integration over the center of gravity, each diagram involves a factor A^Sg.Sg'- In the 
j^^ order term there occurs a factor 





A 



gig'i---gjg'j 





For a plane wave 
obtain 




For our Gaussian beam we 
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Figure 7.4: Intensity distribution of the total transmission, in units of ^J'igjl'K^ versus the 
normalized intensity for an incoming wave with Gaussian profile. With ^ = 2,4, 16, and 64 
(upper to lower curve). 



It is thus convenient to identify Aq = \t^Po with the effective area of a Gaussian beam. 
(This definition is different from the one used in previous chapters, were Aq = t^Pq)- 
As compared to the plane wave case, the j^^ order term is smaller by a factor 1/j for a 
Gaussian profile. This implies for the generating function of the connected diagrams 



$co„(x)=W ^logM, 1 + ^ + . (7.24) 




For small Sa (and large g) there is again a log-normal saddle point. For large Sa the 
dominant shape of the decay is given by the singularity at x = —g and again yields 
P{-Sa) ~ exp( — (/Sa)- In Fig- 7.4 we present the distribution function for different values 
of g. The shape of the distribution is quite similar to the plane case. A difference in the 
figures is that for a Gaussian beam profile the effects seem smaller. One could say that 
the Gaussian profile reduces the fiuctuations, but it is mainly a question of definition of 
the number of modes N and the area of a Gaussian intensity profile. In the graphs the 
distribution for the Gaussian profile at g almost overlaps the plane wave distribution 
for 2 X g. Although the analytic forms are different, the difference is small. 

By expansion of the generating function we recover the results for the second and 
third cumulant, obtained previously 

(■5a)cum = ~, (■5a)cum = ~^{^a)cum (7.25) 

for a plane wave. For a Gaussian beam we recover 

3g' ^ 5 

The results coincide with the results of chapter 6 and describe the experiments of Ref. 
[105] well. 



■5a)cum — „ 7 (■5a)cum — _ {^a)cum- (7.26) 
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7.4 Speckle-intensity distribution 

We apply the same mapping to obtain the distribution of the angular transmission 
coefficient. The angular and total transmission can be related to each in fairly simple 
way, because the interference processes are dominated by the same type of diagrams: the 
loop-less connected diagrams. This is in contrast to the conductance that is dominated 
by far more complicated loop diagrams. The different distribution for total transmission 
and speckle is the consequence of a counting argument only. In the plane wave situation 
the average angular transmission reads {Tab) = ^a^hQ- Let us count the number of 
connected loop-less diagrams that contribute to T^^ = tahA,a.^ah' ' '^ab^la- Now not only 
massless outgoing diffusons, but all pairings into outgoing diffusons contribute. This 
yields an extra combinatorial factor j\ in the j^^ moment: 

{TL)con = jKj-i)K4{^4tUy)- (7.27) 

For the normalized angular transmission coefficient Sab = Tab I {Tab) , we introduce the 
following generating function of the connected diagrams: 



*co„(:^) ^ E ^ ^\ (^ab)con (7.28) 

J = l 



It is easy to see that 

*con(a:) = $con(a:), (7.29) 

with $con given by Eq. (7.15) for plane wave incidence and by Eq. (7.24) for a broad 
Gaussian beam, respectively. In contrast to the total-transmission distribution, the 
cumulants are not only given by the connected diagrams. As example one can look again 
at Fig. 6.1; all the diagrams there contribute to the third moment of both speckle and 
total transmission, but for the speckle pattern there are 3! diagrams for each topology. 
Kogan et al. [118] showed that the summation of the disconnected diagrams can be done 
elegantly by performing an additional integral: Consider again the jth moment (s^^). 
First there is the factor j\ from the different pairings of outgoing amplitudes. Consider 
a disconnected diagram made up of several connected ones. The number of connected 
parts with i incoming (and i outgoing) diffusons is labeled m,-. For instance Fig. 6.2 is 
characterized with nii = 1 and m2 = 1, as it contains one non-interacting diffuson and 
a part connecting 2 diffusons. The total number of diffusons is fixed so m obeys: 

j2^m,=j (7.30) 

i=l 

The number of different topologies that can be obtained is 

n\ 



- (7.31) 
mi\m2\ . . . m„l 

But now we have over-counted as permuting the diffusons of a connected part yields no 
different contribution. Therefore we have to divide with (z!)™' for each i. Summarizing 
we thus find for the moments of P[sab) 

Wab) = j! E Pimu m2, . . . , m„){Sab)Z'n{slb)Zl ■ ■ ■ Wab^^n (7-32) 
mi ,1712, ...,mj 
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Here the primed sum means that Eq.(7.30) should be fulfilled. The P is the combina- 
torial prefactor 

P(mi, m2, . . . , m„) = — — — — — ; ; (7.33) 

^ ' ' ' ^ (l!)™i(2!)™2 ...(j!)™.>„J^2!---m,! ^ ^ 

This factor is one for the connected diagram [nij = 1). Now we know the moments the 
distribution follows from them as 

P{sa,) = r ^-^{si,) exp(x3„,) (7.34) 

It is useful to write the factorial in Sab in Eq.(7.32) as an integral 

j\ = / dvv-' ex-p{—v) (7.35) 

JO 

Yielding 



P{sab) = n ^e^s^bxlv) r -exp{-v)J2 -4^M, (7.36) 

J-ioo ZTTl Jo V TTn ?! 



dx , , , [^Av , (-xY 

j=o 

With M, = Also 



h ' " ™.! I 1! ) V 2! ) ■■■ 

= exp(E^(3^„,)_) (7.37) 
j=o 3 ■ 

Or in words, the sum over the disconnected diagrams replaces the exponent of the sum 
of connected diagrams. One gets 

P{Sab) = / — exp \-XV-Wc^on[X 

Jo V J-ioo ZTTl \ V 

dx 



/too (j^j- 
— Ko(2y^:^)exp(-*con(a:)). (7.38) 
-ioo TTl 



With Kq the modified Bessel function. We plotted the speckle intensity distribution in 
Fig. 7.5 for an incoming plane wave. For large g and moderate Sab we have ^con{x) ~ x 
and we recover the Rayleigh law: P[sab) = exp(— Safc)- The leading correction is found 
by expanding m 1/g 



P(sab) = e" 



■Sob 



l + ^(5afc-45„fc + 2) 



(7.39) 



This was derived previously by Shnerb and Kaveh [122]; here we have related the pref- 
actor of the correction term to the conductivity. Genack and Garcia fitted their data 
to this relation and found g = 14.6. Our Eq. (7.38) describes these data very well for 
g = UA. 
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Figure 7.5: Intensity distribution of speckles versus the normalized intensity Sab for an in- 
coming plane wave. With g = 2,4, and 8 (upper to lower curve at Sab = 5). The dashed line 
corresponds to the Rayleigh law {g = oo). 



For large Sab one can again apply steepest descent, which yields 



by Kogan et al. Their findings are based on truncating ^con{x) after order x"^ . This 
corresponds to including only the simplest connected diagram: the C2 diagram. With 
this diagram also all corresponding disconnected diagrams are constructed. Although 
their theory works quite well, one introduces a extra free parameter for each higher order 
correction, whereas we find the distribution as a function of a single parameter, i.e. g, 
only. Taking the full generating function into account, we find a qualitatively different 
saddle point. The stretched exponential was also seen in experiments. Already in 1989 
Garcia and Genack observed it in their microwave experiments. But unfortunately their 
dynamic range was rather small. The fit of stretched exponential is therefore rather 
unprecise. The maximum intensity was 5 times the average, here one does not expect 
to see the tail already. Rather one observes a cross-over behavior here. However, more 
accurate measurements by Genack and Polkosnik are planned. Nevertheless, in the total 
transmission experiments, which were the inspiration for this research, the higher order 
corrections are easier to extract from the data. 

Using Eq. (7.24) we derive the speckle distribution due to an incoming beam with 
Gaussian profile. A Gaussian profile leads to a different distribution with the same 
asymptotic behavior; it is plotted in Fig. 7.6. 




(7.40) 



This stretched exponential tail differs from the form P{sab 




SabYg] asserted 
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Figure 7.6: Intensity distribution of speckles versus the normalized intensity Sab for an in- 
coming wave with Gaussian profile. With g = 2,4, and 8 (upper to lower curve at Sab = 5). 

7.4.1 Joint distribution 

It turns out that the speckle distribution and the total transmission distribution of a 
certain incoming direction are related.^ One can expect this, a large total transmission 
for a given incoming direction, will also be reflected in the individual speckles. The 
moments of the joint distribution are 

{■4b-i)con = {k + l-l)lk\{T!:+') (7.41) 

The combinatorial factor k\ is for the speckle pattern, the factor [k -\- I — 1)! is the 
number of possible pairings of incoming beam. Deflning a as Sab = CafcSa, one obtains 
(^afcO = k\{m- 1)! (T™) with m = k + L Thus P{a,b,s,) = exp(-a„fc)P(5„). Thus 
one sees that (Tab and Sa are independent variables, yet Sab and Sa are dependent. Their 
joint distribution is 

P{Sab, Sa) = ^ / — eXp{xSa - $co„ • 7.42 

Sa J ZlTl 

We have plotted this distribution in Fig. 7.7. Integration over Sa or Sab yields again 
the distribution P(safc), or P[sa), respectively. As the r.h.s in Eq. (7.42) contains only 
connected diagrams, the same holds thus for the l.h.s. Nevertheless, by analytic continu- 
ation for negative / and the subsequent integration over Sa we flnd the correct expression 

for P{Sab). 

7.5 Concluding remarks 

In the above, we derived the full distribution function of both total transmission and 
angular transmission coefficient by mapping it to the distribution of eigenvalues of Tab- 

^This calculation is an unpublished result of Th. M. Nieuwenhuizen 
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Figure 7.7: The joint distribution function of Sa and Sab, for a Gaussian beam profile, for 
g = 2. The maximum as a function of Sa shifts as Sab changes. 

But at least for the lower one we can use this mapping the other way around. First, 
the experiments seem to confirm the validity of the eigenvalue distribution. For a long 
time the prove for this eigenvalue distribution was done by the sometimes error-ing ran- 
dom matrix theory for quasi ID systems. Only recently Nazarov showed its validity in 
arbitrary dimension. Probably the mesoscopic systems using multiple scattered electro- 
magnetic waves, i.e. microwaves or light, are still the cleanest to prove this fundamental 
distribution. 

Secondly, from our detailed evaluation of the low order diagrams, i.e. the second 
and third cumulant, insight in the eigenvalue distribution can be obtained. It is clear 
from those calculations that also the eigenvalue distribution is based on loop-less con- 
nected diagrams. Therefore the eigenvalue distribution cannot directly be used for the 
conductance. The explicit evaluation of second and third cumulant made it also clear, 
that the results are rather insensitive to absorption and surface effects. 

Note that mapping from the eigenvalue distribution to the speckle and total trans- 
mission distribution is very general. Most interesting would be the inclusion of loops 
into the eigenvalue distribution, which would clarify the situation for small values of g 
but also for the tails. In the important paper of Altshuler, Kravtsov and Lerner, [88] 
it was seen that although loop diagrams may be of lower order, they do become impor- 
tant at high moments (see below for an explanation). Therefore, our results for the tail 
change maybe. Explicit inclusion of the loops would solve this question. 

7.6 Higher moments of the conductance 

It would be interesting to apply the technique above also to the conductance. The 
distribution function of the conductance has been studied extensively, mainly with the 
advanced models mentioned in the introduction, chapter 1. For a review see Altshuler, 
Kravstov, and Lerner in Ref. [23]. Using the non-linear sigma model they find for the 
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lower cumulants (n < g) the estimate 

{g-Um-^g'-". (7.43) 

But at large n there is a topological growth of the diagrams: Introducing an extra 
loop in a diagram does introduce a factor 1/g^ making the diagram smaller, but on the 
same the number of topologically different diagrams increases with roughly a factor 
for large n. It is calculated that the tail of the distribution becomes essentially log- 
normal [134]. Their calculation seems to violate the one parameter scaling theory [12], 
as for the higher moments not only the dimensionless conductance (/, but also extra 
parameters are necessary (however, see also Shapiro [135]). The question whether this 
is really the case is not resolved yet. Surprisingly enough, despite the vast number of 
measurements of the UCF, i.e. ((/^), there are hardly any experimental data on the full 
distribution of g. The question of the precise distribution function is still open, as the 
pre-factors in Eq. (7.43) are unknown. Recently Macedo showed in the quasi ID case, 
that {g^)cum cx (/~^, which would mean that some "miraculous" cancelation occurs for 
the term proportional to 1/g [136]. 

Let us return to our diagrammatic approach. We need to construct connected dia- 
grams with 2n amplitudes on the incoming side going to the outgoing side. The paring 
into diffusons should be the same on incoming and outgoing side. From Eq. (7.43) one 
observes that the nth relative cumulant should be proportional to g^~^". As in relative 
cumulants, each Hikami-box brings a factor (/~^, the diagrams for the nth cumulant 
should consist of 2n — 2 Hikami-boxes. Reconsider the UCF diagrams of chapter 5. 
A possible method for generating these diagrams is: take the diagram for the third 
cumulant (chapter 6) and join an incoming to an outgoing diffuson, thus creating one 
extra internal diffuson. The joining of the diffusons can be done using the A'-operator 
of section 2.6. Apart from forbidden diagrams with a loop of amplitudes^ a subset of 
the C3 diagrams is generated. Yet we do not find diagrams with most crossed diagrams 
internally instead of diffusons, i.e. the diagrams of lower line in Fig. 5.1, by this method. 
In two and three dimensions, the extra internal momentum integral has to be inserted. 
Unfortunately, this research is not concluded yet. 



A 

List of symbols 



In general three-dimensional vectors are denoted as lower case boldface letters, and 
two-dimensional vectors with normal face capitals. 
In the diagrams: 

- Dots like: • denote a t-matrix, except in chapter 5 where they denote a potential. 

- Diffusons are often depicted as close parallel lines 

- Dashed lines connect identical scatterers 



Symbol 


Meaning 


Introduction 


eo 


dielectric constant inside slab 


Eq. (2.4) 


ei 


dielectric constant outside slab 


Eq. (2.4) 




diffuson prefactor 


Eq. (2.67) 


K 


absorption rate 


Eq. (2.34) 


11 


cosine angle w.r.t 2;— axis 


Eq. (2.48) 


P 


transverse coordinate 


Eq. (2.47) 


Po 


beam diameter 


Eq. (4.11) 


Tl 


limit intensity 


Eq. (2.60) 


id 


frequency of light 


Eq. (2.4) 





reduced frequency difference in diffuson 


Eq. (2.36) 
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Symbol 


Meaning 


Introduction 


A 


beam area 


Eq. (2.47) 


a 


albedo 


Eq. (2.28) 


ao 


scatterer radius 


Eq. (2.6) 


Ci 


short range correlation 


Eq. (4.1) 


C2 


long range correlation 


Eq. (4.1) 


C3 


UCF correlation 


Eq. (4.1) 


G 


dressed Green's function 


Eq. (2.22) 


Go 


bare Green's function 


Eq. (2.9) 


9 


dimensionless conductance 


Eq. (6.2) 


H, 


Hikami four point vertex 


Eq. (3.5) 


He 


Hikami six point vertex 


Eq. (3.7) 


/i4, he 


vertex-prefactor 


Eq. (5.8) 


Ik,i 


loop integral 


Eq. (2.31) 


K 


extra scatterer operator 


Eq. (2.82) 


k 


wave vector 


Eq. (2.7) 


L 


sample thickness 


sec. 2.1 


C 


diffuse intensity 


sec. 2.4 


£ 


mean free path 


Eq. (2.26) 


M 


intensity-decay rate 


Eq. (2.40) 


m 


index of refraction ratio 


Eq. (2.7) 


N 


number of modes 


Eq. (2.78) 


n 


scatterer density 


Eq. (4.11) 


P 


3D microscopic momentum p ~ A; 


Eq. (2.30) 


Q 


2D transverse momentum 


Eq. (2.40) 


q 


3D mesoscopic momentum q£ ~ 1 


Eq. (2.30) 


Rill) 


boundary reflection coefficient 


Eq. (2.44) 




normalized Ta 


Eq. (7.3) 


Sab 


normalized Tab 


Eq. (7.4) 


Til,) 


boundary transmission coefficient 


Eq. (2.51) 


Tab 


angular resolved transmission 


Eq. (2.79) 


T 


total transmission 


Eq. (2.80) 


t 


t— matrix 


Eq. (2.19) 


Z 


depth in the sample 


chapter 2 




extrapolation length 


sec. 2.4, (2.41) 


{{0')) 


normalized cumulant 
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Mesoscopische verschijnselen in 
veelvoudige lichtverstrooiing 



Onze ogen zijn onze belangrijkste zintuigen en het licht dat zij zien kan veel interes- 
sante verstrooiingsprocessen ondergaan. Denk bijvoorbeeld aan een zonsondergang of 
aan de kleurschakeringen van een vlinder. Een van die verstrooiingsprocessen is uit 
esthetisch oogpunt niet bijzonder aantrekkelijk, maar voor natuurkundigen een nadere 
studie waard, namelijk licht dat vele malen verstrooit op een coUectie verstrooiers. Het 
veelvoudig verstrooide licht wordt dan diffuus. Een alledaags voorbeeld hiervan is licht 
dat zich voortplant door mist of door melk. In ons laboratorium echter worden de exper- 
imenten gedaan door kleine plakjes witte verf te beschijnen met een laser. Een laser en 
wat witte verf, dat klinkt simpel, maar nauwkeurige metingen zijn echt niet eenvoudig. 
Juist de hoge kwaliteit van de experimenten bleek van groot nut om ook de kleinste 
effecten te kunnen zien. 

Wat is nu zo interessant aan dat veelvuldig verstrooide licht? AUe interessante in- 
formatie wordt toch immers door elkaar gehusseld? Laat ik dat proberen uit te leggen: 
In heldere lucht of glas volgt het licht een rechte weg. In die witte verf wordt door 
alle verstrooiingen de gevolgde weg kronkelig en de afgelegde weg wordt erg lang. De 
zogenaamde diffusiebenadering geeft nu een redelijke beschrijving. Maar er is meer aan 
de hand. Als verstrooiing sterker wordt (de mist dikker wordt) ontstaat er een kans dat 
lichtbundels elkaar ontmoeten. Omdat licht een golfverschijnsel is, kunnen de bundels 
dan interfereren. Nu is deze kans op interferentie meestal erg klein en wordt in de dif- 
fusiebenadering dan ook verwaarloosd. Maar deze effecten zijn in de plakjes witte verf 
wel degelijk waarneembaar en voor mij zijn ze het zout in de pap. Men spreekt ook 
wel over mesoscopische effecten, aanduidende dat men kijkt naar lengteschalen ergens 
tussen enerzijds de macroscopic (de diffusie, waar alles uitmiddelt) en anderzijds de mi- 
croscopic (de individuele verstrooiingsprocessen). Overigens komen deze mesoscopische 
verschijnselen ook voor bij andere diffuse golven zoals geluidsgolven en microgolven. 
Zelfs het transport van elektronen door een stukje verontreinigd metaal valt onder deze 
beschrijving. Sinds de introductie van de quantummechanica weten we immers dat ook 
elektronen een golfkarakter hebben. De theorie in dit werk is op al die systemen binnen 
zekere grenzen toepasbaar. De optische experimenten hebben echter het voordeel zeer 
zuiver en nauwkeurig te zijn. 

Welke invloed hebben die mesoscopische effecten op het transport? We onderscheiden 
twee mogelijkheden voor de interferenties. Ten eerste is het mogelijk dat een lichtbundel, 
na de nodige omzwervingen, zichzelf weer tegenkomt en interfereert. Die terugkeerkans 
blijkt groter te zijn dan men op grond van de diffusietheorie zou verwachten. Omdat 
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licht vaker terugkeert, is de effectieve voortplantingssnelheid kleiner. Dit effect staat 
bekend als zwakke localizatie. Voor elektronen is het zelfs mogelijk de verstrooiing zo 
sterk te maken, dat transport onmogelijk wordt. Op dat moment zijn de elektronen 
niet meer mobiel en de weerstand van het metaal wordt bijzonder groot. Dit wordt ook 
wel sterke of Anderson-localizatie genoemd. Voor licht is Anderson-localizatie nog nooit 
gezien, omdat het niet eenvoudig blijkt licht zo sterk te verstrooien. Optisch wel goed 
waarneembaar is de tweede mogelijkheid voor interferentie. Dit is het process waarbij 
de ene bundel die de witte verf doorkruist, interfereert een andere bundel. Vanwege 
de mogelijke interferentie zuUen de bundels zich niet meer onafhankelijk voortplanten, 
maar verstrengeld raken. Dit komt tot uitdrukking in hun correlatiefunctie, welke de 
mate van relatie tussen de twee intensiteiten aangeeft. Over dit soort processen gaat 
deze studie. 

Het aardige van optische systemen is dat men niet snel is uitgeexperimenteerd: In 
tegenstelling tot electronische systemen, waar men eigenlijk alleen de weerstand kan 
meten, zijn er voor optische systemen drie transmissiegrootheden. AUereerst kan men 
hoekopgelost meten. Voor zo'n meting beschijnt men de witte verf met een evenwijdige 
bundel uit een zekere hoek, en meet de uitgaande intensiteit in een zekere richting. 
Dit geeft de hoekopgeloste transmissie. Daarnaast is het ook mogelijk al het uitgaande 
licht te verzamelen om zo de totale transmissie te verkrijgen. Ten slotte is het ook 
nog mogelijk met vele richtingen in te schijnen en ook weer al het uitgaande licht op 
te vangen. In dat geval meet men de geleiding. Het blijkt dat de correlatiefuncties 
van deze drie grootheden zeer verschillend van elkaar zijn, zowel wat betreft orde van 
grootte, als wat betreft eigenschappen. 

Na een inleiding, begin ik dit proefschrift met de diffusievergelijking en de stralings- 
transportvergelijking. Zij bepalen hoe de diffuse intensiteit zich gedraagt in de plak 
witte verf. Ik zal bespreken hoe het precieze gedrag nabij de rand van de plak moet 
worden meegenomen, want dat blijkt van belang te zijn in de experimenten. Ik sta ook 
stil bij de vraag of het mogelijk is, ondanks de diffusie, objecten in het verstrooiende 
medium waar te nemen. Dit sluit aan bij onderzoek in de medische fysica naar de 
mogelijkheid om dergelijke technieken tomografisch te gebruiken voor het opsporen van 
bloedstolsels of tumoren. 

In hoofdstuk 3 introduceer ik een andere belangrijk element voor de berekeningen: de 
Hikami-box. Deze beschrijft de interferentie tussen twee bundels. De berekeningsmeth- 
ode die ik gebruik is een diagrammatische techniek. Die komt erop neer dat men alle 
mogelijke verstrooiingsprocessen van een zekere klasse systematisch tekent, uitrekent en 
optelt. Het diffuse transport wordt in diverse afbeeldingen aangegeven twee parallelle 
lijnen, terwijl de Hikami-box wordt aangegeven met een grijs vierkantje. Zo, nu be- 
grijpt de lezer hopelijk die diagrammetjes, parallelle lijnen: simpel diffuus transport, 
vierkantjes: een interessant interferentieproces. 

Nu ik alle ingredienten heb, begin ik met de berekening van de mesoscopische ef- 
fecten. De correlatiefuncties van de drie verschillende transmissiegrootheden worden in 
de hoofdstukken 4 en 5 uitgerekend. Ik bereken de correlatiefunctie van twee bundels 
die beide de witte verf doorkruisen, maar een net iets andere invalshoek of een net iets 
andere kleur hebben. Waar mogelijk vergelijk ik onze uitkomsten met experimentele 
gegevens. Hoewel de diagrammatische techniek het voordeel heeft recht-door-zee te zijn 
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en uiteindelijk altijd de juiste uitkomst levert, komt in hoofdstuk 5 een nadeel van deze 
techniek naar voren. Het is namelijk soms moeilijk alle relevante processen uit te vissen 
en hun aantal kan erg groot worden. De berekening wordt nogal arbeidsintensief. Ik 
laat zien dat, ondanks scepsis in de literatuur, deze harde noot toch kraakbaar is. 

In hoofdstuk 6 beschouw ik de samenhang tussen drie bundels voor het geval men de 
totale transmissie meet. Voor dit effect moeten dus drie bundels met elkaar interfereren, 
hetgeen niet zo waarschijnlijk lijkt. Maar hoewel dit inderdaad een erg klein effect 
blijkt, werd het onlangs experimented gemeten. Ik vind goede overeenstemming met de 
experimenten. 

Door deze berekeningen uit te breiden naar nog meer bundels is het mogelijk de 
voUedige verdeling van de uitgaande intensiteiten uit te rekenen voor hoekopgeloste en 
totale transmissie metingen. Dit gebeurt in het zevende en laatste hoofdstuk. Ook 
hiermee zijn experimentele gegevens, dit maal van microgolfexperimenten, in overeen- 
stemming. Het blijkt dat de fluctuaties rond het gemiddelde veel groter zijn dan men 
op grond van de diffusiebenadering zou verwachten. Dit is ook wel te begrijpen: in de 
diffusiebenadering zijn de lichtbundels onafhankelijk en dus zuUen variaties in de ene 
bundel al snel door een andere bundel uitgemiddelt worden. Maar in werkelijkheid zijn, 
ten gevolge van al die interferentie processen, de bundels niet onafhankelijk. Door dit 
gemeenschappelijke gedrag zal er veel minder middeling plaatsvinden en zuUen dus de 
fluctuaties groot zijn. 

In het laatste hoofdstuk worden de conclusies van dit werk het duidelijkst: 1) In- 
terferentie processen leiden tot interessante effecten bij veelvoudig verstrooide golven. 
Er blijken namelijk sterke correlaties en grote fluctuaties in de doorgelaten bundels te 
ontstaan. 2) Optische systemen zijn bijzonder geschikt om dit waar te nemen. 
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